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ABBREVIATIONS 


ATM Atmospheric Pressure 

BPMSL British Petroleum Mean Sea Level 

CFDD Cumulative Freezing Degree-Day 

COVID Coronavirus 

CRIB Colville River Ice Bridge 

CWAT Community Winter Access Trail (approximately 300 miles of groomed snow 
roads connecting Utgiagvik, Atqasuk, Wainwright, and Nuiqsut, to gravel 
infrastructure) 

DS Downstream 

ft Feet 

ft? Square Feet 

ft?/sec Cubic Feet Per Second 

GMT2 Greater Moose’s Tooth Drill Site 2 

GPM Gallons Per Minute 

IE Innovative Civil Engineering, Design, and Consult 

MBI Michael Baker International 

MT7 Moose’s Tooth 7 (occasionally used interchangeably with GMT2) 

NAD83 North American Datum 1983 

OPIB Ocean Point Ice Bridge 

ak Pressure Transducer 

QO Discharge 

SPMT Self-Propelled Modular Transporter 
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STA Survey Station (feet) 
US Upstream 
WSE Water Surface Elevation 
DEFINITIONS 
Albedo A ratio of the reflected insolation to the incident insolation for any of the 
natural or manmade materials 
Insolation Exposure to the sun’s rays 
Overflow When water is present on the surface of the ice 
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1 EXECUTIVE SUMMARY 


An investigation of the proposed Ocean Point Ice Bridge (OPIB) was conducted in the winter 
of 2019-2020. An updated ice bridge design is included in Appendix E that is suitable to 
support the Willow Development maximum module net weight of 3,200 tons loaded on a 
Self-Propelled Modular Transport (SPMT) with maximum allowable gross weight of 
4,200 tons. The construction quantities are presented in Table 1.1. 


TABLE 1.1: ICE CONSTRUCTION QUANTITIES OCEAN POINT ICE BRIDGE 
Ocean Point Ice Bridge Ice Quantity (yd°) 
**TOTAL = 36,500 


“The water equivalent quantity is 6.6 million gallons 


Observations from the winter of 2019-2020 have confirmed that there is a significant potential 
for overflow at the OPIB. Furthermore, there is a decrease in water discharge throughout the 
course of the winter. These two factors were the reason for the revised OPIB design. 


Three potential overflow events were documented at Ocean Point in the winter of 2019-2020 
during January, February, and March. The overflow event in March temporarily closed the 
Community Winter Access Trail (CWAT). The average depth of the overflow was estimated 
to be less than 1 foot. 


Two direct discharge measurements were conducted at Ocean Point by Michael Baker 
International (MBI). The first and discharge, 135 ft3/s (60,000 GPM) was measured on 
December 31, 2019. This had decreased to 9 ft3/s (4,000 GPM) by February 25, 2020 
(93% decrease). A third direct discharge measurement was planned for April 14, 2020 but it 
was cancelled due to COVID travel related protective measures. The discharge probably 
continued to decrease through April and prior to runoff. 


Future data collection efforts should be conducted at Ocean Point to better understand how 
the water discharge varies throughout the winter. It is important to verify the frequency and 


“magnitude of overflow events and the mechanisms that create them. Investigations should 


focus on the period between mid-February to mid-April. This is the planned timeframe for 
construction and use of the OPIB (2024-2025). 
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2 INTRODUCTION 


During the winter of 2018-2019 data were collected at 6 locations along an 8.2-mile reach of 
the Colville River in the vicinity of Ocean Point. The primary objective of this investigation 
was to identify potential locations across the Colville River where a heavy haul ice bridge 
could be constructed with a load capacity greater than the Colville River Ice Bridge (CRIB), 
and unlike the CRIB can be traversed by Self-Propelled Modular Transports (SPMT). Two 
potential locations were identified that met these criteria: 


1. Transect #1 - Henceforth known as the Upstream Site 
2. Transect #5 /#6 - Henceforth known as the Downstream Site 


Figure 2.1 shows the approximate positions of the two locations overlaid on the ‘Willow 
Optimization Option 3’ that is presented in Appendix A. A follow-on study was planned for 
the winter of 2019-2020. The Upstream Site was selected for study because of its location on 
existing permitted routes and its ease of access. 


The objectives for the 2019-2020 investigation were to determine the ice and hydrological 
conditions of the Upstream Site. The Downstream Site is assumed to have similar 
hydrological, bathymetric, and topographic characteristics as the Upstream Site - the details 
of which are contained within this report. 


UPSTREAM 


FIGURE 2.1: PROPOSED OCEAN POINT ICE BRIDGE CROSSINGS 
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2.1 PURPOSE OF 2019-2020 OCEAN POINT RECONNAISSANCE 


The purpose of the 2019-2020 effort was to collect information about the Upstream Site for the 
Ocean Point Ice Bridge (OPIB). The data would be used to document and better understand 
how the crossing evolves throughout the course of the winter. Additionally, the new 
information has been used to revise the ice bridge design, construction quantities. To 
determine the maximum load capacity of the ice bridge the following data were collected: 


1. Cumulative Freezing Degree-Days (CFDD) 5. Water surface elevation (WSE) 
- 2. Natural ice growth 6. Water Discharge (Q) (MBI 2020) 
- 3. Span of free water 7. Overflow observations 
- 4. Local weather observations 


The new information is presented in Section 4 and 5. 
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3 SUPPOSITIONS 


The following suppositions should be considered while reviewing this report. 


1 


O1 H W bh 


. The first year of moving SPMTs with modules will likely be during the winter of 


2024-2025. 


. The module move time span will be mid-February to mid-April. 

. The ice structures are subject to revision based on new information. 

. There is a potential for overflow at the proposed OPIB site. 

. This reports focus is limited to the Upstream Site for the OPIB. The Downstream Site 


is assumed to have similar physical characteristics (ice growth, water depths, 
discharge, construction volumes, topography, bathymetry, geometry, etc.). 


. The use of culverts is not addressed in this report and their potential usefulness has 


been published in the Innovative Civil Engineering, Design, and Consult (ICE) white 
paper: “20200602 - Culverts In Ice Bridges”. 


. The proposed OPIB is not expected to be a grounded ice bridge and is not designed to 


be a grounded ice bridge. 


. Loaded and unloaded SPMTs are capable of traversing grades up to a maximum 5% 


longitudinal and a maximum of 1.5% transverse. 
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4 DATA AND OBSERVATIONS 


The focus of the winter 2019-2020 data collection was at the proposed Upstream Site for the 
OPIB (N70.05348 W151.37277) (Figure 2.1). Three field investigations were made between 
December 31, 2019 and April 14, 2020. The timing of the field investigations was selection to 
coincide with early winter (December 31s), midwinter (February 25"), and prior the end of a 
typical ice road season (April 14). Access to the Upstream Site was overland from 
GMT2/MT7 (Figure 2.1). A Rolligon and Hagglund were used on the first 2 visits and a 
Rolligon and Tucker on the last visit. 


MBI engineers accompanied ICE engineers for the first two field visits (December 31st and 
February 25'"). COVID related travel restrictions prohibited MBI from attending the last field 
visit (April 14¢). 


Ice profiles were surveyed during each field visit (Appendix B). The following data were 
recorded during each ice profile: 


- -Span of ice - Span of free water under the ice 
- -Ice thickness - Water depth 
- -Snow depth 


The surface of the ice was the basis of elevation for each of the ice profiles instead of British 
Petroleum Mean Sea Level (BPMSL). It is common practice to reference ice profiles to the 
water surface elevation. BPMSL ice elevations are not necessary during the early phases of ice 
bridge design. 


The data collected from the ice profiles are necessary for calculating the following 
information: 

1. Crossing Cross-sectional area 

2. Construction Quantities 

3. Direct Discharge (Calculated by MBI) 


Figure 4.1 presents a sketch of the layout plan at the Upstream Site. In addition to ice profiles, 
pressure transducers (PT) were installed approximately 100 ft upstream (US) and 
downstream (DS) of the proposed Upstream Site centerline alignment. The PTs were installed 
on the channel bottom in the deepest part of the cross-section. The PTs measured absolute 
pressure which was translated into water depth by atmospheric (ATM) pressure corrections. 
The sample frequency of the PTs was set to 12-hour intervals (noon and midnight). 


The CWAT was about 85 ft upstream and paralleled the ice profile alignments. The CWAT 
was not in place during the December 31** field visit. 
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SKETCH OF UPSTREAM SITE PLAN 


UPSTREAM PT 


PROPOSED 
CENTERLINE 


ICE PROFILE ~~ 
ALIGNMENT 
DOWNSTREAM PT 


ee LOQD’ hp a lia 


SPAN OF ICE 


oe 


EDGE OF ICE 


EDGE OF ICE 


NOTES: 
- ALL DIMENSIONS ARE APPROXIMATE 
2. CWAT WAS NOT [NM PLACE DURING DECEMBER 31ST FIELD VISIT 


FIGURE 4.1: SKETCH OF UPSTREAM SITE PLAN 
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4.1 CUMULATIVE FREEZING DEGREE-DAYS 


Cumulative Freezing Degree-Days (CFDD) are calculated as a sum of average daily degrees 
below freezing for a specified time period and are frequently used to measure and compare 
the coldness of winter from year to year. The annual CFDD has ranged from a high of 9,300°F 
(2011-2012 winter season) to a low of 6,200°F (2017-2018 winter season). Generally speaking, 
the higher the CFDD then the colder a winter was. However, a higher CFDD doesn’t 


necessarily equate to a longer winter. Figure 4.2 presents the historical CFDD from 
2002 to 2020. 


NPRA NE TUNDRA MONITORING STATION 
CUMULATIVE FREEZING DEGREE-DAYS 
WINTERS 2002 - 2020 


——NE-2 2011-2012 Coest | 
somes NV 2012-2013 
—— NWI 2003-2004 
——— NWT 2005-2006 
—— NV-1 2004-2002 
NE -2 2008-2008 
—— N41 2007-2008 
—— NH 2006-2007 
——— NT 200-2005 
—— Ni 2005-2010 
wee NED 29135-2020 
NWI 201204 
—— NW 22-2003 
Ni 201 S-20t4 
ee NW 2924-2015 
—— NE2 DESDE 
aoe NE? 2015-2017 
ame NE? 2018-2018 
NE2 2017-2016 Warmest 


S017 CFDD 
25 May 2020 
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FIGURE 4.2: CUMULATIVE FREEZING DEGREE-DAYS 2002 - 2020 


The air temperature CFDD index relates to: 
1. Natural ice growth 
2. Ice construction rates 
3. Refreezing of seasonal thawed tundra 
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4.2 SPAN OF ICE 


The span of ice is the distance between the two ice edges. The span of ice at the proposed 
Upstream Site crossing was approximately 1,000 ft (Figure 4.1). It was difficult to establish the 
precise edges of the ice the interfingering of ice and sediments made it difficult to establish 
the exact ice edge location. Generally, the river ice was blown clean of snowdrifts (snow 
depths <0.1 ft). 


4.3. NATURAL ICE GROWTH 


Table 4.1 presents the ice natural ice thickness at the proposed Upstream Site OPIB during the 
three field visits. Also presented is the ice growth rate. The ice growth rate was likely the 
highest between 25 February and 15 April. Historically, natural ice growth after mid-April is 
relatively low. This is due to the increase in average daily temperatures and the increase in 
daily solar radiation. The reduction in natural ice growth makes ice bridge repairs and 
construction difficult. 


TABLE 4.1: NATURAL ICE THICKNESS AND GROWTH RATE 


Average Floating Ice | Growth Rate From Previous Growth Rate From Previous 
Thickness (feet ) Field Visit (feet per day) Field Visit (feet per week) 


2k 


“Day 1 is set to the day that CFDD > 1 (September 24, 2019) 


Natural river ice growth rates depend on air and water temperatures, water velocity, 
overflow, wind speeds and direction, ice snow cover, and ice thickness. Construction 
activities and techniques can modify ice growth to control production rates 


4.4 SPAN OF FREE WATER 


The span of free water under the ice at a given crossing is the distance measured between the 
edges of water below the ice. Figure 4.3 provides an illustration of how the span of free water 
is measured. 


Table 4.2 presents the span of free water under the ice during each of the three field visits. 
Generally, as the ice thickness increases the span of free water and the maximum free water 
column height decreases. 


During the April 14 field visit, the natural ice had come into contact with the river bottom 
about mid-way across the channel. This created two independent spans of free water 
measuring 117 ft and 94 ft wide (Appendix B - April 14, 2020 Ocean Point Crossing Profile). 


4 June 1, 2020 
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Ocean Point Crossing Profile 
December 31, 2019 
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FIGURE 4.3: ICE PROFILE AND SPAN OF FREE WATER AREA 


TABLE 4.2: SPAN OF FREE WATER 


: Nf Span Of Free Water Under | Maximum Free Water Column 
Field Visit Date : 
The Ice (feet) Height Under Ice (feet) 


50 .O 


*This is the combination of two independent spans of 117 ft and 94 ft that were 


separated by 94 ft of grounded ice 


4.5 DISCHARGE MEASUREMENTS 


MBI conducted a direct discharge measurement during the first two field visits (Appendix C). 
A summary of the MBI discharge data are presented in Table 4.3 and Figure 4.4. 


TABLE 4.3: MBI DISCHARGE MEASUREMENTS 


; ay Discharge Discharge Average Velocity 
Field Visit Date 3 
(ft"/s) (GPM) (ft/s) 


April 14, 2020 *Data not obtained due to travel restrictions 
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The average water velocities and discharge decreased as winter progressed. Unlike the 
Mackenzie and Yukon Rivers, the Colville River is classified as an Arctic River. This means 
winter flow stops since the watershed is frozen. However, some flow continues from ground 
water as overflow. 


4.6 WINTER WATER DEPTHS 


Generally, water depth in a river decreases as the tributary discharge decreases with freeze- 
up. However, in the case of Ocean Point the overall water depth increased slightly as winter 
progressed. 


Figure 4.4 illustrates the PT water levels and the atmospheric pressure. Potential overflow 
events were likely recorded by the pressure increases during February and early March 
(Section 4.8). The water levels increased by late March at a rate of about 0.2 ft in 30 days. 
Overall water depths were less than 5.0 ft with the maximum 5.6 ft in February. 
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Water Depths Near The Proposed OPIB 


(ATM Adjusted) 
Jan - April 2020 
Qi). 5S Sa et al 3S 
5.8 H February 25, 2020 April 14, 2020 | 
Discharge = 135 ft?/s (60,540 GPM) ees 
5.7 December 31, 2019 
SAS ta om at 
= 
- 
o 
6 
= 
-_ 
° 
‘ s 
2 
a 
[=] 
38 | | Notes: : 
u T - There is approximately 200 feet between pressure transducers 32.8 
37. Se - Upstream PT N70.05386 W151.37283; Downstream PT N70.05328 W151.37259 20 u7: 
HET SE SAAT peice - Atmospheric pressure data obtained from the DS-2P weather station 
3.6 BELTS LCT ANCOR, [os a Shen T - Pressure transducers setto record at 12 hour intervals 32.6 
Raerees ATM Pressure 
3.5.4 7 7 T Sass ste = t ; 32.5 
S S > S > S S S S S S S S iS) 
\y “vy \v Wy “vy “vy (Vv “v “v wy w “vy iy “vy iv “vy 
Cr ee eS ed Le a TU ee ee 
YY Y YY iA wv ov wv No S 2 9 Ne > > 
) Date 
) FIGURE 4.4: PT WATER DEPTHS NEAR THE PROPOSED OPIB 
) 
) TABLE 4.4: MAXIMUM WATER DEPTHS 


f c Maximum Water Depth 
) Field Visit Date 
| (feet) 


8 
) February 25, 2020 6.1 
) April 14, 2020 67 
| 


Possible explanations of the water depths whether rising or falling are a complex system of 
t many unknown variables. A study of all the possible variations is beyond this report. 
} Increasing ice thickness can constrict the flow and cause an increase in recorded pressure. 
During the coldest weeks of the winter, ground waters are emitted in larger volumes, but the 
) waters freeze before moving far from the effluent. This makes the casual observations see 


| what appears to be low flow. 


4.7 WATER TEMPERATURES 


In addition to pressure, water temperatures were recorded. Figure 4.5 presents the river 


bottom water temperatures and the average daily air temperature as recorded by the Drill 
} Site 2P weather station (plotted on the right vertical axis). What is noteworthy in the data are 
the following: 
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1. There are spike-increases in the river bottom temperatures on February 28 and 
March 24th. These are likely related to overflow events and are discussed in further 
detail in Section 4.8. 

2. For about 5 days starting on March 19% there was an event that caused the temperature 
to rise and fall -/+ 0.2°F at both PT locations. This is also likely related to an overflow 
event (Section 4.8) in addition to a potential phenomenon that is explained below in 
Section 4.7.3. 

3. Daily temperature oscillations are visible in the data beginning on March 19. These 
undulations are likely due to solar diurnal rhythms. The first annual sunrise at the 
OPIB latitude began on January 17 and by March17 the daily sun exposure 
exceeded 12 hours per day. 


River Bottom Water Temperature - Upstream and Downstrem PT Ocean Point 
Jan - April 2020 


Field Vist 


February 25, 2020 Field Vist 
i i = ia April 14, 2020 


__ Air Temperature [°F] 


Water Temperature [°F] 


——»-— Upstream Water Temperature 


Average Water Temperature 


—+— Downstream Water Temperature 


----e--=- Average Daily Air Temperature 


FIGURE 4.5: RIVER WATER TEMPERATURE — OCEAN POINT JANUARY TO APRIL 2020 


4.8 OVERFLOW CONSIDERATIONS, CHARACTERISTICS, AND OBSERVATIONS 


A unique design consideration for the proposed Upstream Site OPIB is the potential for 
overflow. There is a high likelihood for overflow on the Colville River during any given 
winter. The frequency, location, triggers, and magnitude of overflow events can be difficult to 
predict and measure. Considerable attention was made during the data collection to record 
overflow events with instrumentation and observations. There were at least three potential 
overflow events during the winter of 2019-2020. 
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4.8.1 GENERAL OVERFLOW CONSIDERATIONS AT OCEAN POINT 

Early in winter (November - December) overflow events are not typical. The thin ice (<2 ft) 
rises and falls with the changes in water levels. Discharge declines as the rains give water to 
snowfall and the watershed freezes. When overflow events do occur, they tend to be of lesser 
magnitude. Normally, when an overflow event occurs it breaches along the edges of the river 
where the ice has not become firmly grounded. 


By the middle of winter (January - February) the potential for overflow events increase in 
frequency and magnitude. The river ice has become firmly grounded along the edges of the 
river and the cross-sectional area of free water is substantially reduced (80% reduction from 
the OPIB December 31s to the February 25th ice profile). This can lead to higher pressures, 
constrictions, and increased water velocity under the ice. Eventually the pressure becomes 
great enough to form cracks in the ice. The flowing water finds pathways through the cracks 
to the surface of the ice. The overflow from these events tends to flow in all directions on the 
surface of the river ice. These events can be difficult to observe if there is snow on the surface 
of the ice. Furthermore, the presence of snow can increase the amount of time required for the 
overflow to freeze. 


Toward the end of winter (March - April), the potential for overflow is similar to that of the 
middle of winter with the additional contribution from increased solar radiation intensity and 
warmer air temperatures. Overflow from snowmelt due to solar radiation and warmer 
temperatures tend to result in ponding and minimal flow in any particular direction. 
Generally, there is little potential for runoff during this time period. However, runoff should 
be expected during the month of May. 


4.8.2 OVERFLOW OBSERVATIONS SUMMARY 

Table 4.5 summarizes the three potential overflow events recorded by the PTs. Figure 4.6 
presents the relative water depth differential between the US and DS PTs. This difference is 
the DS PT water depth subtracted from the US water depth after each PT water level has been 
set to zero for its initial reading. Positive values on the graph indicate that the water levels 
were greater at the US PT than the DS PT. A water depth differential along the horizontal 
dashed line means that the water levels at the US and DS PTs are back to their initial 
(December 31, 2020) readings. Negative values on the graph indicate that the DS PT recorded 
water levels greater than the US PT. 
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TABLE 4.5: SUMMARY OF POTENTIAL OVERFLOW EVENTS 


Evidence of this potential event was not observed 
January 18, 2020 ; rg 

during the February 25 field visit 
February 28, 2020 The largest potential event recorded 


CWAT reported event and travel suspended due to 
March 24, 2020 0 
flooding 


NOTE: The CWAT was located approximately 85 ft upstream of the ice profile centerline 


Water Depth Differential Between Upstream/Downstream PT 
(ATM Adjusted) 
Jan - April 2020 
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FIGURE 4.6: WATER DEPTH DIFFERENTIAL BETWEEN UPSTREAM AND DOWNSTREAM PTS 


Once the potential February overflow event ended the water temperatures at the US and DS 
PTs returned to their typical differentials. It is likely that flow resumed under the ice. The 
overflow likely flowed toward the banks of the river and froze. 


Evidence of this overflow event was observed during the April 14 field visit. While 
conducting the ice profile there was up to 0.8 ft of frozen overflow observed at the eastern 
and western edges of the profile. Frozen overflow is differentiated from natural ice in that it is 
opaque and rests atop natural ice that is clear. Frozen overflow was not observed between 
STA 3+00 and 6+00 of the ice profile. 
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An overflow event was observed and reported by CWAT users around March 24th. The 
CWAT snow bridge was located about 100 ft upstream from the ice profile centerline. The 
report stated that the snow bridge was flooded and deemed unsafe. This event was likely 
intensified by the warming air temperatures and the increased intensity of solar radiation. 
The overflow from this event was frozen by the time of the April 14 field visit. 


4.9 OCEAN POINT WEATHER STATION INSTALLATION 


A weather station was installed on April 14, 2020 2.4 miles north of the proposed Upstream 
Site for the OPIB (N70.08730 W151.35590). Table 4.6 presents a summary of the parameters 
recorded by the Ocean Point weather station. 


TABLE 4.6: OCEAN POINT WEATHER STATION INSTRUMENTATION 


*Planned installation for fall 2020 


The weather station summary report is included in Appendix D. 
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5 OCEAN POINT ICE BRIDGE DESIGN 


The OPIB design and construction quantities have been updated based on the data and 
observations from the winter of 2019-2020 investigation (Section 4). Further revisions will be 
issued as more information is obtained. The revised ice bridge design (Appendix E) includes 
new design features based on recent discoveries (overflow, ice profiles, etc.). The subsections 
below provide explanations for the updates and revisions of the OPIB design. 


5.1 OVERFLOW MONITORING 


Any overflow, if it occurs, should be monitored beginning in December during the first year 
of heavy haul use (winter 2024-2025). Monitoring stations with remote capabilities would be 
installed prior to construction of the ice bridge and removed before breakup. Overflow 
monitoring will assist with mitigation management outlined in Section 5.2. 


5.2, OVERFLOW MITIGATION MANAGEMENT 


Overflow mitigation should be included in the OPIB design. Overflow events may impact 
the construction and use of the OPIB (Section 4.8). 


Overflow mitigation actions can be classified into two broad categories; passive and active. 
Both mitigation requirements should be incorporated into any OPIB design. 


5.2.1 PASSIVE OVERFLOW MITIGATION 

Passive overflow mitigation are preventative actions put into place to prevent overflow at an 
ice structure. Passive mitigation should not require constant attention and maintenance to 
remain effective; although routine maintenance may be required. Passive mitigation includes 
berms, passageways, and other diversion structures. Passive mitigation designs may require 
some cleanup and recovery depending on the impacts and magnitude of an overflow event. 


5.2.2 ACTIVE OVERFLOW MITIGATION 

Active overflow mitigation operates in combination with passive mitigation. Active overflow 
mitigation are preventative actions established if the passive mitigation systems become 
insufficient. Due to the likelihood of overflow events, active overflow mitigation will be 
necessary for the heavy haul ice bridge over the Colville at Ocean Point. 


Two available forms of active mitigation are high-volumetric flow pumps and designated 
rapid response heavy equipment. The pumps would be staged at the site to move water from 
the upstream side of the ice structures to the downstream side if passive mitigation becomes 
insufficient. Rapid response heavy equipment would be used to clear new pathways for 
water to flow away from the ice structures. Furthermore, the heavy equipment would be 
utilized to repair, maintain, and restore passive mitigation structures. 


te June 1, 2020 
pls Page 5-1 


Willow Development — Ocean Point Ice Bridge Revision 
tenement 


All active mitigation measures must be readily available to assist with overflow management. 
A comprehensive ‘overflow management plan’ should be developed as more data is collected 
in the vicinity of the OPIB. 


5.3. RAMP DESIGNS 


The revised ramp designs of the proposed OPIB have been modified to conform to the 
maximum design specification of 5%. The 3% grade of the central ramps will ensure that the 
ramps have overall lengths greater than that of the maximum specified SPMT length (220 ft) 
to avoid high-centering. 


The ramps are 65 ft wide based on the maximum ice fill thickness. The east ramp width may 
increase depending upon the topographic survey of the east bank. 


5.4 ICE BRIDGE CONSTRUCTION QUANTITIES 


Table 5.1 presents a summary of the ice material quantities required to construct the OPIB 
based on the most recent information. The average natural floating ice thickness is expected 
to be between 4 and 5 ft by February. 


TABLE 5.1: SUMMARY OF MATERIAL QUANTITIES 


*The water equivalent quantity is 6.6 million gallons 


Most of the construction material is required to build the ice ramps and the elevated ice 
bridge over the span of free water. In addition to supporting the loads, the elevated ice bridge 
performs other key functions, which include: 


1. Elevate the driving surface away from the overflow zone. 

2. Increase the maximum allowable width of the free water under the bridge by 
increasing the pressure area at the base of the ice sheet. 

3. Reduce the amount of snow that accumulates on the driving surfaces. 


5.5 EMERGENCY BYPASS / ACCESS ROAD AND RAMPS 


An emergency bypass/access roads and ramps should be constructed on the downstream 
side of the elevated ice bridge and ramps. These roads will provide access around the SPMTs 
while they are navigating ice structures when crossing the OPIB. Specific design details of the 
emergency bypass/access roads will be provided later. 
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rr 
6 CONCLUSIONS & RECOMMENDATIONS 


Conclusions and recommendations are subject to change as new information becomes 
available. 


6.1 CONCLUSIONS 


1. The proposed Upstream Site OPIB will be a non-grounded ice bridge with a 
capacity that is suitable for the 3,200 ton module loaded onto a SPMT with a 4,200- 
ton allowable gross weight. 

2. There is a high probability that at least one overflow event will occur in the vicinity 
of Ocean Point each winter. 

3. More information regarding water discharge and water levels in the vicinity of the 
Ocean Point between the months of February and mid-April is needed. This is the 
timeframe of construction and heavy haul use. 


6.2 RECOMMENDATIONS 


1. MBI should conduct weekly discharge measurements at the OPIB from mid- 
February to mid-April during the winter of 2020-2021. 

2. Ice profiles at the OPIB should be performed on a monthly basis during the winter 
of 2020-2021, with weekly ice profiles between mid-February to mid-April. 

3. Delay any geotechnical investigation until 2023-2024; this task may not be 
necessary. 

4. Install remote monitoring sites at the OPIB to collect water level, air temperature, 
and air pressure throughout the winter of 2020-2021. 

5. Perform a topographic survey from the top of bank to the top of bank along the 
heavy haul ice road alignment at the Upstream Site during the winter of 2020-2021. 
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Ocean Point Crossing Profile 
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Ocean Point Crossing Profile (Looking Downstream) 


April 14, 2020 
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1.0 INTRODUCTION 


Michael Baker International (Michael Baker) collected water resources data for Conoco Phillips Alaska, Inc. 
(COPA) in support of the Willow Project (Willow). Two proposed crossings of the Colville River were 
investigated at Ocean Point. Data was collected during three field events occurring between Fall 2019 and 
Spring 2020. This report summarizes the methods and results of that effort. 


ICE Design & Consult (ICE), UMIAQ, LLC (UMIAQ), Soloy Helicopters, Peak Oilfield Services Company 
(Peak) and CPAI Alpine Helicopter and Field Environmental Coordinators provided support during the field 
program and contributed to a safe and productive field season. 


2.0 LOCATION 


Two transects near Ocean Point were investigated: Ocean Point Downstream (Downstream) and Ocean 
Point Upstream (Upstream, Figure 1). These are two among six transects investigated during the 2018- 
2019 ice road season. The Downstream and Upstream transects were selected based on shallow water 
depths relative to the other transects investigated. Ocean Point Downstream is approximately 6.5 miles 
direct or 8.3 river miles (RM) downstream of Ocean Point Upstream. There is one minor tributary that enters 
the Colville River from the south between the two locations. This tributary is a meandering beaded stream 
that drains multiple lakes. There is also a tributary that enters the Colville River from the northwest at the 
Downstream transect. This tributary is a paleochannel which drains a series of lakes that formed in 
abandoned meanders after the reach of river between the Downstream and Upstream transects migrated 
south. This area is inundated during spring breakup-induced flooding. 


Ocean Point Downstream (Figure 2, also historically referred to as “Transect #6”, the “east crossing”, the 
“downstream crossing”, or “Ocean Point North”) is the alternate proposed crossing location. Ocean Point 
Upstream (Figure 3, also historically referred to as “Transect #1”, the “Rolligon crossing”, the “west 


crossing”, the “upstream crossing’, “Ocean Point South”) is an historic crossing location. It was the location 
of a snow road during the 2018-2019 season and is the preferred proposed crossing location. 


Table 1 provides a summary of dates and data collected at the locations investigated. Table 2 provides a 
Table 1: Field Events 


summary of measurements collected. 


: g 
Water Quality 


Water Surface Elevation Survey 


Bank Active Layer Investigation 


Notes: 
1. Crossing characterized as similar to Ocean Point Upstream so winter investigations were not performed. 
2. Planned field event was cancelled. 


Data Collection 


9/5/2019 
12/31/2019 
2/25/2020 
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Table 2: Data Collected 


water depth, under ice? feet 


flow width 


flow cross-sectional area square feet 


velocity feet per second 


ft/s 


dissolved oxygen percent saturation 
dissolved oxygen milligrams per liter 


water surface elevation feet North American Vertical ft NAVD8&8& 


Datum of 1988 
Bank Active Layer thawed soil depth feet 
Notes: 


1. Data collected only in winter 


Water Surface 
Elevation 
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3.0 METHODS 


Field sampling methods were based on United States Geological Survey (USGS 2006a and 2006b) 
methods. Safety precautions were followed using the North Slope Water Resources 2019 Health, Safety, 
and Environment Plan (Michael Baker 2019a) and the 2019-2020 Winter Hydrology Programs — Job Safety 
Analysis (Michael Baker 2019b). 


3.1 OPEN WATER - FALL 2019 


Open water tasks were completed over the course of multiple days during one field event and were timed 
to occur late in the fall season, prior to freeze up. Sites were accessed via helicopter and inflatable rafts 
with outboard motors. The effort was conducted by a three-person field crew. The data collected includes 
cross-sectional river bottom profiles, discharge, velocity, water depth, water surface elevation, site 
conditions, and in-situ water quality. Soil active layer depths were also investigated for both banks of each 
crossing. 


Discharge, velocity, and cross-sectional river bottom profiles were measured using a RiverRay acoustic 
doppler current profiler (ADCP) (Photo 3.1). The ADCP was mounted in a trimaran. The trimaran was 
tethered via boom to the side of a 13-ft Achilles inflatable raft powered by outboard motor (Photo 3.2). 


Water surface elevation (WSE) at the time of 
discharge and water quality measurements were 
determined using temporary benchmarks installed 
by UMIAQ surveyors and Michael Baker. UMIAQ 
benchmarks are aluminum cap survey control 
referenced to OPUS NAVD88 elevation. Michael 
Baker benchmarks are aluminum cap survey 
control tied by level loop technique to the UMIAQ 
control. 


In-situ water quality parameters including 
temperature, conductivity (C), and salinity were 
recorded using the YSI ProPlus meter; dissolved Photo 3.1: Discharge measurement at Ocean Point 
oxygen (DO) was measured using the YS! ProODO Upstream Transect; 9/5/19 

meter. Specific conductance (SC) referenced to 25 . 
degrees Celsius was _ calculated based on 
temperature, conductivity, and a _ conversion 
coefficient of 0.0196 based on empirical data. 
Measurements were collected at the deepest 
portion of each cross-section; two depths were 
investigated at each to confirm that parameters 
were consistent throughout the water column. 


Photo 3.2: In-situ water quality sampling at Ocean 
Point Downstream Transect; 9/5/19 
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The soil active layer was_ investigated to 
characterize the depth of permafrost below the 
surface in the transition zones adjacent to the river. 
This investigation was performed on the right and 
left banks of each transect, perpendicular to the 
channel. Probing was performed using a 5-foot long 
T-bar probe driven by hand (Photo 3.3). Spacing 
was approximately 5-foot increments and at major 
grade breaks, between water’s edge and the top of 
the riverbank. Results were provided directly to 
Golder as a separate deliverable and included in 
this report along with discharge and water quality in 
Attachment A. 


Photo 3.3: Investigating soil active layer at Ocean 
Point Downstream Transect; right/east bank; 
9/14/19 


3.2 ICE COVER - WINTER 2019 - 2020 


Three field events were planned to investigate the trend in flow quantity and water quality over the course 
of the ice-cover season. Freeze-up typically initiates in mid-October and breakup typically initiates in mid- 
May. Ice cover field events were one day apiece. The first was performed early in the season, the second 
in the middle of the season, and a third was planned at the end of the season. The third field event was 
cancelled due to circumstances related to COVID-19 and changing project priorities. Data was collected at 
one transect, Ocean Point Upstream. This included under-ice cross-sectional river bottom profiles, 
discharge, velocity, water depth, ice thickness, water surface elevation, site conditions related to overflow, 
and in-situ water quality. ; 


A one-person Michael Baker field crew conducted both events, supported by an ICE engineer who 
performed crossing profiling. UMIAQ and Peak provided transportation to the sampling locations and 
general field support. The sites were accessed by Hagglund and Rolligon. 


Water measurements were facilitated by 
mechanically drilling through the river surface 
ice cover. Thermal drill probing was performed 
by ICE to identify the extents of under-ice water 
bounded at the left and right by ice grounded 
against the channel bed. Investigation of soils 
or groundwater within the channel bed was not 
performed. Discharge was determined using 
USGS mid-section techniques. Velocity was 
measured using a handheld Hach flow meter 
(Photo 3.4) and a handheld Sontek flow meter. 
These were attached to a fixed rod and lowered 
to 0.6 the water depth below the ice. In-situ 
water quality parameters investigated were the 
same as those in the fall. Measurements 
were collected at multiple depths from one 
location in the deepest portion of the cross- _ photo 3.4: Attaching Hach flow meter to fixed rod at Ocean 
section. Point Upstream; 2/25/20 


Previously submitted ice cover season field data is provided in Attachment C and Attachment D. 
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4.0 RESULTS AND CONCLUSIONS 


A summary of Colville River Ocean Point water resources information is provided below. Previously 
submitted trip reports and field data are provided in Attachment B through Attachment D. 


4.1 TRANSECT LOCATIONS 


The Ocean Point Downstream transect is located where a tributary enters the Colville River on the left/west 
bank. It is unknown what, if any, flow is contributed from this tributary during the ice-cover season. Bankfull 
width at this location is approximately 2,500 feet based on aerial imagery. This reach of the Colville River 
is relatively straight. Bars are exposed during low water. The low water channel lies approximately central 
to the cross section and the thalweg lies toward the left/west portion of the low water channel. Both banks 
are steep. The right/east bank is steeper than the left/west with sloughing and block failure; evidence of 
both thermal and mechanical erosion. The tops of both banks and overbanks are vegetated. Vegetation is 
present on the upper left/west bank below the top. Vegetation is present on the right/east bank. 


The Ocean Point Upstream transect is located where the Colville River is conveyed within a single channel 
with a bankfull width of approximately 3,200 feet based on aerial imagery. The Colville River transitions 
from relatively straight to a wide bend at this location. The left/east bank is on the cut-bank inside of the 
bend and the right/west bank is on the point-bar outside of the bend. Bars are exposed during low water. 
The low water channel is located closer to the left/east bank and the thalweg lies toward the in the right/west 
portion of the low water channel. Both banks are steep. The left/east bank is steeper than the right/west 
bank. The tops of both banks and overbanks are vegetated. Vegetation is present on the upper left/east 
bank below the top. Vegetation is present on the right/west bank. 


No springs were observed at the bank of either transect during the open water field event. No overflow, 
aufeis, or evidence of any other notable hydraulic occurrence was observed at Ocean Point Upstream 
during the ice-cover field events. Open water data collected at Ocean Point Downstream were compared 
against those collected at Ocean Point Upstream. With respect to discharge quantity and water quality 
parameters, values were similar between the two locations. 


Ice-cover investigation was performed only at Ocean 
Point Upstream (Photo 4.1). This was determined 
considering: the comparison of open water discharge 
quantity and water quality values, the better suitability of 
the Upstream transect to the Downstream transect based 
on the potential for undesirable geomorphological and 
hydraulic influences (i.e. actively eroding and sloughing 
bank and a tributary at Ocean Point Downstream), the 
remoteness of Ocean Point relative to facilities and the 
challenges of winter accessibility, the time available for 
Bi ele investigation, and the historical use of the upstream 
ee aoe | transect. 


Photo 4.1: Field investigation at Ocean Point 
Upstream, looking toward the left/east bank; 
2/25/20 
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4.2 BANK ACTIVE LAYER 


The active layer was investigated in September, at the end of the thawing season. Elevations were not 
surveyed. Depths of thawed soil, as measured by probing to refusal, ranged from 1.75 feet to greater than 
5 feet. Thaw depths were shallower at the tops of banks and deeper approaching the channel. A summary 
of thaw depths is provided in Table 3. Approximate locations probed are provided in Figure 2 and Figure 3. 


Table 3: Colville River Ocean Point Bank Thaw Depth Summary 


Thaw depths of 5 feet indicated frozen ground was not encountered. These were assumed to be 5 feet for 
calculated averages though actual values are greater. 


4.3 PHYSICAL WATER MEASUREMENTS 


Colville River open water discharge measured at the Ocean Point Downstream transect was within 200 cfs, 
or 0.7% difference, of the discharge measured at the Ocean Point Upstream transect. The errors associated 
with each discharge measurement transect (two at the Downstream transect and four at the Upstream 
transect) were less than +/-1.5%. The total average error associated with discharge measured at each 
location was 0.0%, with a standard deviation of 0.3% at Ocean Point Downstream and 1.1% at Ocean Point 
Upstream. 


Measuring discharge under ice cover is subject to limitations not applicable to open water measurements. 
Unlike open water where it is obvious where the edge of water exists, it is not possible to see the extents 
of the cross-sectional area of flow under the ice. Further, it is not possible to profile the entire cross- section. 
It is assumed that the cross-sectional area is relatively uniform upstream of, downstream of, and between 
measurement stations. However, the potential exists for “unseen” grounded or relatively shallow areas 
which would influence measured velocity direction and magnitude if occurring upstream or downstream of 
a measurement station. Grounded areas between measurement stations would reduce the estimated cross- 
sectional area of flow and resulting discharge. Colville River discharge measured at Ocean Point Upstream 
during the ice-cover season was significantly less than discharge measured during the open water season. 
Discharge decreased as the ice-cover season progressed. 


This decreasing trend is also apparent in the Colville River at Umiat where a continuous gage station is 
operated by the USGS (USGS 2020). This location is approximately 70 RM upstream of Ocean Point 
(Figure 4). The drainage area between is expansive, including multiple large tributaries as well as an 
unknown quantity of groundwater springs. Despite this, general comparison regarding seasonal discharge 
trends can be made. 
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The USGS Colville River gage at Umiat (USGS 15875000 COLVILLE R AT UMIAT AK) was established in 
1953. Hydrologic stage and discharge data are available from this site. Values provided are historical and 
current; they are alternatively measured, instantaneous, time-averaged, and _ statistical. Direct 
measurements validate calculated results. River hydraulics and environmental factors differ between the 
open water and ice-cover seasons, which necessitates different approaches to data collection and 
calculation during each. 


During the open water season, instantaneous stage at Umiat is measured and provided. Instantaneous 
discharge is determined based on stage using a stage-discharge rating curve. The rating curve was 
developed by plotting measured stage events against measured concurrent discharge events. The 
accuracy of the rating curve is directly proportional to the accuracy and quantity of the measurements used 
to plot it. Numerous factors affect discharge measurements including temporal site conditions, equipment 
and technique used, and experience of the hydrologist (USGS 2010). USGS evaluates discharge 
measurements qualitatively by the ratings “excellent” (within 2%), “good” (within 5%), “fair” (within 8%), and 
“poor” (greater than 8%). Since 2002, 153 direct discharge measurements have been made by the USGS 
at the Umiat gage site. Of these, 102 have occurred during the open water season. Of those, 24% were 
rated “good”, 61% were rated “fair’, and 18% were rated “poor”. None were rated “excellent” and one was 
not rated. Open water time-averaged and statistical values, i.e. daily means, mean of daily means, peaks, 
etc. for stage and discharge are determined based on instantaneous and measured values. 


Instantaneous, time-averaged, and statistical stage values at Umiat are not provided during the ice-cover 
season. Time-averaged and statistical discharge values are provided, however. Measured stage and 
discharge values are also provided. There have been 37 direct measurements performed under the 
influence of surface ice cover. Of those, 5% were rated “good”, 16% were rated “fair”, and 76% were rated 
“poor’. None were rated “excellent” and one was not rated. Daily mean discharge is determined not based 
on stage, but instead on storage depletion modeling based on time and using a low-flow value immediately 
prior to freeze-up as the controlling factor. Umiat daily mean discharges for the 2019-2020 ice-cover season 
have yet to be validated and made available. The mean of daily mean values, however, are available for 
comparison to measurements collected at Ocean Point. The period of record informing those is between 
October 1, 2001 and September 30, 2019. 


Colville River discharges measured at Ocean Point are provided in Table 4. Colville River discharges at 
Umiat are provided for comparison. 


- = ° 
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Notes: 

1. USGS Gage 15875000 COLVILLE R AT UMIAT AK 

2. Between Colville River Umiat and Ocean Point Upstream 

3. Daily mean discharge record; mean of daily mean discharge is 20,100 cfs 

4. No daily mean discharge record yet available for this date; value is mean of daily mean between 
10/01/2001 and 9/30/2020 


Table 4: Colville River Discharge 


(cfs) 
Umiat? 
mean discharge (cfs) 


easured discharge 
(cfs) 


Ocean Point 
Downstream 
Ocean Point 
Downstream 


Ocean Point Upstream 
measured discharge 
cean Point Upstream 
Percent difference? 

% 
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The lowest annual mean of daily mean discharges for the Colville River at Umiat is 2.9 cfs. This occurs 
between April 13 and April 21. Discharge is similarly low during the month of April and relatively low 


throughout the winter 


season. Discharge increases by orders of magnitude as breakup processes initiate 


in mid-May, peaking at the end of May. These data are provided in Table 5 and graphically in Chart 1 and 


Chart 2. 


Table 5: Mean of Daily Mean Discharge Values for Colville River at Umiat 


00060, Discharge, cubic feet per second, 


Mean of daily mean values for each day for water year of record in, ft3/s (Calculation Period 2001-10-01 -> 2019-09-30 
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Chart 1: USGS Mean of Daily Mean Discharge Values for Colville River at Umiat - Annual 
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Chart 2: USGS Mean of Daily Mean Discharge Values for Colville River at Umiat - December through May 
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It is reasonable to assume that Colville River seasonal discharge trends at Ocean Point are similar to those 
at Umiat. Direct measurements at Ocean Point support this conclusion. Lacking evidence suggesting 
otherwise, it is further reasonable to assume that mid-April discharge at Ocean Point, had it been measured, 
would have been lower than the 9 cfs measured in February. 


A summary of additional physical parameters measured during the project performance period are provided 
in Table 6 through Table 9. 


Table 6: Colville River Ocean Point Discharge and Water Surface Elevation Summary 


Ocean Point Downstream Ocean Point Upstream 
Ocean Point Downstream water surface elevation Ocean Point Upstream water surface elevation 
measured discharge (cfs) (ft NAVD88) measured discharge (cfs) (ft NAVD88) 


(0 890 300 
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Table 7: Colville River Ocean Point Velocity Summary 


Ocean Point Downstream Ocean Point Downstream Ocean Point Upstream Ocean Point Upstream 
measured velocity (ft/s) measured velocity (ft/s) measured velocity (ft/s) measured velocity (ft/s) 
maximum average maximum average 
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Table 8: Colville River Ocean Point Water Depth Summary 


Ocean Point Downstream 
effective water depth (ft/s) 
maximum 
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rune | 
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Ocean Point Downstream Ocean Point Upstream Ocean Point Upstream 
effective water depth (ft/s) | effective water depth (ft/s) | effective water depth (ft/s) 
average maximum average 
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Table 9: Colville River Ocean Point Flow Width Summary 

Ocean Point Downstream Ocean Point Upstream 
flow width (ft) flow width (ft) 
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4.4 Water Quality Measurements 


Salinity and conductivity measurements throughout the monitoring period suggest this location is upstream 
of coastal influence. While values increased between the open water and ice-cover seasons, results are 
indicative of freshwater rather than a saline environment. Increases are likely attributable to concentration 
as a result of the freshwater freezing process, which readily excludes entrained materials. Temperature 
decreased between the open water and ice-cover season. Dissolved oxygen also decreased. This is typical 
of water bodies under the influence of ice cover, which prevents the introduction and mixing of atmospheric 
oxygen. 


Table 10: Colville River Ocean Point Water Quality Summary 


total depth (ft) 
temperature (°C) 
conductivity (uS/cm) 
specific conductance 
dissolved oxygen 
dissolved oxygen (%) 
salinity (ppt) 


DOWNSTREAM 


9/5/2019 10.0 


12/31/2019 
2/25/2020 


UPSTREAM 


9/5/2019 
12/31/2019 
2/25/2020 
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User-Manual- RevC.pdf 
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AttachmentA Photos - Riverbanks 


Note: Photo locations are provided in Figure 2 and Figure 3. Photo elevations are unknown. 
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| Photo A.1: Ocean Point Downstream transect left/west bank, looking upstream/southwest toward 
| tributary; 9/5/19 
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Photo A.2: Ocean Point Downstream transect right/east bank; looking upstream/southwest; 9/4/19 
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Photo A.3: Ocean Point Downstream transect right/east bank; looking toward channel/west; 9/4/19 


Photo A.4: Ocean Point Downstream transect right/east bank; looking downstream /northeast; 9 /4/19 
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Photo A.5: Ocean Point Upstream transect left/east bank 
9/5/19 
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exposed bar; looking upstream/northwest; 


Photo A.6: Ocean Point Upstream transect left/east bank; looking downs 
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tream/southeast; 9/5/19 
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Photo A.7: Ocean Point Upstream transect right/west bank; looking upstream/northwest; 9/4/19 
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Photo A.9: Ocean Point Upstream transect right/west bank exposed bar; looking upstream/northwest; 
9/5/19 
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Acronyms & Abbreviations 


Ae: 
ADCP 
cfs 

CPAI 

fps 

ft 
uS/cm 
mg/L 
Michael Baker 
NAVD88 
OPUS 


9/24/19 


degrees Celsius 

acoustic doppler current profiler 
cubic feet per second 
ConocoPhillips Alaska, Inc. 

feet per second 

feet 

microsiemens per centimeter 
milligrams per liter 

Michael Baker International 
North American vertical datum of 1988 
Online Positioning User Service 
parts per thousand 

discharge 

velocity 


water surface elevation 
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1. INTRODUCTION 


Michael Baker International (Michael Baker) collected open water data for Conoco Phillips Alaska, Inc. (CPAI) for 
the Willow ice road project. This field event was the first of multiple trips between Fall 2019 and Spring 2020; it 
occurred between September 4 and September 6, 2019. These efforts are designed to support the characterization 
of two proposed ice road crossings of the Colville River near Ocean Point. This document presents a summary of 
this field effort and the preliminary results of the data collection. 


2. LOCATIONS 


Two transects near Ocean Point were investigated: Transect #1 and Transect #6 (Figure 1). These are two among 
six transects investigated during the 2018-2019 ice road season. Both were selected based on shallow water 
depths relative to the other areas. Transect #1 (also referred to as the “rolligon crossing”, the “west crossing”, 
or the “upstream crossing”) is an historic crossing location. It was the location of a snow road during the 2018- 
2019 season and is the preferred proposed heavy haul crossing location. Transect #6 (also referred to as the 
“east crossing” or the “downstream crossing”) is the alternate proposed crossing location. 
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3. METHODS 


Michael Baker collected field data at two proposed crossing locations across the Colville River near Ocean Point. 
Data included cross-sectional river bottom profiles, discharge, velocity, water depth, water surface elevation, site 
conditions, and general in-situ water quality parameters. Soil active layer depths were also investigated for both 
banks of each crossing. 


Discharge, velocity, and cross-sectional river bottom profiles 
were measured using a RiverRay acoustic doppler current 
profiler (ADCP) (Photo 1). The ADCP was mounted in a 
trimaran. The trimaran was tethered via boom to the side of a 
Zodiac inflatable raft powered by outboard motor (Photo 2). 


In-situ water quality parameters including temperature, 
conductivity (C), and salinity were recorded using the YS| 
ProPlus meter; dissolved oxygen (DO) was measured using the 
YSI ProODO meter. Specific conductance (SC) referenced to 25 5 , 
degrees Celsius was calculated based on temperature, Photo 1. Discharge measurement in progress 
ae t oe at Transect #1; Sept 5, 2019 

conductivity, and a conversion coefficient of 0.0196 based on 

empirical data. Measurements were collected at the deepest 
portion of each cross-section; two depths were investigated at 
each to confirm that parameters were consistent throughout 
the water column. 


Water surface elevation (WSE) at the time of discharge and 
water quality measurements were determined using temporary 
benchmarks installed by UMIAQ surveyors and Michael Baker. 
UMIAQ benchmarks are aluminum cap survey control 
referenced to OPUS NAVD88 elevation. Michael Baker 
benchmarks are aluminum cap survey control tied by level loop 
technique to the UMIAQ control. 


Photo 2. In-situ water quality sampling at 
Transect #6; Sept 5, 2019 


The soil active layer was investigated to characterize the depth 

of permafrost below the surface in the transition zones adjacent to the river. This investigation was performed on 
the right and left banks of each transect, perpendicular to the channel. Probing was performed using a 5-foot 
long T-bar probe driven by hand (Photo 3). Spacing was 
approximately 5-foot increments and at major grade breaks, 
between water's edge and the top of the riverbank. Results were 
provided directly to Golder as a separate deliverable and included 
in this report as Appendix A. 


These tasks were completed over the course of multiple days 
during one field event and were timed to occur late in the fall 
season, prior to freeze up. 


Photo 3. Investigating soil active layer at 
Transect #6; left bank. Sept 4, 2019 
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4. RESULTS 


Results of the field effort are presented below. 


4.1. Colville River Ocean Point Measured Discharge 


TRANSECT #1 

Location: Transect #1/Rolligon Route 

Date & Time: September 5, 2019 2:50 PM 

Equipment: RiverRay ADCP attached to the side of an inflatable raft with outboard motor 
WSE (ft NAVD88): 11.67 

Discharge ([Q] cfs): 29,000 

Velocity ([v] fps): (oi 

Measurement Rating: Good 

Measurement Notes At the time of the measurement, open-channel conditions were present. Wind 


was negligible and surface waves were not present. Prior to deployment, 
diagnostic tests were performed, and the internal compass calibrated. 


Table 1. Colville River Ocean Point Transect #1 Measured Discharge Summary 


Measurement Starting TotalQ DeltaQ Measured atphe Measured oes soe ee cee 


Transect # Bank (cfs) (%) Q (cfs) aie yi 4 Width (ft) ‘Cane joe eee epripael, 
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A. Colville River Ocean Point Transect#1/Rolligon Route: Measurement Transect 001 Raw 
Data Output 
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B. Colville River Ocean Point Transect#1/Rolligon Route: Measurement Transect 002 Raw 
Data Output 
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C. Colville River Ocean Point Transect#1/Rolligon Route: Measurement Transect 003 Raw 
Data Output 


5,000 {1/5} 
1th) 


Distance North (Ref BT 


D. Colville River Ocean Point Transect#1/Rolligon Route: Measurement Transect 004 Raw 
’ Data Output 


5,000 |R/s} 


Distance ae! 5 (Ref BT) {ft} 
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TRANSECT #6 

Location: Transect #6 

Date & Time: September 5, 2019 4:50 PM 

Equipment: RiverRay ADCP attached to the side of an inflatable raft with outboard motor 
WSE (ft NAVD88): 8.82 

Discharge ([Q] cfs): 28,900 

Velocity ([v] fps): 2.8 

Measurement Rating: Fair 

Measurement Notes: At the time of the measurement, open-channel conditions were present. Wind 


was negligible and surface waves were not present. Prior to deployment, 
diagnostic tests were performed, and the internal compass calibrated. 


A tributary enters the Colville River at this crossing location. Discharge was 
measured just downstream of the tributary. The deeper channel bathymetry 
influenced by this tributary is evident in the left bank side of the profile in the 
figures below. 


Table 2. Colville River Ocean Point Transect #6 Measured Discharge Summary 


Delta Measured Total Measured Total 


Measurement Starting TotalQ DeltaQ Measured Measured 


: A loci V i 
Transect # Bank (cfs) (%) Q (cfs) epepek ) ‘ Width (ft) Giants (ft2) Taran mae 
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A. Colville River Ocean Point Transect#6: Measurement Transect 000 Raw Data Output 
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4.2. Colville River Ocean Point Measured Water Quality 


Table 3. Measured Water Quality Parameters 


Total = 2b Specific 


Date Depth Depth Temp Conductivity Conductance 


& 


Dissolved Oxygen _ Salinity 
(% 


Time (ft) (ft) (°C)  (uS/cm) 


10 
| 9/05/19 
| 3:50pm 
| 9.8 
pe | a0 


é 4 
9/05/19 
oe 
Notes: ce Ee aS 


(1) sample depth i is ‘measured from the water Sureice Seay 

(2) ‘Temperature, conductivity, and salinity were measured sh: S| Pro 

(3) Dissolved oxygen was measured using a YSI ProDO meter. : Kee Boe ae, 
(4) Specific conductance (Spee ae eno) was > obtained using a Conve Ol coefficient of 0. 01 96 based 
onempirical data. ee Re COI Ran Ge oe CUE ECOG 
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Active Layer Investigation Results 


Transect 1 - Left Bank 
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Approximate locations of active layer probe transects are provided in Figure 2 and Figure 3 of main report. 


Active layer bank OES were er surveyed. 
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Transect 6 - Right Bank 


Distance Depth 


*Thaw depths of 5 feet indicate frozen 
ground was not encountered. 


Transect 6 - Left Bank 
Distance | Depth 


| Edge of water 


170 | 


| tundra near bluff | 


180 


Approximate locations of active layer probe transects are provided in Figure 2 and Figure 3 of main report. 


Active layer bank elevations were not Aue 
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Transect 1 -Left Bank Transect 1 -Right Top of bluff 


Bank 


Top of bluff 


Beach ~ 3.1' to 3.3' 


Beach Edge of water Edge of water 


Small creek draining from tundra 


thaw >5' 


STA 300 STA 0 


Transect 6 -Left Transect 6 -Right 
Top of bluff 
Bank 


Top of bluff, thick willows 


Bluff, sloughing and eroding 
STA 175 


Notes: 
Approximate locations of active layer probe transects are provided in Figure 2 and Figure 3 of main report. 
Active layer bank elevations were not surveyed. 
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Location Name: Colville River at Ocean Point - South (Transect #1) Date Collected: 12/31/2019 
Field Party: C. Lematta, M. Hendee (ICE) Computed By: G. Yager Checked By: H. Runa 
Start Time: 11:20 Finish Time: 14:10 Weather: winds 7mph WNW, Partly Cloudy Temp: =e 
Channel Characteristics: Effective Width: 650 ft Average Velocity: 0.15 fps 

Effective Area: 880 sq ft Discharge: 135 cfs 
Measurement Details: Method: Midsection; 0.6 depth Number of Sections: 12 
Crossing: Wading Cable [Under Ice] Boat Meter: HACH FH950 
Side of bridge: Upstream Downstream N/A N/A ft above bottom of weight 


GAGE READINGS Weight: N/A lbs 
N/A Count: N/A 
Spin Test: N/A revolutions 
after N/A minutes 


Measurement Rated: Excellent | Good Fair based on "Descriptions" 
Descriptions: 


From Field Notes: 


All water columns were less than 2.5 ft. deep. Measurement began on the East/Left Edge of Water (Sta 0+00). Velocity measurements 
were inititally collected with a Sontek acoustic doppler velocity meter but results were inconsistent and unreliable. Measurements were 
then collected with a Hach electromagnetic velocity meter beginning at Sta 5+00 and remeasured at Sta 5+50 and 6+00. Large quantites 
of sediment encountered while drilling through ice results in dulling of all bits beyond ability for use at station 3+50. The thermal drill does 
not create holes large enough for velocity meter probes. Ice and water depths were measured between stations 3+50 and 0+00 and 

velocity at those stations is estimated. 


Calculation Notes: 


The average measured velocity was extrapolated to stations where velocity measurements were not acquired (Sta 3+50 to 0+00). This 
resulted in a computed discharge of 135 cfs. To provide a range of uncertainty, if the minimum velocity was extrapolated to these 
stations, the computed discharge would be 111 cfs, and if the maximum velocity was extrapolated, the computed discharge would be 176 
cfs. 


Colville River at Ocean Point Transect #1 
Date Collected: 12/31/2019 
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Total Discharge: 134.8 cfs 
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Project Name: Willow Ice Road Support Field Personnel: C. Lematta 
Project Number: 174311 Date Collected: December 31, 2019 


Velocity Measurement 


Location: Colville River Ocean Point South (Transect #1) Method: under ice; 0.6 depth 
Weather: -20°F, 7 mph wind Meter: HACH FH950 


: Ice Under Ice Water : 
Station peering! Thickness Depth verouy 

(NAD83) (ft 6 (ft/s) 
OF007 >" East/Lett Bank, N70.05340 W151.36929" [7 15)e5| “grounded [| 
eee I See | a).O. kee | Megrounded 42 |im spe 
ce DE) el Dc © 
LS EN ET Ea ee a cd © = Se 
fe nnn | m0 ay ea |e DID Ween | WK 28 He By Pe OB | 
_ Se ee ae a ee 
A ee ee NO 0ba70 Wb 1 G2 eS eae sgrounded iad eo 
| 7#50__ | West/Right Bank; N70.05319W151.37564 | 0.7 ~~ | grounded |= 


Notes: All water columns were less than 2.5 feet deep. Measurement began on the East/Left Edge of Water 
(Station 0+00). 


Velocity measurements were initially collected with the Sontek acoustic doppler velocity meter, but 
results were inconsistent and unreliable. Measurements were then collected with a Hach 
electromagnetic velocity meter. 


Large quantities of sediment encountered while drilling though ice results in dulling of all bits beyond 
ability for use between station 3+50 and 0+00. Thermal drill does not create holes large enough for 
velcity meter probes. Ice and water depths were measured between stations 3+50 and 0+00 and 


velocity is estimated for each station, indicated by an" * ". 
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GPS coordinates given in NAD83. Mud not measured. 
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Ocean Point Crossing Profile 
December 31, 2019 
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- Mud not measured at this crossing 
- Water elevation assumed 0.0' 
- Stations are based on West Bank as Station 0+00 
Channel Bottom Sap ihe 6 : 
ae ie eee - Stationing is approximate 
- Datum is NAD8&3 


Bottom of Ice 
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Colville River at Ocean Point- Transect #1 Michael Baker 
Water Quality INTERNATIONAL 
Sample Date: December 31, 2019 


Water Ice Specific DO 
Conductance (% 


Freeboard, Sample Depth Temp Conductivity 
(ft) (uS/cm) 


Salinity 


Location & Time Depth Thickness 
(ppt) 


Notes: 
(1) Sample location coordinates referenced to NAD83 datum. 

(2) Freeboard is the distance from the top of ice to the water surface. 

(3) Sample depth is measured from the water surface. 

(4) Temperature, salinity, and conductivity were measured using a YS! Pro1030 meter. 

(5) Specific conductance (referenced to 25°C) was obtained using a conversion coefficient of 0.0196 based on empirical data. 
(6) Dissolved oxygen was measured using a YSI ProODO meter. 

(7) Time shown indicates the start of the measurement. 

(8) Temperature measurements have an accuracy of +/- 0.2°C 
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Michael Baker 
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Discharge Measurement Notes 


Location Name: Colville River at Ocean Point - South (Transect #1) Date Collected: 2/25/2020 
Field Party: D.Roe, J. Varga (UMIAQ) Computed By: D. Roe Checked By: H. Runa 
Start Time: O25 Finish Time: 16:00 Weather: winds 5 mph, Sunny Temp: 4) 4 
Channel Characteristics: Effective Width: 304 ft Average Velocity: 0.04 fps 

Effective Area: 228 sq ft Discharge: 9 cfs 
Measurement Details: Method: Midsection; 0.6 depth Number of Sections: 13 
Crossing: Wading Cable Boat Meter: HACH FH950 
Side of bridge: Upstream Downstream N/A N/A _ ft above bottom of weight 


GAGE READINGS Weight: N/A Ibs 
gee ee 2) CAGE READINGS Sara eee 
Start Finish Change 


Sta 2+85 | 5.48 ft NAVD88 RTK survey Count: N/A 
| Spin Test: N/A revolutions 
| after N/A minutes 
Measurement Rated: Excellent Good Poor based on "Descriptions" 


Descriptions: 


From Field Notes: 


Negative freeboard occurred between Sta 5+12 and 5+72, averaging 0.2' above the top of ice surface. Positive occurred between Sta 2+85 
and 4+52, averaging 0.2' below the ice surface. Ice was grounded in the middle of the channel from Sta 4+72 to 4+92. Depth of ice was not 
recorded at these locations but is estimated at approximately 5.0'. Discharge was measured beginning where water was encountered at Sta 
5+72; observations were spaced every 20' using RTK GPS to where grounded ice was encountered at Sta 2+85. 


Calculation Notes: 
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Colville River at Ocean Point Transect #1 
Date Collected: 02/25/2020 
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Total Discharge: 


9.2 cfs 


Project Name: Willow Ice Road Support Field Personnel: D. Roe 
Project Number: 174311 Date Collected: February 25, 2020 


Velocity Measurement 


Location: Colville River Ocean Point South (Transect #1) Method: under ice; 0.6 depth 
Weather: -5°F, 5 mph wind Meter: HACH FH950 
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West/Right Bank; N70.053197 W151.375716 


Survey: Water surface elevation surveyed at station 2+85 = 5.48 ft NAVD88. 


Notes: All water columns were less than 1.3 feet deep. Measurement began on the East/Left Edge of Water 
(Station 0+00). 


Velocity measurements were attempted using the Sontek acoustic doppler velocity meter, but results 
were inconsistent and unreliable. Measurements were then collected with a Hach electromagnetic 
velocity meter instead. 


Negative freeboard occurred between station 5+12 and 5+72, averaging 0.2' above the top of ice 


surface. Postive freeboard occurred between stations 2+85 and 4+52, averaging 0.2' below the ice 
surface. 
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Grounded 0.3 West Bank; N70.05309 W151.37573; Crossed Blue Tipped Lath Placed Here 


Avg floating ice thickness = 


General Comments: 
GPS coordinates given in NAD83. Mud not measured. Water depth at upstream transducer = 5.6'; Water depth at downstream transducer = 5.6'; the snow birdge has been installed 
between the transducers and upstream of the ice profile survey; RTK survey was conducted on the edge of ice, edge of water, transducer locatoins, and at each velocity measurement 


location. 
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Ocean Point Crossing Profile 
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NOTES: 


- Mud not measured at this crossing 
- Water elevation assumed 0.0' 
- Stations are based on East Bank as Station 0+00 
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West End of Crossing 
N70.05309 W151.37573 


Colville River at Ocean Point- Transect #1 

Water Quality INTERNATIONAL 

Sample Date: February 25, 2020 

Specific DO 

Conductance (% 
(uS/cm) 


Water Ice 
Location & Time Depth Thickness 
(ft) (ft) 


Be 
Notes: 


(1) Sample location coordinates referenced to NAD83 datum. 

(2) Freeboard is the distance from the top of ice to the water surface. 

(3) Sample depth is measured from the water surface. 

(4) Temperature, salinity, and conductivity were measured using a YSI Pro1030 meter. 

(5) Specific conductance (referenced to 25°C) was obtained using a conversion coefficient of 0.0196 based on empirical 
data. 

(6) Dissolved oxygen was measured using a YSI ProODO meter. 

(7) Time shown indicates the start of the measurement. 

(8) Temperature measurements have an accuracy of +/-0.2 degrees Celcius 


Salinity 
(ppt) 


Freeboard Sample Depth Temp Conductivity 
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OCEAN POINT WEATHER STATION 


SYSTEM DESCRIPTION 
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OCEAN POINT WEATHER STATION 


OVERVIEW 


The Ocean Point weather station was installed on 14 April 2020 for the purpose of monitoring weather 
and climate parameters. The station is located 12 miles SW of Nuiqsut, Alaska on the north side of the 
Colville River where it bends to the south at Ocean Point. The vicinity map is shown in Figure 1. The 
geographical coordinates of the weather station are N70.08730, W151.35590. The site elevation is 
approximately 129 feet above mean sea level. The location map is shown in Figure 2. 


The weather station monitors, in real-time, wind speed and direction, solar radiation, air temperature, 
relative humidity, snow depth, barometric pressure, and ground temperature. 
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Figure 1 — Ocean Point Weather Station Vicinity Map 
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Figure 2 — Ocean Point Weather Station Location Map 


The Ocean Point weather station is configured with the following instruments that are depicted in Figure 
3; 


Wind monitor (wind speed and direction) 
Pyranometer (solar radiation) 

Air temperature sensor 

Relative humidity sensor 

Snow depth sensor 


Barometric pressure sensor 


seer latte SIT 


Ground temperature sensor array 


Page 3 


eam sn ctettee i el i ee ee ee ee ee ee a i, ee ee, ee rr a 


ALICE 


DESIGN & CONSULT 


The sensor information is processed by a CR1000X data processor and transmitted to a remote server 
via satellite modem. Power is supplied by a single 12-volt deep cycle battery that is recharged by a 50- 


watt solar panel. 


Lightning Rod 


| 


iridium Antenna 


Datalogger, 
Barometer, Modem 


Pyranometer 


Relative Humidity 
Sensor 


Figure 3 — Ocean Point Weather Station Configuration 
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INSTRUMENTATION SPECIFICATIONS 


Data Processor 
Campbell Scientific Inc. CR1000X 
Serial number: 3236 


Wind Monitor: 
R. M. Young Model 05108-45 Alpine Version Wind Monitor 
Serial number: WM164717 
Operating temperature range: -50°C to +60°C 
Wind speed range: 0 to 224 mph 
Accuracy: +/- 0.6 mph or 1% of reading 
Starting Threshold: 2.2 mph 
Wind Direction Range: 0 to 360° 
Accuracy: +/- 3° 
The wind monitor is mounted at a height of 7 feet above the ground 


Pyranometer 
Hukseflux LP02 
Serial number: 47571 
Light spectrum waveband: 285 to 3000 nm 
Maximum irradiance: 2000 W/m‘42 
Sensitivity (nominal): 15 yV/(W/m‘’2) 
Operating temperature range: -40°C to +80°C 
Temperature dependence: <0.15% per °C 
ISO classification: Second class 


Snow Depth Sensor 
Campbell Scientific, Inc. SR50A Sonic Ranging Sensor 
Serial number: 335867 
Operating temperature range: -45°C to +50°C 
Measurement range: 1.6 ft to 32.8 ft 
Resolution: 0.01 in 
Accuracy: +/- 0.4 in or 0.4% of distance to target (whichever is greater) 
Height above ground (HAG) = 64.25 in 
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Air Temperature / Relative Humidity Sensor 
HygroVUE5 Digital Temperature and Relative Humidity Sensor 
Sensor element: SHT35 
Temperature range: -40°C to +70°C 
Tolerance: +/- 0.4°C (over -40°C to +70°C range) 
Response time: 130s (wind speed of 1 m/s) 
Resolution: 0.001°C 
Long-term drift: < +/- 0.03°C per year 
Humidity range: 0% to 100% RH 
Accuracy (at 25°C): +/- 1.8% (over 0% to 80%), +/- 3% (over 80% to 100%) 
Additional errors (-40°C to +60°C): +/- 1% 
Short-term hysteresis: < 1% 
Response time: 8s (wind speed of 1 m/s at +25°C) 


Barometer 

Setra CS100 

Serial number: 7325117 

Pressure range: 600 mBar to 1100 mBar 

Accuracy: +/- 0.5 mBar (+20°C) 
+/- 1.0 mBar (0°C to 40°C) 
+/- 1.5 mBar (-20°C to +50°C) 
+/- 2.0 mBar (-40°C to +60°C) 

Linearity: +/- 0.4 mBar 

Hysteresis: +/- 0.05 mBar 

Repeatability: +/- 0.03 mBar 


Modem 
Iridium 9522B 


Ground Temperature Array 
BeadedStream Digital Acquisition Cable 
Serial number: 3484 
Operating temperature range: -55°C to +125°C 
Sensor accuracy: +/- 0.1°C from -10°C to +30°C 
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DATA PROCESSING and TRANSMISSION 


The following parameters are sampled at 30 second intervals: 
Wind speed and direction 
Solar radiation 
Air temperature 


The following parameters are sampled at 10 minute intervals: 
Ground temperature 


The following parameters are sampled hourly: 
Snow Depth 


Barometric Pressure 


The parameters are transmitted hourly via the Iridium modem to a host computer operated by Polar 
Alpine Inc. The information is presented on a private web page. 
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APPENDIX A — POWER SUPPLY INFORMATION 
The Ocean Point weather station is powered by the following equipment: 
Sun Xtender® PVX-840T Solar Battery 


Campbell Scientific, Inc. SP50 Solar Panel 
Campbell Scientific, Inc. CH150 Charging Regulator 
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ABBREVIATIONS 


ATM Atmospheric Pressure 
BPMSL British Petroleum Mean Sea Level 
GEDD Cumulative Freezing Degree-Day 
DS Downstream 
DS-2P Drill Site 2P approximately 21.5 miles east of the proposed OPIB 
ft Feet 
ft/sec Feet per second 
ft? Square Feet 
ft?/sec Cubic Feet Per Second 
GMT2 Greater Moose’s Tooth Drill Site 2 
gpm Gallons Per Minute 
ICE Innovative Civil Engineering, Design, and Consult 
MBI Michael Baker International 
NAD83 North American Datum 1983 
OPIB Ocean Point Ice Bridge 
psi pounds per square inch 
PF Pressure Transducer 
Q Discharge (cubic feet per second) 
SPMT Self-Propelled Modular Transporter 
DLA Survey Station (feet) 
¥US Upstream 
WSE Water Surface Elevation 
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DEFINITIONS 

Hagglund A low ground pressure rubber tracked articulated vehicle 

Overflow When water is present on the surface of the ice 

Rolligon An extremely low ground pressure pneumatic rubber tire all-terrain vehicle (as 

low as 4 psi) with an indirect drive that is approved for summer tundra travel 
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1 EXECUTIVE SUMMARY 


A third winter investigation of the proposed Ocean Point Ice Bridge (OPIB) location was 
conducted in the winter of 2020-2021. For the first time a topographic survey along the 
alignment and across the river from bank-to-bank was conducted. With this new information 
the ice bridge design was updated. Prior to the topographic survey, the elevations of the river 
banks were estimated. The updated ice bridge design is included in Appendix D and the 
construction quantities are presented in Table 1.1. The ice bridge is designed to support the 
Willow Development maximum module net weight of 3,200 tons loaded on a Self-Propelled 
Modular Transport (SPMT) with maximum allowable gross weight of 4,200 tons. 


TABLE 1.1: ICE CONSTRUCTION QUANTITIES OCEAN POINT ICE BRIDGE 
Ocean Point Ice Bridge Ice Quantity (yd°) 
**TOTAL = 64,700 


**The water equivalent quantity is 11.7 million gallons 


Observations from the winter of 2020 - 2021 have confirmed that the crossing has the 
potential to naturally ground. Overflow was not observed or reported to have occurred at the 
proposed crossing site. 


Two direct discharge measurements were conducted by Michael Baker International (MBI). 
The first discharge of 13.8 ft?/s (6,200. gpm) was measured on February 17, 2021. This had 
decreased to 0.7 ft/s (300 gpm) by March 10, 2021 (95% decrease). Discharge measurements 
were planned for March 234, April 7, and 21st but they were not possible since the crossing 
had become naturally grounded. However, there was a possible trace amount of discharge 
through the crossing via the channel bottom and small pockets of water below the ice. 


Future data collection efforts should be conducted at Ocean Point to better understand how 
the water discharge varies through the winter and the historical channel morphological data 
record. It is important to verify the frequency and magnitude of overflow events and the 
mechanisms that cause them. Continuing investigations should focus on the time period 
between mid-February to mid-April. This is the forecast timeframe for the construction and 
use of the OPIB (2024-2025) for the Willow module move. 
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2 INTRODUCTION 
The objectives of the 2020-2021 investigation were to: 


1. Determine the ice and hydrological conditions at the proposed OPIB location 
2. Update the OPIB design pending any new discoveries 


The following data was collected for this effort: 


1. Cumulative Freezing Degree-Days (CFDD) 5. Water surface elevation (WSE) 
2. Natural ice growth 6. Water discharge (Q) (MBI 2020) 
3. Span of free water 7. Overflow observations 


4. Local weather observations 


KA Preliminary August 6, 2021 
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3 SUPPOSITIONS 


The following suppositions should be considered while reviewing this report. 


1. The first year SPMTs moving with modules will occur during the winter of 2024-2025. 

2. The module move is scheduled to take place between mid-February and mid-April. 

3. The ice structure designs are subject to revision based on new information. 

4. There is a potential for overflow at the proposed OPIB site (ICE 2020). 

5. The proposed OPIB will not be a totally grounded ice bridge. 

6. Loaded and unloaded SPMTs are capable of negotiating grades up to a maximum 5% 
longitudinal and with a maximum 1.5% lateral slope (side-to-side). 

7. SPMTs are limited to a maximum grade change of +/-1.7% over 96.5 feet as measured 
along the travel surface. 


= Preliminary August 6, 2021 
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4 DATA AND OBSERVATIONS 


An area map of the proposed OPIB location is included in Appendix A. Five field 
investigations were made between February 17 and April 21,2021. MBI engineers 
accompanied ICE engineers for each of the field visits. The seasonal timing of the field 
investigations was selection to coincide with the proposed Willow Module Move schedule. A 
Rolligon and Hagglund were used to access the project site from Drill Site 2P (DS-2P). 


Ice profiles were conducted during each field visit (Appendix B). The following data were 
recorded during each ice profile: 


- -Span of ice - Span of free water under the ice 
- -Ice thickness - Water depth 
- -Snow depth 


The surface of the ice was the basis of elevation for each of the ice profiles instead of British 
Petroleum Mean Sea Level (BPMSL). It is common practice to reference ice profiles to the 
WSE. BPMSL ice elevations are not necessary during the early phases of ice bridge design. 


The data collected from the ice profiles are necessary for calculating the following 
information: 

1. Crossing Cross-sectional area 

2. Construction Quantities 

3. Direct Discharge (Measured by MBI) 


Pressure transducers (PT) were installed approximately 200 ft upstream (US) and 
downstream (DS) of the proposed OPIB centerline alignment (Appendix D). The PTs were 
installed on the channel bottom in the deepest part of the cross-section. The PTs measured 
absolute pressure which was translated into water depth by atmospheric (ATM) pressure 
corrections. The sample frequency of the PTs was set to 15-minute intervals commencing on 
February 17 at noon and terminating on April 21st at noon. 
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4.1 CUMULATIVE FREEZING DEGREE-DAYS 


Cumulative Freezing Degree-Days (CFDD) are calculated as a sum of average daily degrees 
below freezing for a specified time period and are frequently used to measure and compare 
the coldness of winter from year to year. The annual CFDD (referenced to Fahrenheit 
degrees) has ranged from a high of 9,300 (2011-2012 winter season) to a low of 6,200 (2017- 
2018 winter season). The higher the CFDD the colder the winter. However, a higher CFDD 
doesn’t necessarily equate to a longer winter. Figure 4.1 presents the historical CFDD from 
2002 to 2021. The 2020-2021 winter was on the colder end of the data set spectrum. 


NPRA TUNDRA MONITORING STATION 
CUMULATIVE FREEZING DEGREE-DAYS 
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FIGURE 4.1: CUMULATIVE FREEZING DEGREE-DAYS 2002 - 2021 


The air temperature CFDD index correlates with: 
1. Natural ice growth 
2. Ice construction rates 
3. Refreezing of seasonal thawed tundra 
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4.2 SPAN OF ICE 


The span of ice is the surface distance between the outer most edges of ice across a water 
body. The initial length of the span of ice at the proposed OPIB crossing location on 
February 17 was approximately 1,000 ft. It was difficult to establish the precise edges of the 
ice since the interfingering of ice and sediments made it difficult to establish the exact edge 
location. Generally, the river ice was blown clean and snow depths tended to be less than 
Oeit 


4.3 NATURAL ICE GROWTH 


Table 4.1 presents the ice natural ice thickness and ice growth rate at the proposed OPIB site 
over the course of the five field visits. The ice growth rate is not shown after March 234 since 
the crossing became grounded. Natural river ice growth rates depend on air and water 
temperatures, water velocity, overflow, wind speed and direction, snow cover, and ice 
thickness. Historically, natural ice growth after mid-April is relatively low. This is primarily 
due to the increase in average daily temperatures and solar radiation. As a result, the 
reduction in natural ice growth may make ice bridge repairs and construction difficult during 


this timeframe. 


TABLE 4.1: NATURAL ICE THICKNESS AND GROWTH RATE - 2021 


Thickness (feet) Field Visit (feet per day) Field Visit (feet per week) 
We : 
*Day 1 is set to the first day that CFDD > 1 (September 10, 2020) 
tCrossing is grounded 


4.4 SPAN OF FREE WATER 


The span of free water is the distance between the edges of water below the ice. Figure 4.2 
provides an illustration of the span of free water. 


Table 4.2 presents the length of the free water span during each of the five field visits. As the 
ice thickness increases the span of free water length and column thickness decreases. 


By the time of the April 7* field visit, the natural ice had grounded across the alignment. It is 
possible that small pockets of flow were present under the ice but were not detected with the 
50-foot ice profile intervals. 
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Ocean Point Crossing Profile 
February17, 2021 
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FIGURE 4.2: ICE PROFILE AND SPAN OF FREE WATER AREA 
TABLE 4.2: SPAN OF FREE WATER 
’ Me Span Of Free Water Length Maximum Free Water 
Field Visit Date ; 
(feet) Column Height (feet) 
*Single location along ice profile was not grounded; additional ice profiles were 
conducted 200 feet upstream and downstream parallel to the alignment and were 
grounded 
tCrossing is grounded 
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4.5 DISCHARGE MEASUREMENTS 


MBI conducted discharge measurements during the first two field visits (Appendix C). A 
summary of the MBI discharge data are presented in Table 4.3 and are overlaid in Figure 4.3. 


TABLE 4.3: MBI DISCHARGE MEASUREMENTS 


5 Discharge Discharge Average Velocity 
Field Visit Date 3 
(ft'/s) (gpm) (ft/s) 
March 10, 2021 0.7 300 0.01 


The average water velocities and discharge decreased as winter progressed. Unlike the 
Mackenzie and Yukon Rivers, the Colville River is classified as an Arctic River. This means 


winter flow stops since the watershed is frozen. However, some flow continues from 


groundwater and springs. 


4.6 WINTER WATER DEPTHS 


Generally, water depth in a river decreases as the tributary discharge decreases with the onset 
of freeze-up. This was the case with Ocean Point location for the winter of 2020-2021. 
However, this was not the case during the winter of 2019-2020 (ICE 2020). 


Figure 4.3 presents the water levels at the OPIB crossing as measured by the PTs. The water 
levels began to decrease about 4 days after the PT installation on February 17. The US PT 
froze around March 17+ and then thawed around April 13. The data recorded during this 
timeframe is not representative of the water level but the pressure induced onto the PT by the 
expanding freezing water - ice. Similarly, the DS PT froze around April 2™4 and was removed 
from the channel before it had a chance to thaw in-situ. The reason for the differences in 
freeze dates may be attributed to differences in water velocities, relative snow cover, water 
chemistry, and distribution of flow within the braided channel. Overall water depths 
measured by the PTs were less than 5.0 ft with the maximum 5.2 ft recorded during the 
installation of the PTs on February 17%. 


Table 4.4 presents the maximum water depths at the OPIB location measured during the ice 
profiles. The maximum water depths decreased with each successive field visit. Once the ice 
became grounded water depths could not be measured. A study of all the possible water 
variations is beyond this report. 
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2021 Ocean Point Water Level 
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FIGURE 4.3: WATER DEPTHS AT THE PROPOSED OCEAN POINT CROSSING 


TABLE 4.4: MAXIMUM WATER DEPTHS ALONG ICE PROFILE ALIGNMENT - 2021 


: 55: Maximum Water Depth 
Field Visit Date 
(feet) 


= 
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4.7 WATER AND AIR TEMPERATURES 


Water and air temperatures were recorded at Ocean Point. Figure 4.4 presents the river 
bottom water temperatures recorded at 15 minutes intervals. Figure 4.5 presents the average 
daily air temperature as recorded at the DS-2P weather station. The DS-2P weather station is 
located approximately 21.5 miles east of Ocean Point. What is noteworthy within the data set 
are the following: 
1. The below freezing temperatures recorded by the transducers. This is consistent with 
the grounded ice profiles conducted after March 23", 
2. The thermal drill created temperature spikes as a result of hot water being used to drill 
through the ice are evident on Figure 4.3. 
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3. Average daily air temperatures stayed below freezing between February 17% and 
April 21. 


2021 Ocean Point Water Temperature 
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FIGURE 4.4: RIVER WATER AND ICE TEMPERATURES AT OCEAN POINT 


August 6, 2021 
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FIGURE 4.5: AVERAGE DAILY AIR TEMPERATURES DRILL SITE 2P WEATHER STATION 


4.8 OVERFLOW CONSIDERATIONS, CHARACTERISTICS, AND OBSERVATIONS 


Unlike the 3 overflow events during the previous winter (ICE 2020), overflow events were not 
apparent within the PT data record (Figure 4.3). Furthermore, evidence of overflow events 
was not observed during the 5 field visits. 


A unique design consideration for the proposed OPIB is the potential for overflow. There is a 
high likelihood for overflow on the Colville River during any given winter. The frequency, 
location, triggers, and magnitude of overflow events can be difficult to predict and measure. 
Considerable attention was made during the data collection to record overflow events with 
instrumentation and observations. 


Generally, in early in winter (November - December) overflow events are not typical. The 
thin ice (<2 ft) rises and falls with the changes in water levels. Discharge declines as the rains 
transition to snowfall and the watershed freezes. When overflow events do occur, they tend 
to be of lesser magnitude. Normally, when an overflow event occurs it breaches along the 
edges of the river where the ice has not become firmly grounded. 


By the middle of winter (January - February) the potential for overflow events increases in 
frequency and magnitude. The river ice becomes firmly grounded along the edges of the river 
and the cross-sectional area of free water is substantially reduced. This can lead to higher 
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pressures, constrictions, and increased water velocity under the ice. Eventually the pressure 
becomes great enough to form cracks in the ice. The flowing water finds pathways through 
the cracks to the surface of the ice. The overflow from these events tends to flow in all 
directions on the surface of the river ice. These events can be difficult to observe if there is 
snow on the surface of the ice. Furthermore, the presence of snow can increase the amount of 
time required for the overflow to freeze. 


Toward the end of winter (March - April), the potential for overflow is similar to that of the 
middle of winter with the additional contribution from increased solar radiation and warmer 
air temperatures. Overflow from snowmelt tends to result in ponding and minimal flow. 
Generally, there is no potential for runoff during this time period. The Colville River runoff 
normally occurs May. 
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5 OCEAN POINT ICE BRIDGE DESIGN 


The OPIB design and construction quantities have been updated based on the data and 
observations from the winter of 2020-2021 investigation (Section 4) and the topographic 
survey conducted by Umiaq survey in August 2021. Future revisions will be issued as more 
information is obtained about the crossing. The revised ice bridge design (Appendix E) 
includes two main ramps into and out of the river floodplain. 


5.1 ICE BRIDGE CONSTRUCTION QUANTITIES 


Table 5.1 presents a summary of the in-place ice volume of the completed OPIB based on the 
most recent information. The ice quantities represent ice or water that must be mined, hauled, 
and placed for the completion of the ice bridge. The total takes in to account the typical 
natural ice thickness that will be present at the crossing at the commencement of construction. 
The average maximum natural floating ice thickness is expected to be between 4 and 5 ft by 
February 1*. 


TABLE 5.1: SUMMARY OF OPIB ICE QUANTITIES 


Ice Structure 


*TOTAL = 64,700 


*The water equivalent quantity is 11.7 million gallons 


The increase in construction quantities from the previous bridge design is driven primarily by 
the recently surveyed elevations of the tops of the east and west banks. 


5.2 EMERGENCY BYPASS ROAD AND RAMPS 


An emergency bypass road and ramps will be constructed on the downstream side of the ice 
bridge and ramps. These roads will provide emergency access around the SPMTs while they 
are crossing the OPIB. The emergency bypass roads are not designed for supply traffic or ‘go- 
arounds’ while the SPMTs are traversing the OPIB. Specific design details of the emergency 
bypass road will be provided later. 
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6 CONCLUSIONS & RECOMMENDATIONS 


Conclusions and recommendations are subject to change as new information becomes 


available. 


6.1 CONCLUSIONS 


ah 


4. 


The proposed Upstream Site OPIB will be a non-grounded ice bridge with a 
capacity that is suitable for the 3,200-ton module loaded onto a SPMT with a 4,200- 


ton allowable gross weight. 


. There is a high probability that at least one overflow event will occur in the vicinity 


of Ocean Point each winter. 


. More information is needed regarding water discharge and water levels in the 


vicinity of the Ocean Point between the months of February and mid-April. 
The crossing has the potential to naturally ground. 


6.2. RECOMMENDATIONS 


Ls 


MBI should conduct weekly discharge measurements at the OPIB from mid- 
February to mid-April during the winter of 2021-2022. 


. Ice profiles should be conducted every two weeks at the OPIB from mid-February 


to mid-April during the winter of 2021-2022. 


. Delay any geotechnical investigation at the crossing until 2023-2024; this task may 


not be necessary. 


. Continue to install remote monitoring sites at the OPIB to collect water level, air 


temperature, and air pressure throughout the winter of 2021-2022. 


. Expand the topographic survey of the river banks US and DS of the proposed OPIB 


crossing alignment during the late fall of 2022 (early September). 


. Collect bathymetry at the OPIB during the late fall of 2022. 
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Ocean Point Crossing Profile 
March 10, 2021 
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Ocean Point Crossing Profile 
March 23, 2021 
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Ocean Point Crossing Profile 
April 7, 2021 
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ACRONYMS & ABBREVIATIONS 


% sat percent saturation 

Ae: degrees Celsius 

cfs cubic feet per second 

COPA ConocoPhillips Alaska, Inc. 

DO dissolved oxygen 

ft feet 

ft/s feet per second 

ICE ICE Design & Consult 

US/cm microsiemens per centimeter 

mg/L milligrams per liter 

Michael Baker Michael Baker International 

NAVD88& North American Vertical Datum of 1988 
Ocean Point the South transect, Transect #1, the Rolligon crossing, the west crossing, the 
Upstream upstream crossing, Ocean Point South 
Peak Peak Oilfield Services Company 

ppt parts per thousand 

Q discharge 

UMIAQ UMIAQ, LLC 

USGS United States Geological Survey 
Willow Willow Project 


WSE water surface elevation 
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1.0 INTRODUCTION 


Michael Baker International (Michael Baker) collected water resources data for Conoco Phillips Alaska, 
Inc. (COPA) in support of the Willow Project (Willow). The proposed ice road crossing of the Colville River 
was investigated at Ocean Point. During winter of 2021, data was collected in four field events. This 
report summarizes the methods and results of that effort. 


ICE Design & Consult (ICE), Peak, and UMIAQ, LLC (UMIAQ) provided support during the field program 
and contributed to a safe and productive field season. 


2.0 CROSSING LOCATION 


The Ocean Point Upstream transect near Ocean Point was investigated this year and is shown in Figure 
1. This transect was selected based on shallow water depths relative to the other transects investigated 
in previous years. Ocean Point Upstream (also historically referred to as “Transect #1”, the “Rolligon 
crossing”, the “west crossing”, the “upstream crossing”, “Ocean Point South”) is an historic ice and snow 
road crossing location. It was the location of a snow road during the 2018-2019 season and an ice road 
crossing for Cruz Construction during the 2020-21 season. This is the preferred proposed crossing 
location. 


Table 1 provides a summary of dates and data collected at the transect. Table 2 provides a summary of 
measurements collected. 


Table 1: Field Events 


Ocean Point 


Data Collection 
discharge 


Grounded Grounded 


water quality 
Channel Ice Channel Ice 


water surface elevation 
survey 


1. Discharge estimated based on limited field measurements 
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Data Collected 


cf 


. specific conductance microSiemens per centimeter 
Water Quality — 
dissolved oxygen percent saturation 
dissolved oxygen milligrams per liter 


Water Surface : : : 
: water surface elevation | feet North American Vertical Datum of 1988 
Elevation 
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3.0 METHODS 


Field sampling methods were based on United States Geological Survey (USGS 2006a and 2006b) 
methods. Safety precautions were followed using the North Slope Water Resources 2020 Health, Safety, 
and Environment Plan (Michael Baker 2020a) and the 2020-2021 Winter Hydrology Programs — Job 
Safety Analysis (Michael Baker 2020b). 


Measuring discharge under ice cover is subject to limitations not applicable to open water measurements. 
Unlike open water where it is obvious where the edge of water exists, it is not possible to see the extents 
of the cross-sectional area of flow under the ice. Further, it is not possible to profile the entire measurable 
cross-section since velocity measurements are limited to only where holes are drilled through the ice. It 
is assumed that the cross-sectional area is reasonably uniform upstream, downstream, and between 
measurement stations. However, the potential exists for “unseen” grounded or relatively shallow areas 
which would influence measured velocity direction and magnitude if occurring upstream or downstream 
of a measurement station. Grounded areas between measurement stations would reduce the estimated 
cross-sectional area of flow and resulting discharge. 


Four field trips were performed to investigate the trend in discharge and water quality over the course of 
the ice-cover season. The ice-cover season typically initiates with freeze-up in mid-October and ends 
with spring breakup in mid-May. Ice-cover field events were one day apiece. The trips ranged from Mid- 
February to the end of April. 


A one-person Michael Baker field crew conducted all events, supported by an ICE engineer who 
performed crossing bathymetric profiling. UMIAQ and Peak provided transportation to the sampling 
locations and general field support. The sites were accessed by Hagglund and Rolligon. 


Thermal drill probing was performed by ICE to 
identify the extents of under-ice water bounded 
by ice grounded against the channel bed. Water 
measurements were facilitated by mechanically 
drilling through the river ice. Investigation of soils 
or groundwater within the channel bed was not 
performed. Discharge was determined using 
USGS mid-section techniques. Velocity was 
measured using a handheld Hach flow meter. 
This was attached to a fixed rod and lowered to 
0.6 the water depth below the ice. In-situ water 
quality measurements were collected at the 
deepest section. Field crew used a YSI ProSolo 
meter to collect temperature, conductivity, 
salinity, and dissolved oxygen. Measurements — teenie 
were taken at multiple depths throughout the Photo 3.1: Crew setting up GPS at Ocean Point; 
water column, if possible. AEM 


Previously submitted ice cover season field data 
is provided in Attachment A. 
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4.0 RESULTS AND CONCLUSIONS 


A summary of Colville River Ocean Point water resources information collected during the 2020-2021 
winter field season is provided below. No overflow, aufeis, or evidence of any other notable hydraulic 
occurrence was observed at the transect during the ice-cover field events. Discharge decreased as the 
ice-cover season progressed. 


4.1 DISCHARGE MEASUREMENTS 


The first field event occurred on February 17. ICE profiled the crossing and Michael Baker measured 
velocity through 3 holes before both transport vehicles experienced mechanical issues in the cold 
temperatures. The trip was aborted before all data was collected. The average velocity was applied to 
the profile ICE collected to get an estimated discharge of 13.8 cfs. 


The second trip occur on March 10". The field crew was able to collect discharge and water quality 
measurements at the crossing. This included under-ice cross-sectional bathymetric profiles, discharge, 
velocity, water depth, ice thickness, water surface elevation, site conditions related to overflow, and in- 
situ water quality. The channel ice was grounded out in the middle of the channel creating 2 flow paths. 
The discharge was measured at 0.7 cfs. 


On the third and fourth trip, April 8 and April 21°, the channel ice was completely grounded, leaving no 
liquid water. Colville River discharges measured at Ocean Point are provided in Table 3. 


Table 3: Colville River Discharge Summary 


Ocean Point Upstream 


Date ies = | ae ee Mere bee eae 
2/17/2021 EX 
3/10/2021 
520: ees TT 
4/21/2021 
Notes ye a ee ae 


4.2 WATER QUALITY MEASUREMENTS 


Water quality measurements were collected on March 10". Slightly elevated salinity and conductivity 
measurements suggest this location may have had minor coastal influence this year. Though this year’s 
values show an increase from 2019-2020 ice-cover seasons (average conductivity was 257 uS/cm), 
results are at the upper limits of freshwater. This year, the brackish water moved up the Colville River 
sooner than past years. The dissolved oxygen measurement was typical of water bodies under the 
influence of ice cover because ice prevents the introduction and mixing of atmospheric oxygen into the 
water. 
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Table 4: Colville River Ocean Point Water Quality Summary 


Ocean Point Water Quality 


2/17/2021 


3/10/2021 
[Vas ar 


4/8/2021 


4/21/2021 


Notes: 


1. Aborted measurement due to vehicle issues. 


2. Channel ice was grounded out. 
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Michael Baker Discharge Measurement Notes 


INTERNATIONAL 


Field Party: S. Orizotti, Roy Baldwin (UMIAQ) Computed By: D. Roe Checked By: D. Roe 


Start Time: 12:15 Finish Time: 14:30 Weather: winds 10 mph Temp: -34 °F 


Effective Area: 495 sq ft Discharge: 13.8 cfs 


Measurement Details: Method: | _....Widsection; 0.6 depth Number of Sections: 10 


PSS CR ea ao Sapp lela Lie Aa lia asl Ie 


Side of bridge: Upstream Downstream N/A N/A ft above bottom of weight 


GAGE READINGS Weight: N/A _Ibs 


Spin Test: _N/A__ revolutions 


after ___N/A__ minutes 
Measurement Rated: Excellent Good Fair based on "Descriptions" 
Descriptions: 


From Field Notes: 


and 1+50 before vehicle issues halted work. Measurement was aborted and field crew left for safety. Ice design was able to profile the 
Grossing before aborting. The average velocity of 0.03 was applied to the rest of crossing to estimate discharge. Survey was not completed due ' 


vehicle issues. 


Distance 
from 
initial 
point 


Depth | Thickness 


Freeboard 


(ft) 


Colville River at Ocean Point Transect #1 


Measurement 
Depth Below 
Top of Ice 


(ft) 


Date Collected: 02/17/2021 


Effective | Section | Effective 
v1 V2 
(ft) (ft) (fps) | (fps) 


grounded 


grounded 


i) 
oO 
o 


37.5 
45.0 
45.0 
65.0 
pe 


grounded 


VELOCITY 


V3 


(fps) 


Average V 


(fps) 


Average V 
(fps) 


Project Name: Ocean Point Discharge Field Personnel: S. Orizotti 
Project Number: 183014 Date Collected: February 17, 2021 


Velocity Measurement 


Location: Colville River Transect #1 near Ocean Point Method: under ice; 0.6 depth 
Weather: -34°F, 10 mph wind Meter: HACH FH950 
: Ice Under Ice Water , 
: Location : Velocity 
Station (NAD83) Thickness Depth 
(ft) (ft) 


[TOUTS TAREE Se RESID ccc er 
(et ae 
eee ee 
Mn) ee Soe RUST RIN 
ee i eer 
[eS a ee 
Sa ea eee cand 
(a 
ae ree eae ee ee 
eas SRE 
ee ere ee 
| I ae a 
Ue Ge See ts ee 
be oes) Mine daimng 8) aA 8 


Survey: Vehicle issues arose before survey could commence. 


ea 


Notes: Velocities were taken at station 1+00, 1+25 and 1+50. The average velocity was 0.03 ft/sec in that 
section and was applied to the other stations to estimate the discharge. 
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Depth Below Ice Surface (feet) 


3.0 


Left Bank 


Grounded 
Ice 


Colville River at Ocean Point - Transect #1 


Measured 
Velocities Estimated Velocity 
<< <—$__» |< $$ —— 


1+00 2+0 3+00 4+00 5+00 
Station (feet) 


Right Bank 


Grounded Ice 


——= Channel Bed 


FoF ice 
6+00 7+00 


Michael Baker Discharge Measurement Notes 


INTERNATIONAL 


Location NaMe@? nnn Olville River at Ocean Point = Transect ey ae Date Collected: | 3/10/2021. 
IS OEY Fe eccms K: Braun, J. Varga (UMA) snes Computed By: |S. Orizotti Checked By: |. DROS nme 
Start Time: 11:47 Finish Time: 13:00 Weather: mun WINGS 15 Mp, SUNY. asnesese US To esi ra ale 
Channel Characteristics: Effective Width: 118 ft Average Velocity: | 0.01 fps 
Effective Area: 0. 8 ft Discharge: 0.7. cfs 
Measurement Details: Methods eee Midsection; 0.6 depth. Number of Sections: | 6. 
Crossing: Wading Cable Boat Meter: HACH FH95O 
Side of bridge: Upstream Downstream NAS eee ee N/A _ ft above bottom of weight 
Weight: N/A Ibs 
Sta 1+00 4.18 ft NAVD88 é RTK survey SL AR oA, Seth cote baa a 
SPU test te NA revolutions 
oleiiae 2. NUN Ge minutes 
Measurement Rated: Excellent Good Fair based on "Descriptions" 
Descriptions: 
From Field Notes: _Ice grounded out from stations 0+00 to 1+00 and 3+00 to 7+00. Ice was not grounded out at station 2*50. 


Colville River at Ocean Point Transect #1 
Date Collected: 03/10/2021 


VELOCITY 


‘ 
(ft) (fps) | (fps) | (fps) 


grounded 


fea = 


Distance 
from 
initial 


Adjusted {Discharge 
Average V 


Average 
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(fps) (fps) 


Pema | = 
2+25 4.8 grounded 
Pew [ee[ ee [= | s [@ [ wo] [ow fon pou] om [om ] | 
| aes | | grounded 


Del 5.1 grounded 


i 
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oO 


grounded 
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grounded 
grounded 


grounded 


NO 
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—s 
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i i, ee ee, a ee i ee em ee ee ce ee ee ee ee ee 
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Project Name: Ocean Point Discharge Field Personnel: K. Braun 
Project Number: 183014 Date Collected: March 10, 2021 


Velocity Measurement 


Location: Colville River Transect #1 near Ocean Point Method: under ice; 0.6 depth 
Weather: 20°F, 15 mph wind Meter: HACH FH950 
eeotion Ice Under Ice Water Velocity 
Station (NADB3) Thickness Depth 
(ft) (ft) 
ee ee ee ce nS kr grounded «| 
ee re eaten a Tene OES eos ee |e ees |aetgrounded — 
er eee eer De Se OLB ee 
en mn REI tat eee es ale warn et | AON Pee 0.8. 8 
ee eee rr ee ee OO ee Oe 
ee Ce Ae ee Meare ne a ES Ph A Sm 0.4 
ee ee he Oe A Oe 
tea an | Semen ase Matalin WE | OS ONE erk grounded | 
a a SB grounded or 
eee ee ne ANNE Witare aime | 2 A_B |W grounded sit. pun i a P 
ee ee ee OL Or DOR 
2 Se Saas aaa a ee 
ee ee ee. ae groundeds) of ai en 
er ee ee ae remeber MORRO ees | ar grounded ey fice | 
eo grounded) [| 
Ea a ee 
i ee ee grounded | fy 
ee eee man maison Ne on Wits | ear 4 27) Pet | grounded af se | 
ee ee ee groimded pc ek | 
es es | MEG waa eu ai eves) Pent 2.2 oes) grounded ven} Oe 
ere | east ot Bank N70:053813.W151,969598 i) Tt | grounded’. fe | 


Survey: Water surface elevation surveyed at station 1+00 = 4.18 ft NAVD88. Average ice elevation was 5.4 ft 
and the average water surface elevation was 4.1 ft. 


Notes: All water columns were less than 0.8 feet deep. Velocity measurements were collected with a Hach 


electromagnetic velocity meter. Postive freeboard occurred between stations 1+00 through 1+75 and 
at 2+50, averaging 1.3' below the ice surface. Water quality parameters were taken at Sta 1+00. 
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Depth Below Ice Surface (feet) 


=e 0 ie 


| 
Left Bank 
2.0 4 

) 

i 

| 


Grounded 


Ice 
Saeed 


0.0 +: 


Ne 


Measured 
Velocities 


Colville River at Ocean Point - Transect #1 


Grounded 
Ice 


Measured Grounded Ice 


2+00 3+00 4+00 5+00 
Station (feet) 


Right Bank 


—— Channel Bed 


SEES nNRORRRERORRERA DR SERRERIRRERRREEEERREEREEEEEEEEEEEEEEEE | 


6+00 7+00 


Colville River at Ocean Point- Transect #1 Michael Baker 
Water Quality INTERNATIONAL 
Sample Date: | March 10, 2021 
Water Ice Specific DO 


Freeboard S le Depth T ductivi ini 
Location & Time Depth Thickness fraceaimes ek eteatine aid Drea mee Conductance a! (% ny 


(°C) (uS/cm) ImB/L) 7 coturation) 


30.4 


Notes: 

(1) Sample location coordinates referenced to NAD83 datum. 

(2) Freeboard is the distance from the top of ice to the water surface. 

(3) Sample depth is measured from the water surface. 

(4) Temperature, salinity, dissolved oxygen, and conductivity were measured using a YSI ProSolo meter. 

(5) Specific conductance (referenced to 25°C) was obtained using a conversion coefficient of 0.0196 based on empirical data. 
(6) Time shown indicates the start of the measurement. 

(7) Temperature measurements have an accuracy of +/- 0.2°C 
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Willow Development — Ocean Point Ice Bridge Revision 


Appendix D 
- OCEAN POINT ICE BRIDGE DESIGN DRAWINGS 
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1.0 INTRODUCTION 


Michael Baker International (Michael Baker) collected water resources data for Conoco Phillips Alaska, 
Inc. (CPAI) in support of the Willow Project (Willow). The proposed ice road crossing of the Colville River 
was investigated at Ocean Point. During winter of 2022, data was collected in six field events. This report 
summarizes the methods and results of that effort. 


ICE Design & Consult (ICE), Northern Energy Services (NES), and UMIAQ, LLC (UMIAQ) provided 
support during the field program and contributed to a safe and productive field season. 


2.0 CROSSING LOCATION 


The Ocean Point Upstream transect near Ocean Point was investigated this year and is shown in Figure 
1. This transect was selected based on shallow water depths relative to the other transects investigated 
in previous years. Ocean Point Upstream (also historically referred to as “Transect #1”, the “Rolligon 
crossing’, the “west crossing”, the “upstream crossing”, “Ocean Point South”) is an historic ice and snow 
road crossing location. It was the location of a snow road during the 2018-2019 season, an ice road 
crossing for Cruz Construction during the 2020-21 season and a snow road crossing for Cruz 
Construction during 2021-22. This is the preferred proposed crossing location. 


Table 1 provides a summary of dates and data collected at the transect. Table 2 provides a summary of 
measurements collected. 


Table 1: Field Events 


Data Collection 
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Table 2: Data Collected 


Data Collected eee Unitsiaae 
water depth, under ice 


Data Type 


ice thickness 
ft 
velocity 


; specific conductance 
Water Quality 


salinity ppt 


dissolved oxygen % sat 
dissolved oxygen 


ft NAVD88 


Water Surface 
Elevation 


water surface elevation 
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Figure 1: Colville River Ocean Point Upstream Transect 
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3.0 METHODS 


Field sampling methods were based on United States Geological Survey (USGS 2006a and 2006b) 
methods. Safety precautions were followed using the North Slope Water Resources 2021 Health, Safety, 
and Environment Plan (Michael Baker 2021a) and the 2021-2022 Winter Hydrology Programs — Job 
Safety Analysis (Michael Baker 2021b). 


Measuring discharge under ice cover is subject to limitations not applicable to open water measurements. 
Unlike open water where it is obvious where the edge of water exists, it is not possible to see the extents 
of the cross-sectional area of flow under the ice. Further, it is not possible to profile the entire measurable 
cross-section since velocity measurements are limited to only where holes are drilled through the ice. It 
is assumed that the cross-sectional area is reasonably uniform upstream, downstream, and between 
measurement stations. However, the potential exists for “unseen” grounded or relatively shallow areas 
which would influence measured velocity direction and magnitude if occurring upstream or downstream 
of a measurement station. Grounded areas between measurement stations would reduce the estimated 
cross-sectional area of flow and resulting discharge. 


Six field trips were performed to investigate the trend in discharge and water quality over the course of 
the ice-cover season. The ice-cover season typically initiates with freeze-up in mid-October and ends 
with spring breakup in mid-May. Ice-cover field events were one day apiece. The trips ranged from Mid- 
February to the end of April. 


A one-person Michael Baker field crew conducted all events, supported by an ICE engineer who 
performed crossing bathymetric profiling. UMIAQ and NES provided transportation to the sampling 
locations and general field support. The sites were accessed by Hagglund and Rolligon. 


Thermal drill probing was performed by ICE Photo 3.2: Crew measuring ice thickness; 3/15/22 
to identify the extents of under-ice water 


bounded by ice grounded against the 
channel bed. Water measurements were 
facilitated by mechanically drilling through 
the river ice. Investigation of soils or 
groundwater within the channel bed was not 
performed. Discharge was _ determined 
using USGS _ mid-section techniques. 
Velocity was measured using a handheld 
Hach flow meter. The meter was attached to 
a fixed rod and lowered to 60% of the water 
depth below the ice. In-situ water quality 
measurements were collected at the 
deepest section. Field crew used a YSI 
ProSolo meter to collect temperature, 
conductivity, salinity, and dissolved oxygen. 
Measurements were taken at multiple 
depths throughout the water column, if 
possible. 


Previously submitted ice cover season field data is provided in Attachment A including ICE profiles of the 
crossing for all six trips. 
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4.0 RESULTS AND CONCLUSIONS 


A summary of Colville River water resources information collected at Ocean Point during the 2021-2022 
winter field season is provided below. No overflow, aufeis, or evidence of any other notable hydraulic 
occurrence was observed at the transect during the ice-cover field events. Discharge decreased as the 
ice-cover season progressed. 


4.1 DISCHARGE MEASUREMENTS 


The first field event occurred on February 16. ICE profiled the crossing and Michael Baker measured 
velocity through 4 holes before the velocity meter battery malfunctioned due to the cold temperatures. 
Two main flow paths were identified at the crossing: the main channel and the western path. Freeboard 
in the western path (Sta. 0+50 to 1+00) and the main channel were 1.0 foot and 0.0 feet, respectively. 
The western path was likely grounded out upstream and water in this channel is backwater from 
downstream where it connects to the main flow path. To estimate the discharge, the main flow path 
velocities were applied to the missing velocity measurements with the same general depth range. The 
estimated discharge was 9.8 cubic-feet-per-second (cfs). The crew collected in-situ water quality 
measurements at the deepest location. 


The second trip occurred on March 3. The field crew collected discharge and water quality measurements 
at the crossing. This included under-ice cross-sectional bathymetric profiles, velocity, water depth, ice 
thickness, water surface elevation, site conditions related to overflow, and in-situ water quality. The 
channel ice was grounded in the middle of the channel creating 2 flow paths. The western flow path (3+30 
to 3+60) was deep enough to collect flow measurements, but the eastern path (4+40 to 4+50) was too 
shallow for the velocity meter to fit under the ice. Discharge was measured at 0.04 cfs. 


The third trip occurred on March 15. The field crew collected discharge and water quality measurements 
at the crossing. This included under-ice cross-sectional bathymetric profiles, velocity, water depth, ice 
thickness, water surface elevation, site conditions related to overflow, and in-situ water quality. The 
eastern flow path had frozen in and the western flow path was flowing water. Velocities had increased 
slightly from the previous trip. Discharge was measured at 0.20 cfs. 


The last three trips occurred on March 30, April 12 and April 27. The channel ice was completely 
grounded, leaving no liquid water. A summary of Colville River discharges measured at Ocean Point are 
provided in Table 3. 
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Table 3: Colville River Discharge Summary 


Average 
Average Ice Water Effective | Average Measured 
Thickness Depth Width Velocity Discharge 
Under Ice (ft/s) (cfs) 


Notes: ND (no data). 

@ Channel was grounded by ice. 

> Velocity meter malfunctioned and only a partial measurement was collected. Velocities were applied to the profile 
to estimate the discharge. 


4.2 WATER QUALITY MEASUREMENTS 


Water quality measurements were collected on February 16, March 3, and March 15. Slightly elevated 
salinity and conductivity measurements suggest this location may have had minor coastal influence 
similar to last year. The dissolved oxygen measurements were low which is typical because the ice cover 
prevents the introduction and mixing of atmospheric oxygen into the water. Table 4 presents the results 
of the water quality measurements on Colville River at Ocean Point. 
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Table 4: Colville River Ocean Point Water Quality Summary 


total depth (ft) 
temperature (°C) 
conductivity (uS/cm) 
dissolved oxygen 
dissolved oxygen (%) 
salinity (ppt) 


specific conductance 


sOe1 


4/21/2021 


2/16/2022 
ND? | 


3: 
ND 
D 
ay, 
2 
3 
ND 
ND 


Notes: ND (no data). 
* Channel ice was grounded. 
> Field crew had to abandon the trip before water quality measurements were obtained. 
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Michael Baker Discharge Measurement Notes 


INTERNATIONAL 


Location Name: Colville River at Ocean Point - Transect #1 Date Collected: = 2/16/2022 
Field Party: Checked By: DIRS coum 
Start Time: 14:50 Finish Time: 16:00 Weather: winds 15 mph, Cloudy Temp: 2a) ls 


Measurement Details: 


Crossing: Wading Cable Boat Meter: = HACHFH9SO 
Side of bridge: Upstream Downstream NA er rey 8 Tyee N/A __ ft above bottom of weight 
GAGE READINGS Weight! NIAW eye Ibs 
Sta 3+50 4.6 ft NAVD88 OUND il tbh oi he NIAT Tons sa 
SpinTest: NIAR ee revolutions 
ACh ae NAC SI minutes 
Measurement Rated: Excellent Good Fair based on "Descriptions" 


Descriptions: 


From Field Notes: Ice grounded out from stations 0+00 to 0+50 and 6+00 to 7+50. 


Colville River at Ocean Point Transect #1 
Date Collected: 02/16/2022 


Distance VELOCITY 


from Measurement Effective | Section | Effective 
5 j Discharge 
mikial ey Dept Thickness) hi corcs aloe oot ce Depth mide pikes ial | Saale 
2 of Ice Average V 
(ft) (ft) (ft) (ft) (ft) (fps) | (fps) | (fps) (fps) (fps) (cfs) 
grounded 


Too shallow for the velocity meter. 


es 


4+00 


Pee 
Ge RS SS 
5+00 ee Too shallow for the velocity meter. 
Sa Too shallow for the velocity meter. 


grounded 
Total Discharge: 


PIEPER ERR RERBE EER EERE EERE ERE REREZEEREREREEGS ESE & | 


Elevation (ft NAVD88) 


Right Edge of Ice Colville River at Ocean Point - Transect #1 


8.0 { 


Left Edge of Ice 


\ 
\ 
\ 
Grounded| Measured 
ice Velocity oe G Tet 
P > ‘ > nabie to Measured rounde ce 
\ Measure Velocity La 


ome Channel Bed 


C= Ice 


— — —Snow 


i Rios SENS aera oie ine Reale shee: mae 


0+00 1+00 2+00 3+00 4+00 5+00 6+00 


7+00 
Station (feet) 


Colville River at Ocean Point- Transect #1 Michael Baker 
Water Quality INTERNATIONAL 


Sample Date: February 16, 2022 


Water Ice Specific 


Freeboard S le Depth Temp Conductivi DO ini 
Location &Time Depth Thickness agra es ceria acid - J Conductance Sines, 


le) 
(ft) (°C) (uS/cm) =a (mg/L) (% Saturation) (ppt) 


Notes: 

(1) Sample location coordinates referenced to NAD83 datum. 

(2) Freeboard is the distance from the top of ice to the water surface. 

(3) Sample depth is measured from the water surface. 

(4) Temperature, salinity, dissolved oxygen, and conductivity were measured using a YSI ProSolo meter. 

(5) Specific conductance (referenced to 25°C) was obtained using a conversion coefficient of 0.0196 based on empirical data. 
(6) Time shown indicates the start of the measurement. 

(7) Temperature measurements have an accuracy of +/- 0.2°C 


Michael Baker Discharge Measurement Notes 


INTERNATIONAL 


Location Names mn Colville River at Ocean Point: Transect Fay mun Date Collected: 9/3/2022 
Field Party: _.D. Roe, M. Rourick (UMIAQ) Computed By: M. Wharton Checked By: KDB te 
Start Time: 13:35, Finish Time: 14:00 Weather uma WiINdS O MPH, SUNY nnn LN otic 
Channel Characteristics: Effective Width: _—ss 20 ft Average Velocity: _ 0.01 fps 
Eifective Area: 0 ie! 6 sq ft Discharge: _0.04 cfs 
Measurement Details: Method: Midsection; 0.6 depth un Number of Sections: 2. 
Crossing: Wading Cable Boat Meter, HACH FHS50_ 
Side of bridge: Upstream Downstream NAAT aaa ei N/A _ ft above bottom of weight 
Weight: N/A___ibs 
Sta 3+50 __ 3.7 ft NAVD88 - COURIC aut ire Se NVA ok nee let of 
| SPIE! GStaas mam INIAcoo tee revolutions 
SU Bialea adn NVA ey minutes 
Measurement Rated: Excellent Good Fair based on "Descriptions" 
Descriptions: 
From Field Notes: ice grounded out from stations 0+00 to 3+30, 4+00 to 4+30, and 4+60 to 7+50. 


Calculation Notes: Average velocity adjusted by coefficient of 0.92 to account for measurements collected at 0.6 feet of depth 


Colville River at Ocean Point Transect #1 
Date Collected: 03/02/2022 


Distance VELOCITY 
from 


initial 


Measurement 
Depth Below 
Top of Ice 


Effective | Section | Effective 


Freeboard Adjusted Discharge 


Average V 


(ft) (fps) | (fps) | (fps) (cfs) 


grounded 


grounded 
grounded 
grounded 
grounded 
grounded 
grounded 


slush 


oT an 


Total Discharge: 


Elevation (ft NAVD88) 


Colville River at Ocean Point - Transect #1 
Right Edge of Ice 


8.0 - 
| 
Measured 
Velocit 
\ Grounded Ice y Grounded Ice 
SSS a a ee rT al 
aN 
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6.0 4 \ 
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Station (feet) 


Left Edge of Ice 


een Channel Bed 
= ce 


— — — Snow 


7+00 


Colville River at Ocean Point- Transect #1 Michael Baker 
Water Quality INTERNATIONAL 
Sample Date: March 3, 2022 
Specific 
Conductance 


Water Ice 


DO DO Salinit 
Location & Time Depth Thickness coer 


Freeboard Sample Depth Temp Conductivity 


(ft) (ft) (°C) (uS/cm) = (mg/L) (% Saturation) (ppt) 


Notes: 

(1) Sample location coordinates referenced to NAD83 datum. 

(2) Freeboard is the distance from the top of ice to the water surface. 

(3) Sample depth is measured from the water surface. 

(4) Temperature, salinity, dissolved oxygen, and conductivity were measured using a YSI ProSolo meter. 

(5) Specific conductance (referenced to 25°C) was obtained using a conversion coefficient of 0.0196 based on empirical data. 
(6) Time shown indicates the start of the measurement. 

(7) Temperature measurements have an accuracy of +/- 0.2°C 


Michael Baker Discharge Measurement Notes 


INTERNATIONAL 


Location Names mun Colville River at Ocean Point = Transect Fy ma Date Collected: 3/15/2022 
eg ee eM Whatton Computed By: Sho Checked By: | DIRi 
Start Time: 13:35, Finish Time: 14:00) Weather mmm WiNdS O MPA, SUNY mm pemnps FetO ak 

Channel Characteristics: Effective Width: ss. 20 ft Average Velocity: 0.03 fps 
ENectvesArca, eas! 6 sq ft Discharge: 0.20 cfs 
Measurement Details: Method. ae Midsection; 0:6 dept nu Number of Sections: 
Crossing: Wading Cable Boat Meters Soe HACH FHOSO 
Side of bridge: Upstream Downstream SUBSE AW elses WO eee N /A Miu ft above bottom of weight 
Weight N/A__ibs 
COME aia Asean. UD Senet 0 Ce 
Spin Lestic(— N/A revolutions 
after ete. Nae ls minutes 
Measurement Rated: Excellent Good Fair based on "Descriptions" 
Descriptions: 
From Field Notes: Ice grounded out from stations 0+00 to 3+30 and 3+60 to 7+00. 


Colville River at Ocean Point Transect #1 
Date Collected: 03/15/2022 
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Colville River at Ocean Point- Transect #1 Michael Baker 
Water Quality INTERNATIONAL 
Sample Date: March 15, 2022 


Water Ice Specific 


Freeboard S le Depth T Conductivi DO DO Salini 
Location & Time Depth Thickness drs cog ripe et rec Porcoptel arta be’ Conductance pe 


(°C) (uS/cm) (mg/L) (% Saturation) (ppt) 


Notes: 

(1) Sample location coordinates referenced to NAD83 datum. 

(2) Freeboard is the distance from the top of ice to the water surface. 

(3) Sample depth is measured from the water surface. 

(4) Temperature, salinity, dissolved oxygen, and conductivity were measured using a YSI ProSolo meter. 

(5) Specific conductance (referenced to 25°C) was obtained using a conversion coefficient of 0.0196 based on empirical data. 
(6) Time shown indicates the start of the measurement. 

(7) Temperature measurements have an accuracy of +/- 0.2°C 
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NOTES: 
- Mud not measured at this crossing 
- Water elevation assumed 0.0' 


Stations are based on West Bank as Station 0+00 


- West edge of water is STA -0+75 
- Stationing is approximate 
- Datum is WGS84 
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1.0 INUPIAQ AND SCIENTIFIC NAMES 

Some readers may better recognize locations, and common plant and animal names by their Ifupiaq or scientific 
names. The appendix provides Ifupiaq names for places (Table E.1.1), and Iftupiaq and scientific names for plants 
(Table E.1.2), mammals (Table E.1.3), fish (Table E.1.4), and birds (Table E.1.5). If an Ifupiaq name did not 
have a known scientific name, it was labeled as unknown (UNK), and vice versa. Figure E.1.1 shows locations of 
the Imupiag place names. 


Table E.1.1. Place Names 


Location 

Site near the mouth of the Miluveach River 
Anaktuvuk Pass 

Bering Sea 

Arctic foothills 

Colville River 

Kuparuk oil field 


Bering Sea-mi Tagiuq 
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Kuukpaaarugmi niuqtuaviq 
Kuukpaaarugmi qimiqqat 
Kuparuk Pingo 

Kuparuk River 

East Channel of the Colville River 
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Kuukpigruaq Kupigruak Channel 
Milugiak 


Napasalu 


Miluveach River and surrounding area 


Channel connecting Nigliq Channel to the Colville River 


Nigligat ‘Second Nuiqsut’, located on the East Channel of the Colville River, near the mouth 
gus of the Colville River 

Nigliq Channel - Westernmost channel of the Colville River Delta, where Nuiqsut is 
located 


Old village site on Nuekshat Island in the East Channel of the Colville River 


Oliktok Point 


Site on a large island in the East Channel of the Colville River, between the mouths 
of the Miluveach and Kachemach rivers; frequently used for caribou hunting 


Nialiq Channel 


Nuiqsapiaq 


Uuliktug nuvugak 
Pisiktagvik 

Qakimak Kachemach River and surrounding area 
Tagium Sifaa Beaufort Sea-mi Beaufort Sea coast 


Tasiqpak Narvaq Teshekpuk Lake 
Source: (HDR 2015; NSB 2016a, 2016b; OHA 2016; SRB&A 2014, 2016; USACE 2012) 
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Existing Pipeline 

§@H Existing Infrastructure 
NPR-A Special Areas 

"ZZ Colville River Special Area 


Teshekpuk Lake Special Area 


Land Designation 
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National Petroleum Reserve 
in Alaska 


No warranty is made by the Bureau of Land 
Management as to the accuracy, reliability, or 
completeness of these data for individual or 
aggregate use with other data. Original data 

were compiled from various sources. This 
information may not meet National Map Accuracy 
Standards. This product was developed through 
digital means and may be updated without 
notification 


Figure E.1.1 
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Table E.1.2. Plants 


Common Name 


Scientific Name 


Draba micropetala Alpine draba 
Draba pauciflora Fewflower draba 


Tussock cottongrass 
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Gimmie, Ur mie sem | 
See eae Vaccinium vitis-idaea 


_ Note: spp. (species); UNK (unknown) 
Source: MacLean 2014 
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Table E.1.3. Terrestrial and Marine Mammals 


Scientific Name Common Name 
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Scientific Name Common Name 


Enhydra lutris kenyoni 
Erignathus barbatus 
Eschrichtius robustus 
Eubalaena japonica 
Eumetopias jubatus 
Histriophoca fasciata 
Lemmus trimucronatus 
Lagenorhynchus obliquidens 
Lepus americanus 
Megaptera novaeangliae 
Mesoplodon stejnegeri 
Microtus miurus 
Microtus oeconomus 
Monodon monoceros 
Mustela erminea 
Mustela nivalis 
Odobenus rosmarus divergens 
Ondatra zibethicus 
Ovibos moschatus 
Phoca hispida, Pusa hispida 
Phoca largha pallas 
Phocoena phocoena 
Phocoenoides dalli 
Physeter macrocephalus 
Sorex tundrensis 
Spermophilus parryii Arctic ground squirrel 

Ursus maritimus 
Kayugtug, Qiangaq, Qigfiqtag Vulpes vulpes 


K Ziphius cavirostris Cuvier’s beaked whale 


Note: UNK (unknown) 
Source: MacLean 2014 
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Table E.1.4. Fish 
Inupiaq Name Scientific Name Common Name 


ennai 
Coregonus laurettae 
Coregonus pidschian 
Coregonus sardinella 


Cottus cognatus Slimy sculpin 
Dallia pectoralis Alaska blackfish 


Eleginus gracilis 
ilag Kakilasak, Kakalisauraq Gasterosteus aculeatus 
Lethenteron camtschaticum 
Liopsetta glacialis 
Panmaksraq, Panmagrak, Panmagraq Mallotus villosus 
Kanayuq Myoxocephalus quadricornis 
Oncorhynchus gorbuscha 
Oncorhynchus keta 


Oncorhynchus kisutch Coho salmon 
Red salmon (sockeye) 


Iqalugaq 


Milugiaq 
Qaaktaq 
Tiipugq 
Aanaaktiq 
Pikuktuug 


Iqalusaaq 


Kanayuq 


Ttuuqiniq 
Uugaq 
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Oncorhynchus nerka 


: Oncorhynchus tshawytscha King salmon (Chinook) 


Osmerus mordax Rainbow smelt 
Saviguunnaq Prosopium cylindraceum Round whitefish 
Kakalisauraq Pungitius pungitius Ninespine stickleback 


Iqalukpik, Paiktuk, Anayuqaksraq, 


alaioik Salvelinus alpinus 
alli 
Paaer Galak 


Note: UNK (unknown) 
Source: MacLean 2014 
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Table E.1.5. Birds 
Inupiaq Name Scientific Name 


Common Name 


Saqsakiq Acanthis flammea and A. hornemanni 


Kurugaq Anas acuta 

Kurugagnaq 

Qaqlutuug, Alluutaq 

Qainniq 

Kurugaqtaq 

Niglivik, Niglivialuk Anser albifrons 
Antigone candensis 
Arenaria interpres 


Aythya valisineria Canvasback 


Tatirgaq 
Tinmiaqpak 
Tullignaq 
Nipailuktaq 


Bartramia longicauda Upland sandpiper 


Branta canadensis 
Bubo scandiacus 


Niglingaq 


Calcarius lapponicus 
Calidris bairdii 


Qupatuk, Putukiutuk 
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Calidris canutus 
Calidris fuscicollis 
Calidris himantopus 
Calidris mauri 
Calidris melanotos 
Calidris minutilla 
Calidris pusilla 


Suyukpaligauraq 
Suigukpaligauraq 


Sugukpaligauraq 
Puvviaqtuuq 
Livilivillauraq 


Livilivillakpak 


Catharus minimus 
Cepphus grylle 
Charadrius semipalmatus 
Chen caerulescens 
Clangula hyemalis 


w 


Kurraquraq 


Papiktuug 
Aaghaaliq 


Qugruk Cygnus columbianus 
Kirgaviatchauraq Falco columbarius 
Kirgavik Falco peregrinus tundrius 
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Inupiaq Name Scientific Name Common Name 


Falco rusticolus 
Gallinago delicata 
Haliaeetus leucocephalus 


Aatqarruaq 


Tuullik 
Taasiniq 
Malgi 


srauq 


Tinmiaqpak 
Agargiq, Nasaullik 
Niksaaktuniq 


Nauyavaaq 


pet) 
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penn 
Larus glaucescens 
Larus hyperboreus 


Limnodromus scolopaceus Long-billed dowitcher 


Nauyavasrugruk 


Sigukpalik 


Limosa lapponica Bar-tailed godwit 


Luscinia svecica 
Melanitta americana 
Melanitta perspicillata 
Melospiza lincolnii 
Motacilla tschutschensis 
Numenius phaeopus 
Passerculus sandwichensis 
Passerella iliaca 


Phalaropus fulicarius Red phalarope 
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Iktigvik 


Phalaropus lobatus Red-necked phalarope 


Phylloscopus borealis Arctic warbler 


Plectrophenax nivalis 
Pluvialis dominica 
Pluvialis squatarola 
Podiceps grisegena 
Polysticta stelleri 
Rissa tridactyla Black-legged kittiwake 

Somateria fischeri Spectacled eider 

Somateria mollissima 


Qinalik Somateria spectabilis King eider 
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Qavaasuk 


Amauligruaq 


Misapsaq Spizella arborea American tree sparrow 
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Stercorarius longicaudus Long-tailed jaeger 


Migiaqsaayuk Stercorarius parasiticus Parasitic jaeger 
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Inupiaq Name Scientific Name Common Name 


Mtg 
Iqirgagiaq Sabine’s gull 
Nunaktuagruk Zonotrichia leucophrys White-crowned sparrow 


Note: UNK (unknown) 
Source: MacLean 2014 
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Glossary Terms 
Active Layer — The top layer of ground subject to annual thawing and freezing in areas underlain by permafrost. 


Albedo — A measure of how a surface reflects incoming radiation; a surface with a higher albedo reflects more 
radiation than a surface with lower albedo. 


Anthropogenic — Resulting from the influence of human beings on nature. 


Black Carbon — A component of fine particulate matter that is formed from the incomplete combustion of fossil 
fuels and biomass. 


Carbon Dioxide Equivalent (CO2e) — The amount of greenhouse gases that would have an equivalent global 
warming potential as carbon dioxide when measured over a specific timescale. 


Greenhouse Gas (GHG) — Gaseous compounds, such as carbon dioxide, methane, and nitrous oxide, among 
others, that block heat from escaping to space and warm the Earth’s atmosphere. 


Lake Tapping — The sudden drainage of lakes caused by ice melting or dislodging and opening up a drainage 
channel. 


Particulate Matter 2.5 (PM2.s) — Particulate matter less than 2.5 microns in aerodynamic diameter in ambient air; 
this fraction of particulate matter penetrates most deeply into the lungs. 


Positive Forcing — When earth receives more incoming energy from sunlight than it radiates to space. 


Thermokarst — A land surface with karst-like features and hollows produced by melting ice-rich soil or 
permafrost. 
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1.0 AFFECTED ENVIRONMENT 


Climate change is affecting natural systems across the globe with enhanced impacts in the Arctic. The atmosphere 
and oceans have warmed, ice cover is shrinking, and permafrost is melting in high-latitude and high-elevation 
regions. The dominant cause of the observed warming since the mid-twentieth century can be attributed to human 
influences (IPCC 2014, 2021). 
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_ and warm the Earth’s atmosphere. GHGs are necessary for keeping the planet at a habitable temperature, and 

_ without GHGs, Earth’s surface temperature would be around 60 degrees Fahrenheit (°F) cooler than it is now. 

_ Natural biological and geological processes regulate levels of naturally occurring GHGs in the atmosphere; 

_ however, anthropogenic emissions haven driven atmospheric concentrations of GHGs to levels unprecedented in 
_ at least the last 800,000 years (IPCC 2014, 2021). Concentrations of CO2, N20, and CHa, have increased by 47%, 
_ 156%, and 23%, respectively, since 1750, largely due to economic and population growth (IPCC 2021). Ongoing 


Although black carbon is not a GHG, it affects climate in a variety of ways. Black carbon is emitted as a 
combustion byproduct, and the concentration of black carbon can vary spatially, seasonally, and vertically in the 
atmosphere (AMAP 2015; Creamean, Maahn et al. 2018; Stohl, Klimont et al. 2013; Xu, Martin et al. 2017). 
Black carbon affects the climate by absorbing and scattering solar radiation (i.e., sunlight). It can also influence 
clouds by altering the size and number of water droplets and ice crystals in water and ice clouds. Black carbon in 
cloud droplets decreases cloud albedo, which heats and dissipates the clouds. This also alters the temperature 
structure within and around the cloud, changing cloud distribution. 
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_e EO 13990, Protecting Public Health and the Environment and Restoring Science to Tackle the Climate 
Crisis included directives to establish an Interagency Working Group on Social Cost of Carbon to develop 
social costs associated with GHGs for cost-benefit analyses and to rescind the 2019 draft guidance from 
the Council on Environmental Quality (CEQ) entitled “Draft National Environmental Policy Act 
Guidance on Consideration of Greenhouse Gas Emissions” (84 FR 30097). 


InJ anuary 2021, two executive orders (EOs) were issued to address the climate crisis: 
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| Additional discussion of laws and policies relevant to GHGs and climate change is available in the BLM 
_ Specialist Report on Annual Greenhouse Gas Emissions and Climate Trends (2020) (herein referred to as the 
_ BLM Specialist Report). 


1.3. Observed Climate Trends 
1.3.1. Arctic* 
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| at least 1900 (Thoman, Richter-Menge et al. 2020). 


_ Spring snow cover extent, observed by satellites, has been decreasing over arctic land since 2005, especially in 

_ May and June (Derksen, Brown et al. 2017). The North American Arctic snow cover extent in June has been 
_ below the long-term average every year since 2006, and the complete 2020 snow-free period in the Arctic was the 
_ second-longest since recording started in 1998 (Moon, Druckenmiller et al. 2021). With decreased snow cover 
_ extent and shorter snow cover duration in the Arctic, more of the sun’s energy is absorbed by the dark land 

_ surface, warming the surface further. This results in a reinforcing feedback effect that further reduces snow cover 

_ (Melillo, Richmond et al. 2014). 


_ The extent of sea ice in the Arctic is also decreasing. Since the early 1980s, average annual sea ice extent has 
decreased by 3.5% to 4.1% per decade and the annual minimum sea ice extent, which occurs in September, has 

_ decreased at a rate of 11% to 16% per decade (USGCRP 2018). The 15 lowest September sea ice extents in the 

| satellite record (since 1979) have all occurred in the last 15 years (Moon, Druckenmiller et al. 2021). The extent 

_ of very old ice (4 years or older), which is thicker and more resilient to short-term temperature changes, has also 
_ decreased, with old ice only comprising 4.4% of ice cover in the Arctic Ocean in March 2020 compared to 33% in | 
_ 1985 (Perovich, Meier et al. 2020). Similar to decreases in snow cover extent, decreased sea ice extent alsohasa 
_ feedback effect on climate. An increased amount of the sun’s energy is absorbed by the open ocean relative to 

_ oceans covered by ice, leading to an increased rate of sea ice melting. Reductions in sea ice also make the Arctic 

_ more accessible by ships for transportation, oil and gas exploration, and tourism. This can lead to increased GHG 
emissions as well as other risks such as oil spills and drilling or maritime-related accidents (Melillo, Richmond et 

_ al. 2014). Rising air temperatures over land affect the Arctic permafrost layer. Permafrost is material that exists at 
_ or below 32°F (0°C) for at least 2 years, and the active layer is the layer above the permafrost that thaws 
_ seasonally. The northern circumpolar permafrost zone stores 1,700 petagrams (or 1,700 gigatons) of organic 
_ carbon, locked in place due to the slow rate of plant material decomposition in the frozen ground (Schuur, Abbott 
_ et al. 2013). With rising temperatures and decreasing snow cover, the permafrost extent is predicted to decrease 
_ significantly by the year 2100 (Slater and Lawrence 2013). Thawing permafrost releases CO2 and CH, to the 
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| raised through the public comment process.” 


Appendix E.2A Climate and Climate Change Page 2 


Willow Master Development Plan Final Supplemental Environmental Impact Statement 


_ atmosphere and delivers organic-rich soils to the bottoms of lakes, resulting in decomposition that releases 
_ additional CH4. Recent studies (Voigt, Marushchak et al. 2017) suggest that thawing permafrost could also lead to _ 
_ the release of significant amounts of N2O. These emissions can accelerate climate feedback effects (Jones, Irrgang 
_ et al. 2020; Markon, Trainor et al. 2012). 


LA reduction in sea ice has led to increased primary productivity (i.e., the rate at which energy is converted through 

_ photosynthetic and chemosynthetic processes into organic substances) in the Arctic Ocean (Moon, Druckenmiller | 
et al. 2021). Warmer temperatures combined with reduced ice cover have led to the greening of the tundra and 
_ increases in soil moisture and the amount of snow meltwater available. These changes have led to an increased 

_ active layer depth, changes in herbivore activity patterns, and reductions in human usage of the land due to 

| ground being frozen for a shorter period of time (Clement, Bengtson et al. 2013; Epstein, Bhatt et al. 2017). 

_ Although the greening of the tundra can store carbon as biomass, the effect of these changes in the Arctic has 

_ been a net release of carbon into the atmosphere (Epstein, Bhatt et al. 2017; Richter-Menge, Overland et al. 2017). | 


| Black carbon has a magnified impact on climate in the Arctic due to snow and ice albedo feedback. This feedback — 
_ occurs when black carbon settles on top of snow or ice and decreases the reflectivity (albedo) of the surface. This | 
_ allows more heat to be absorbed by the surface, leading to increased melting, which further decreases the albedo. 
| This feedback is prominent in the Arctic because so much of the surface is snow and ice, which have high albedo. 
The IPCC (2021) reports | that there 1 is “high confidence” that snowmelt in the Arctic is enhanced by deposition of 


1.3.2 North Slope 

Similar to the Arctic as a whole, the North Slope has experienced increased average temperatures, decreased sea 
ice and snow cover extent, an expanded growing season, and thawing permafrost. Temperatures in the North 
Slope have been warming at a rate 2.6 times faster than the continental U.S. (USGCRP 2018). Permafrost loss in 
Alaska’s North Slope is already widespread and progressive deep thawing of permafrost in the North Slope region 
may begin in 30 to 40 years (Thoman, Richter-Menge et al. 2020).Over the 35-year record (1982 to 2016), the 
North Slope has shown substantial increases in tundra greenness (Richter- Menge, Overland et al. 2017). A 
warming climate, in addition to regulatory changes and methods for measuring frost depth, has contributed to a 
reduction in the tundra travel season from 200 days in the 1970s to less than 120 days in 2003 (NSB 2014). With 
continued climate warming and precipitation changes, the tundra travel season is expected to shorten further. 
Since the mid-1980s, Alaskan permafrost on the Arctic coast has warmed between 6°F to 8°F at a depth of 3.3 feet 
(1 meter [m]). In 2016, all but one permafrost observational site documented record high temperatures ata depth 
of 65.6 feet (20 m) on the North Slope. Depth temperatures at 65.6 feet (20 m) in this region have been increasing 
at rates between 0.38°F and 1.19°F per decade since 2000. The active layer depth was at a 210-year maximum on 
the North Slope in 2016 (Richter-Menge, Overland et al. 2017). 


1.4 Observed Greenhouse Gas Trends 
ac 4.1 National” 


bana Sinks. In 2019, 6, 558. 3 million metric tons (MMT) of COze were emitted in the U.S. (EPA 2021d). The major | 
_ economic sectors contributing to GHG emissions in the U.S. in 2019 were transportation (28.6%), electricity 

_ generation (25.1%), industry (22.9%), and agriculture (10.2%). CO» from fossil fuel combustion has accounted 

| for approximately 76% of U.S. GHG emissions since 1990, and the U.S. accounted for approximately 15% of 

_ global CO> emissions from fossil fuel combustion in 2018 (EPA 2021d). These fossil fuel combustion CO2 

_ emissions increased by approximately 2.6% between 1990 and 2019 but decreased by approximately 15.6% from 

| 2005 levels. 


|The 2019 U.S. emissions inventory was used as it was available at the time of initiating this analysis. While there 
_ was a large change (-10.6%) in U.S. emissions from 2019 to 2020, this was large due to the impacts of the 

_ coronavirus (COVID-19) pandemic on travel and economic activity (EPA 2022a) and therefore 2020 is not 
_ considered a representative year for analysis. Thus, the results are presented in this Supplemental EIS with respect 
to the 2019 USS. total emissions. 
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1.4.2 Alaska* 

The EPA documents GHG emissions from Alaska in the Alaska Greenhouse Gas Emissions Inventory. Emissions 
are calculated using a top-down approach, where emissions factors are applied to statewide activity data from 
1990 to 2015. In 2015, approximately 40 MMT COvze were emitted in Alaska. This is a decrease of approximately 
8% from 1990 levels and a decrease of approximately 23% from the peak emissions observed in 2005 (ADEC 
2018). 
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_ Nations Environment Programme (UNEP) Emissions Gap Report (2021) estimates global GHG emissions in 2030 
_ would need to be 55% lower than projected 2030 emissions to limit global warming to 1.5°C and 30% lower to | 
_ limit warming to 2°C. UNEP (2021) estimated that current pledges for 2030 reduce the projected 2030 emissions 

_ by only 7.5%. An analysis by Tong, Zhang et al. (2019) indicates that future global CO2 emissions anticipated 

_ from existing and proposed energy infrastructure already exceed the carbon emissions budget needed to limit 

_ global warming to 1.5°C; however, other studies suggest that attaining a 1.5°C warming limit is possible by 

_ replacing existing infrastructure with zero-carbon alternatives at the end of their life spans, enabling us to meet 

_ climate goals (Smith, Forster et al. 2019). For U.S. emissions, the U.S. Energy Information Administration 

_ (USEIA) estimates trends in future U.S. CO2 emissions in the Annual Energy Outlook 2021 Report (2021a). U.S. 

- CO» emissions are predicted to decrease from 2023 to 2035 as a result of a transition away from coal and a rise in 

_ natural gas and renewable energy, but emissions are then projected to trend upward after 2035 due to increasing 

_ population and economic growth, with the rate of increase depending on economic conditions. 


_ Climate projections under both higher (representative concentration pathway [RCP] 8.5) and lower (RCP 4.5) 

_ GHG emission scenarios shows that the state of Alaska should expect warmer annual temperatures, reduced snow 
_ cover and sea ice extents, thinner sea ice, and potential increases in the area burned by wildfire (USGCRP 2018). | 
_ Under RCP 8.5, the interior and northern areas of the state are projected to warm by 10°F to 16°F (BLM 2020). In 
_ coastal areas of the North Slope, the number of nights below freezing is projected to decrease by more than 45 
_ nights per year (BLM 2020). 


_ Climate projections for Alaska indicate that snow cover duration is expected to drop, with a later date of first 

_ snowfall and an earlier snowmelt, and the arctic waters could be virtually ice-free by late summer before 2050 

_ (BLM 2020; Markon, Trainor et al. 2012; Mudryk, Elias Chereque et al. 2020). Models predict permafrost 

_ thawing will continue, with some models predicting that near-surface permafrost will likely disappear on 16% to 

_ 24% of the landscape of Alaska by the end of the 21st century (BLM 2020; USGCRP 2018). This will impact 

_ rural Alaskan communities by likely disrupting sewage systems and community water supplies. The increasing 

_ trend in the length of the Alaska growing season is also projected to continue. This change will reduce water 

_ storage as well as increase the risk and extent of wildfires and insect outbreaks in the region. Warmer 
_ temperatures, wetland drying, and increased summer thunderstorms will likely continue to increase the number of — 
_ wildfires in Alaska (USGCRP 2018). ; 


Warmer temperatures in the Willow Master Development Plan (MDP) Project (Project) area will lead to a deeper 
active layer, which would affect the surrounding ecosystem. A deeper active layer would allow improved water 
drainage and the migration of deeper-rooted plant communities farther north. Changes in plant communities 
would also be driven by the expanded growing season and warmer, drier soils. These vegetation changes would 
promote soil formation as root development and organic matter in the soil profile increase. 


As the active layer deepens, damage from traffic over the surface during non-frozen periods would likely increase 
due to accelerated erosion and subsidence of permafrost. Permafrost thawing could also lead to thermokarst or 
slumping, resulting in increased nutrient loading and suspended sediment in lakes and rivers. Warmer 
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temperatures may lead to an increase in the frequency of lake tapping (sudden drainage) events as degrading ice 
wedges integrate into drainage channels at lower elevations. 


Arctic fish species will be affected by increased water temperatures as air temperatures increase, but this impact is 
difficult to predict. Arctic bird species will be affected by habitat loss as aquatic and semiaquatic habitats are 
converted into drier habitats. A reduction in available habitat would likely cause changes in bird distributions, 
increased competition for resources, and declines in productivity. 


Paleontological resources could be adversely affected by climate change, but the impact is difficult to determine. 
Paleontological sites may more rapidly decompose in a warmer climate, and sites on hillsides, bluff faces, 
riverbanks, and terraces may be destroyed by mass wasting. Erosion may lead to increased exposure of known 
paleontological sites. Many known paleontological sites in the Project area have been exposed due to erosion with 
few negative impacts. 


As with paleontological resources, cultural resources on the North Slope could also be impacted by mass wasting, 
warmer temperatures, and erosion. In addition, as the permafrost thaws and the active layer deepens, cultural 
resources may be incorporated into the active layer. These sites would then be exposed to cryoturbation (frost 
mixing) and subject to vertical disturbances that may cause sites at different vertical layers to become mixed. 
These disturbances can occur in both vertical directions as seasonal frost cracking can cause downward 
movement, and frost heaving and sorting, ice wedging, and involutions can push artifacts upward. 


Climate change may impact the accessibility of mineral material deposits on the North Slope. While the existence 
and location of these deposits will not be affected, the excavation process may be made easier, due to the thawing 
permafrost, or more difficult, as developing deposits in areas with thawed permafrost may require water removal 
or excavation in swampy conditions. 


2.0 ANALYSIS METHODS 
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_ facilitate multicomponent climate policies by allowing emissions of different GHGs and other climate forcing 

_ agents to be expressed in a common unit (CO:e). The global warming potential (GWP) was introduced in the 

_ IPCC’s first assessment report, where it was also used to illustrate the difficulties in comparing components with 
_ differing physical properties using a single metric. Each GHG has a GWP that accounts for the intensity of the 

_ GHG’s heat trapping effect and its longevity in the atmosphere. GHG emissions are reported in units of COze to 
_ account for the varying GWP of pollutants and to allow for more direct comparisons of the global warming 

| impacts of different GHGs. 


The 100-year GWP was adopted by the United Nations Framework Convention on Climate Change (IPCC 2014) 
and its Kyoto Protocol and is now used widely as the default metric. In addition, the EPA uses the 100-year time 
horizon in its Inventory of U.S. Greenhouse Gas Emissions and Sinks: 1990-2019 (EPA 2021d). 


The 100-year GWP is only one of several possible emission metrics and time horizons. The choice of emission 
metric and time horizon depends on the type of application and policy context; hence, no single metric is optimal 
for all policy goals. All metrics have shortcomings and choices contain value judgments, such as the climate 
effect considered and the weighting of effects over time (which explicitly or implicitly discounts impacts over 
time) and the climate policy goal and the degree to which metrics incorporate economic or only physical 
considerations. There are significant uncertainties related to metrics, and the magnitudes of the uncertainties differ 
across the metric type and time horizon. Three such metrics type/time horizon combinations are listed in Table 
E.2.1 and were used in the GHG analysis. In general, the uncertainty increases for metrics along the cause and 
effect chain from emission to effects. 


All Project GHG emissions were converted to units of CO2e for ease of comparison using the GWP values shown 
in Table E.2.1. 
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Table E.2.1. Global Warming Potential Factors* 
Rationale for Time Horizon 
Used by the U.S. Environmental Protection Agency in its GHG inventories and GHG reporting rule 
requirements under 40 CFR 98(a) (EPA 2019). 
Used by the IPCC in its climate change synthesis report of the sixth assessment report (IPCC 2021). 
The IPCC used different CHs GWPs for sources originating from fossil carbon and non-fossil 
carbon. The CHs fossil values were used here as all Project emissions originate from fossil carbon. 
Used by the IPCC in its climate change synthesis report of the sixth assessment report (IPCC 2021). 


The IPCC used different CHs GWPs for sources originating from fossil carbon and non-fossil 
carbon. The CHs fossil values were used here as all Project emissions originate from fossil carbon. 


bocce cedececcnncusdwacccceteccecceccnccsnccencssnesacenacesccnssocanccssccasdamecaeuacsoos 


2.2 Direct Greenhouse Gas Emissions Calculation Methods* 

ConocoPhillips Alaska, Inc. (CPAI) developed a Project emissions inventory (CPAI 2019) of all known emissions 
sources (e.g., vehicles, aircraft, drill rigs, generators) that would be present during the construction and life of the 
Project for Alternative B (Proponent’s Project). BLM reviewed the emissions inventories, and the Alternative B 
inventory was used as the basis for estimating emissions from Alternatives C (Disconnected Infield Roads) and D 
(Disconnected Access). In support of the Supplemental Environmental Impact Statement (EIS), CPAI developed a 
Project emissions inventory for Alternative E (SLR 2022c, presented in Appendix E.3B, Attachment B) which 
was also reviewed by BLM. GHG emissions were calculated for each alternative as part of these inventories to 
estimate the Project’s direct GHG emissions. 
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_ and module delivery. Emissions from these activities would come from stationary combustion sources, mobile on- | 
_ road and nonroad tailpipe combustion sources, fugitive sources, aircraft sources, and marine vessel sources. GHG 
_ emissions quantified from these activities include CO2, CH, and N20. Details on these activities and emissions 

_ are provided in Appendix E.3B, Air Quality T echnical Support Document. An additional discussion is provided 

_ below for methane which is a GHG of particular concern in oil and gas production due to its high GWP. 


_ Methane emissions due to oil and gas production activities are due to several point and non-point (fugitive) 

_ sources during the life of the facility, from construction to production operations. Methane gas, in addition to 

_ other hydrocarbon gases, is naturally entrained within the geologic formations and is extracted, often 

_ accompanying the oil, during the flowback and production phases of the project. Methane emission points related 
_ to production facilities include flares, combustors, blowdowns, tanks, pig launching, and component leaks (e.g., 

_ valves and pumps). Additional sources of methane are related to natural gas-fired engines, heaters, and turbines. 

_ Each of these sources has been identified and evaluated in this Supplemental EIS. 


_ The purpose of flares/combustors is to provide a required minimum amount of methane, ethane, and VOC 

_ destruction and removal efficiency (at least 95%) to prevent large releases (“venting”) into the atmosphere. Vapor 
| recovery units recover gases prior to release and recycle them back into the process. The Project would use 

_ instrument air (no methane) and electrically-driven process control devices rather than produced gas which 

_ minimize fugitive emission leaks. 


_ The amount of methane extracted from a wellsite is finite and it will either be captured for use, destroyed (95% 
destruction removal efficiency), or inadvertently leaked (components). Detailed methane emissions inventories 
_ for the Project were developed by identifying the potential emission points discussed above and quantifying the 
_ number of components (fugitive emissions), then using EPA AP-42 factors, representative gas analyses, and 

_ maximum volumetric flow data (inherent operations) to determine individual estimations for each point. The 
summation of these points is evaluated, and a site-wide methane emissions estimation is determined in each action | 
_ alternative. The sources of methane emissions assessed in the emission inventories are inclusive of well work- 
_ related venting, fugitive, and combustion-related emissions. 


| The Project will also be subject to the proposed rule, 40 CFR Part 60, Subpart OOOOb which aims to reduce 
| methane emissions by implementing expanded performance standards, methane limits, and rigorous monitoring 
_ (via optical gas imaging). The Project will not use methane venting to enhance production. 


| The GWPs shown in Table E.2.1 were used to calculate total COve based on emissions of CO2, CH4, and N20. 
Under Alternatives B, C, and E, the Project would have a 30-year life, while under Alternative D, the Project 
_ would have a 31-year life. For additional information regarding the methods used to estimate direct on-site 
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emissions for each alternative, see Chapter 2 and related attachments in the Willow MDP Supplemental EIS Air 
Quality Technical Support Document provided as Appendix E.3B, Air Quality Technical Support Document. 


Direct emissions of GHGs from business commuting of employees and contractors by air travel were also 

_ estimated. It is anticipated that most employees and contractors would commute from Anchorage to the North 

_ Slope (i.e., Willow or Alpine airstrips) using ConocoPhillips Global Aviation services or equivalent (SLR 2022b) 
_. CPAI projections of the number of flights per year to Willow and Alpine airstrips under each action alternative 

_ (CPAI 2019; SLR 2022b) were used to calculate direct GHG emissions from business commuting of employees 

_ using passenger flights. 


_ The International Civil Aviation Organization (ICAO) has developed a standard for calculating the amount of 

_ CO» emissions generated by a passenger on an aircraft (ICAO 2018). ICAO utilizes a distance-based approach 

_ with up-to-date aircraft-type data. The inputs needed to calculate the emissions are as follows: emission ratio per 

_ metric ton of jet fuel used, total fuel used, passenger to freight weight ratio, number of commercial seats available | 
_ on the flight, and passenger load factor. The total fuel used is calculated based on the ICAO fuel utilization tables, | 
_ which give the average fuel used for a given aircraft type per average trip distance (ICAO 201 8). The distances 
_ from Anchorage to the Willow and Alpine airstrips were calculated using the Google Maps straight distance 

_ calculator, and then an ICAO correction factor was applied to account for the great-circle distance between 

_ airports. ICAO provides CO: emission factors for jet fuel but not for N20, and CHy. Thus, CO2, CH4 and N20 

_ emission factors for conventional jet fuel were obtained from EPA (2021b). The energy density used was the 

_ higher heating value of jet fuel from the Greenhouse Gases, Regulated Emissions, and Energy use in 

_ Technologies Model database (Argonne National Laboratory 2021). All other characteristic values, including 


a a a a Sr NN ae gy a 8S a aT neha a a ae a 


| (GLEEM; Wolvovsky 2021) is used (with updates, as discussed below) to estimate indirect GHG emissions from 

_ domestic transportation, refinement, and oil usage. This model was developed to support the Outer Continental 

| Shelf Oil and Gas Leasing Program 2017-2022, and it represents the best available resource for estimating 

_ indirect GHG emissions from petroleum products refined and consumed domestically. A description of the 
_model’s capabilities and methodology can be found in Wolvovsky (2021). Updates were made to the model inputs | 
_ for the Willow MDP Supplemental EIS to incorporate additional data, as discussed below. For this Supplemental | 
_ EIS, GLEEM was used to estimate the downstream GHG emissions associated with consumption of the oil and 

_ gas produced from the Project as well as the energy substitutes (ranging from other oil sources to renewable 

_ sources). BLM’s EnergySub Model estimates these energy substitutes that could replace production from the 

| Project or, equivalently, be displaced due to the Project (see Appendix E.2B)’. Substitution rates from EnergySub 

_ were rounded to the nearest whole percentage for use in GLEEM. BOEM’s Office of Environmental Programs 

_ developed GLEEM to estimate the full lifecycle emissions from both production and consumption of Outer 

_ Continental Shelf resources. For this Project, only the downstream portion of the model is used, as the upstream 

_ component is derived in combination with an offshore-specific separate model. The use of BOEM’s GLEEM in 

_ the GHG analysis for the Project is limited to the emissions associated with the processing and consumption of oil | 
_ and gas resources and not the upstream emissions from actual production of the resources that were calculated 
separately, as discussed in Section 2.2, Direct Greenhouse Gas Emissions Calculation Methods. 


| The 2021 version of GLEEM was downloaded from BOEM’s website? on February 22, 2022, for use in this 
_ analysis. The following updates were made to GLEEM input data for use in this Supplemental EIS: 


_e All national mineral activity data (e.g., U.S. crude oil refinery inputs and refinery processing gain, U.S. 

petroleum product consumption) in GLEEM were updated to use the latest data available from the 

: USEIA’s Monthly Energy Review Report (USEIA 2022). 

e The national emissions used in GLEEM 2021 for crude oil refining are from the Petroleum Systems 

source category of EPA (2021d), which excludes all combustion emissions of CO2 except for flaring. 
EPA (2021d) includes the CO2 combustion emissions from crude oil refining in the industrial sector 
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emissions of the Fossil Fuel Combustion category. GLEEM was updated to use the total U.S. refinery 
GHG emissions reported under the GHGRP (Subpart Y of 40 CFR 98) that include stationary fuel 
combustion emissions as well as process emissions (e.g., flares, process units, vents, blowdowns, fugitive — 
leaks) (EPA 202 1a). | 
_e The national mineral activity data used in GLEEM is for the year 2020 from USEIA (2021b), while the 
national emissions data used in GLEEM is for the year 2019 from the EPA’s annual GHG emissions 
inventory (2021d). All national mineral activity and emissions data used as GLEEM inputs were updated 
to use a 5-year average of recent years (2015 to 2019) instead of a single year. 
_e GLEEM was updated to assume that all Project oil produced under the action alternatives (and energy 
| substitutes under the No Action Alternative) are combusted. This results in a conservatively high estimate — 
of combustion emissions as some oil and natural gas are used for non-combustible products (e.g., 
fertilizers). 
_e Minor corrections were made to the EPA (2021b) stationary combustion emission factors used to estimate : 
the downstream combustion emissions of propylene, petroleum coke, and industrial coal. 
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| The Project would increase total U.S. crude oil production, which the results from EnergySub indicate would 

_ reduce prices for oil and other energy sources and result in changes in both domestic and foreign energy 

_ consumption. The changes in domestic and foreign oil consumption because of Project production are estimated 

_ using the BLM EnergySub model (Appendix E.2B). The increases in oil consumption domestically and abroad 

_ would result in GHG emissions. Emissions from the change in domestic consumption of crude oil and other 
_ energy sources (e.g., coal, natural gas) under the No Action Alternative are estimated using GLEEM with updates 
_ to model inputs, as described above. Emissions from the change in foreign oil consumption under Alternatives B, 
_C, D, and E are estimated by applying an EPA (2021a) stationary combustion emission factor to the change in 

_ foreign oil consumption estimated by EnergySub. Due to the lack of information on the type and amount of 
_ petroleum products consumed in foreign markets, the highest emission factor (11.91 kilograms of CO: per gallon, 
_ 0.47 gram of CHs per gallon, and 0.09 gram of N2O per gallon) reported by EPA across all petroleum products 
_ (EPA 2021b) are used and all oil are assumed to be combusted for a conservatively high estimate of emissions. 


: In addition to the indirect emissions estimated by GLEEM, indirect GHG emissions from the transport of Willow 
_ oil via pipeline and barge and deliveries of diesel fuel to the Project via barge, rail, and truck were also estimated 
as described below. 05 ee 
2.3.1. Transport of Project Oil to Refineries via the Trans-Alaska Pipeline System and Polar 
Tankers* 
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_ Sales-quality crude oil processed at the Willow Processing Facility (WPF) would be transported through the 
_ Willow Pipeline to a tie-in with the Alpine Sales Pipeline near drill site Colville Delta 4 North. The oil would then 
_ travel through the Alpine Pipeline to the Kuparuk Pipeline and then to the Trans-Alaska Pipeline System (TAPS) 
_ near Deadhorse, Alaska. From there, the oil would travel through TAPS to the Valdez Marine Terminal located in | 
_ southern Alaska, where it would be loaded onto polar tankers to be transported to refineries. To estimate 
additional indirect GHG emissions from the transport of Project oil via TAPS, emissions of CO2, CHs, and N20 

_ from the four active TAPS pump stations (i.e, TAPS Pump Stations 3, 4, 7 and Alyeska Pipeline Pump Station 

_ 01) and the Valdez Marine Terminal were obtained from the EPA’s Facility Level Information on Green House 

_ gases Tool for the period 2015 to 2019 (EPA 2022b). The annual reported GHG emissions from TAPS and the 

_ Valdez Marine terminal were then divided by the total annual TAPS throughput (Alyeska Pipeline Service 

_ Company 2022) to estimate the emissions intensity (i-e., metric tons of CO2, N2O, and CH, per barrel of oil) from 

_ transport. The average emissions intensity from 2015 to 2019 was multiplied by the yearly Willow production 

_ under each action alternative to obtain an estimate of the annual GHG emissions from the transport of Project oil 

_ through TAPS. A similar methodology was used to estimate annual GHG emissions from the transport of Willow 

_ oil on polar tankers from the Valdez Marine Terminal to refineries. Emissions intensities for the polar tankers 
| (e.g., metric tons of CO2, CHa, and N20 per millions of barrels of oil transported) were obtained from CPAI (SLR 
_ 2022b) and multiplied by the annual Project production under each action alternative to obtain annual emissions 
: estimates. 
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|e Barge from Valdez to the Port of Anchorage 
_¢ Rail from the Port of Anchorage to Fairbanks 
/e Truck from Fairbanks to Deadhorse, then to Kuparuk 


_ To estimate the emissions from barge and rail transport of liquid fuels, an EPA (2020) guidance approach was 

_ used that estimates emissions from gross ton-miles. Alaska Railroad (201 5) reported that “it takes just one gallon 
_ of fuel to move a ton of freight the length of the entire Railbelt.” The system map for Alaskan Railroad indicated 
_ that the entire Railbelt from Seward to Fairbanks is 470 miles long (Alaska Railroad 2020), which translates to 

_ 470 ton-miles per gallon of diesel fuel consumed. 


_ TTI (2017) reported a similar freight fuel efficiency for rail and provided barge transport efficiency as well. The 
_ TTI values were used to estimate the emissions for barge (inland towing) and rail freight moves. Note that the 

_ ton-miles are for freight moves and returning empty (deadhead) is incorporated in the overall efficiency 

_ represented here: 


e 647 ton-miles per gallon for inland towing | 
: e 477 ton-miles per gallon for railroads 
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_ EPA (2009) provides locomotive engine emission factors directly in gram per gallon units for different railroad 

_ authorities to account for the expected fleet age distribution and other factors. The Small Railroads category was 
_ conservatively used for the emission factors as a high emissions case because many railroads run older engines 
_and on higher emitting switching duty. The same diesel fuel carbon density was used as for towboats, above (i.e., 
_ 10,217 grams CO: per gallon), and the EPA (2021b) mobile combustion emission factors for diesel locomotives 

_ of 0.8 and 0.26 gram per gallon were used for CH, and N;O, respectively. Multiplying the freight tonnage by the 
_ distance moved (one way) provides ton-miles and dividing by the freight transport efficiency estimates the fuel 

_ consumed by mode. Then the fuel consumption multiplied by the emission factors in gallon units provides the 

_ expected emissions from freight transport. 


_ Estimates from emissions of liquid fuel transport via truck were calculated using the latest version of the EPA’s 

_ Motor Vehicle Emission Simulator (MOVES), MOVES3 (EPA 2021e). MOVES was run in inventory rate mode 

_ for the state of Alaska. Based on information provided by CPAI (SLR 2022b), the vehicle type chosen was diesel 

_ combination long-haul truck. The model was run for the first project year (2023), the fifth year (2027), and the 

_ tenth year (2032), from which emissions levels were assumed to be constant during the remaining Project life. 

_ This is a conservative assumption as equipment turnover over time would likely decrease emissions. The model 

_ emissions and activity were output on an annual basis by vehicle type, fuel type, road type, and calendar year and 

_ aggregated annually across all model years representing the MOVES default national age distribution for diesel 

_ combination long-haul trucks. Year 1 (2023) was used as a conservative surrogate for Years 1 to 4, Year 5 (2027) 

_ was used as a conservative surrogate for Years 5 to 9, and Year 10 (2032) was used as a surrogate for emissions in : 
_ Year 10 through the end of the Project. Running, short-term idling, and extended idling emission factors were 
_ then calculated using output emissions and their respective activity surrogate; extended idle hours for long-term 

_ idling (mandated driver rest), source hours operating for short-term idling (idling of 1 hour or less is expected to 

_ happen during travel milestone stops), and mileage for running exhaust for all GHGs. Emissions were then 
calculated using annual activity under each alternative provided by CPAI (SLR 2022a) alongside calculated 
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3.0 ENVIRONMENTAL CONSEQUENCES 
3.1 Effects of the Project on Climate Change 


3.1.1 Alternative A: No Action* 
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: Under Alternative A, the Project would not be developed and direct and indirect GHG emissions from the Project 
: would not occur and hence not contribute to climate change. Current trends in global, U.S., and Alaska GHG 

_ emissions would continue, unaffected by the Project. Energy demand would potentially be satisfied by non- 

_ Project sources, ranging from other oil sources to renewable sources. The BLM EnergySub Report (Appendix 

_ E.2B, Bureau of Land Management Energy Substitution Model (EnergySub)) presents an estimate of the amount 
_ of Project crude oil production that would be substituted by replacement (“substitute”) energy sources in the No 

_ Action Alternative (Alternative A). As described in Section 2.3, the substitution rates estimated by BLM 

_ EnergySub are used as inputs to BOEM’s GLEEM (Wolvovsky 2021) with updates to estimate emissions from 

_ the substitute energy sources for the Project. 


_ The GHG emissions estimated for the substitute energy sources under the No Action Alternative include the 

_ transport, processing, and downstream combustion emissions as calculated by GLEEM; on-site emissions for the 
_ substitute energy sources were not calculated and thus are not included in the emission presented here. The 

_ substitution rates are a function of the Project production, and thus the substitution rates estimated by BLM 

_ EnergySub vary slightly across the action alternatives due to differences the amount and timing of production 

_ under each alternative. The GHG emissions under Alternative A (No Action) are shown in Tables E.2.2, E.2.3, 

_ and E.2.4 based on the substitution rates under Alternatives B and C, Alternative D, and Alternative E, 

| respectively. 


Table E.2.2. Greenhouse Gas Emissions (metric tons) from Substitute Energy Sources under Alternative A 
(No Action) based on Substitution Rates for Alternative B and Alternative C* 


CO: CH, N2O COze COre 


EVeat 55.0 eG ee Ona 
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20,394,922.0 | 1,561.0 20,496,859 20,579,124 
22,451,640.0 | 1,718.0 224.0 22,563,988 22,654,527 


17,760,706.0 | 1,359.0 177.0 17,849,525 17,921,145 
15,670,219.0 | 1,199.0 156.0 15,748,537 15,811,725 


20,026,812.0 | 1,532.0 200.0 20,127,066 20,207,802 
13,876,925.0 | 1,062.0 138.0 13,946,247 14,002,214 


12,468,631.0 954.0 124.0 12,530,912 12,581,188 
11,242,708.0 860.0 112.0 11,298,912 11,344,234 
10,080,949.0 771.0 100.0 10,131,225 10,171,857 


7,260,986.0 556.0 7,297,211 7,326,512 
6,440,324.0 493.0 6,472,487 6,498,469 


5,481,198.0 419.0 5,508,699 5,530,781 
4,792,250.0 367.0 4,816,291 4,835,632 
4,160,713.0 318.0 4,181,382 4,198,141 


3,269,132.0 250.0 3,285,591 3,298,766 
2,958,429.0 226.0 2,973,081 2,984,991 
2,725,402.0 209.0 2,739,001 2,750,016 


9,027,261.0 691.0 90.0 9,072,423 9,108,839 
8,179,584.0 626.0 8,220,625 8,253,615 
Year 26 2,519,393.0 193.0 2,531,969 2,542,141 


3,694,660.0 283.0 3,713,194 3,728,109 


2,134,392.0 163.0 2,144,982 2,153,573 
2,090,488.0 160.0 2,100,989 2,109,421 
1,911,497.0 146.0 1,921,035 1,928,729 


Appendix E.2A Climate and Climate Change Page 10 


, a ee ee ee, ee eee eee 


Willow Master Development Plan Final Supplemental Environmental Impact Statement 


AR6 100-year GWPs AR6 20-year GWPs 
_Pvears0e 1} 1,762,900.0, || | 135.0 1,778,952 
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Table E.2.3. Greenhouse Gas Emissions (metric tons) from Substitute Energy Sources under Alternative A 
No Action) based on Substitution Rates for Alternative D* 


CO> CH, N20 COre COre 

AR6 100-year GWPs (AR6 20-year GWPs 
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22,451,640.0 
20,026,812.0 
17,760,706.0 
15,670,219.0 
13,876,925.0 
12,468,631.0 
11,242,708.0 
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4,816,291 4,835,632 
4,181,382 4,198,141 
3,713,194 3,728,109 
3,285,591 3,298,766 
2,973,081 2,984,991 
2,739,001 2,750,016 
2,531,969 2,542,141 
2,144,982 
2,100,989 
1,921,035 
LAL S837 


48.0 
41.0 
37.0 
Year 24 33.0 
Year 25 29.0 
Year 26 27.0 
Year 27 25.0 
Year 28 21.0 
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_[Year 31 18.0 
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Table E.2.4 Greenhouse Gas Emissions (metric tons) from Substitute Energy Sources under Alternative A 
(No Action) based on Substitution Rates for Alternative E* 
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13,453,655 


12,070,026 


10,649,083 


4,707,820.0 


3.1.2 Alternative B: Proponent’s Project 
Alternative B direct and indirect COxe emissions are quantified and described in the following sections. Black 
carbon effects on climate are also discussed. 


3.1.2.1 Direct Greenhouse Gas Emissions* 

Direct and indirect emissions of the GHGs CO2, CHa, and N2O will impact the climate. The Project is also 
expected to produce a small amount of sulfur dioxide, a GHG that has an overall cooling effect; however, the 
effect of sulfur dioxide emissions would be negligible. Direct emissions for the Project include, but are not limited 
to, emissions from vehicle traffic, air traffic, power generation, and drill rigs. 


GHGs have long lifetimes (i.e., 10 to 100 years) before they are chemically broken down or otherwise removed 
from the atmosphere through absorption or deposition. Since GHGs are relatively stable, changes in GHG 
emissions have long-lasting effects on the climate. Alternative B direct GHG emissions estimated over the 30- 
year Project lifetime are provided in the main body of this Supplemental EIS (Section 3.2.2.3 of the SEIS). 
Emissions are given in COve units to account for the GWP of pollutants and were calculated using GWP values 
for both 100-year and 20-year time horizons (Table E.2.1). The annual average gross emissions under Alternative 
B shown are shown in Table E.2.2: these emissions do not account for the market substitution effects discussed in 
Section 3.2.2.2, Alternative A: No Action. The annual direct gross emissions under Alternative B are provided in 
Table E.2.6 for each year of the life of the Project. 


Table E.2.5. Annual Average Gross Greenhouse Gas Emissions (metric tons) under Alternative B 
(thousand metric tons per year)* 


; CQve : Re 
GHG COvre CQre 
Pee eee oi N20. ie (10-year AR6 GWP) | (20-year AR6 GWP) 


0.2945 | 0.0018 
8,651 0.614 0.089 8,693 8,694 8,726 

9,415 0.9087 0.0908 | 9,465 9,467 9,515 

Note: AR4 (fourth assessment report of the Intergovernmental Panel on Climate Change [IPCC]); AR6 (sixth assessment report of the IPCC); CH, (methane); 
CO; (carbon dioxide); COve (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential); N2O (nitrous oxide). Year 0 only included 
1 month of construction activity and thus this year was excluded from the average annual emissions. 

*Total values may have small differences due to rounding. 
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Table E.2.6. Ann 


ual Direct Greenhouse Gas Emissions (metric tons) under Alternative B* 
CO2 CHs N20 COzre COre 
AR6 100-year GWPs AR6 20-year GWPs 


Year 7 943,887 960,227 
Year 8 
Year 9 
Year 10 
Year 11 
Pes ON See so os, tm Le A 8098 et a 893,800) ee 


ar 18 863,961.2 356.6 875,098 893,890 
863,961.2 356.6 875,098 893,890 


863,961.2 356.6 875,098 893,890 
| 863,961.2 | 


er 10 | 
oe Aad 
eee) | 


863,961.2 356.6 875,098 893,890 


24 
Year 29 
Year 30 


re 
elgg |g |g 
BRIS 18 JF 
NININ IN 
CO; DN 


3.1.2.2 Indirect and Total Domestic Greenhouse Gas Emissions* 

Indirect emissions are expected to come from transportation, refinement, and downstream consumption of the oil 
extracted by the Project. Natural gas extracted from the Project would be either beneficially used onsite or 
reinjected into the well and would not be transported for consumption. 


sccananamenanceomenarsenpurenumsesensnsesenassussvasscsunennsnenansaue-secenaansucssnancpameumpanunpesneannnsennencacnnarssanosiissesanasinensinennan=nnS=nsnannewepeeerenseeehena tare aeaeenn teannnnn anne aTn nn SAAR DARA ERne en PAPO MAAS ONTO eee ae oh enna Pate nnenenTRseeOnOeeninnnnn ARPA MOOR ON Onna nA AAP OAD OTTO OANA ORONO ONO CESS TAO ROAND AOS OE OSSD 


_ provided in Table E.2.7 for each year of the life of the Project. 


| Table E.2.8, Table E.2.9, and Table E.2.10 show the total (gross and net) domestic GHG over the Project duration 
_ for each action alternative using three sets of GWPS: IPCC AR4 100-year, IPCC AR6 100-year, and IPCC AR6 
_ 20-year, respectively. When applying the 100-year GWPs from the IPCC AR4, Alternative B’s annual average 

_ direct GHG emissions (772 thousand metric tons [TMT] of COve per year) over the 30-year Project life are 

_ approximately 1.9% of the 2015 Alaska GHG inventory. The annual average total gross (i.e., sum of direct and 

_ gross indirect) GHG emissions of 9,465 TMT of COsze per year represents approximately 0.1% of the 2019 U.S. 

_ GHG inventory (note that the indirect emissions are compared to the national totals and not Alaska totals as most 
_ of the indirect use is expected to occur outside Alaska). These emissions would represent approximately 0.3% of 
_ the U.S. net GHG emissions target for 2030. When applying the 100-year GWP from the IPCC AR6, Alternative 
_ B’s annual average direct GHG emissions (774 TMT of COve per year) are approximately 1.9% of the 2015 
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| Alaska GHG inventory. The annual average total gross GHG emissions are 9,467 TMT of COze per year; they 

_ constitute approximately 0.1% of the U.S. GHG inventory and approximately 0.3% of the U.S. net GHG 
_ emissions target for 2030. When applying the 20-year GWPs from the IPCC AR6, Alternative B’s annual average © 
_ direct GHG emissions (789 TMT of COvze per year) are approximately 2.0% of the 2015 Alaska GHG inventory. | 
_ The annual average total gross GHG emissions of 9,515 TMT of COve per year represent approximately 0.1% of 

_ the 2019 U.S. GHG inventory and approximately 0.3% of the U.S. 2030 net GHG emissions target. In all three 


Table E.2.7. Annual Gross Domestic Indirect Greenhouse Gas Emissions (metric tons) under Alternative 


B* 

AR6 100-year GWPs AR6 20-year GWPs 
hes ne Mi ee ee 
PYear2 [1095.55]. Or ae 0 eee ee ee ee ee ie 
Vear3. fe 1868.8 1 020 ee 0.0 ae Es el ar enero! eeeee | oa ee 
Year 4 fe LAUR PO 00 ee ne 
rYearS | | 28394 0 | Oat e ie lOc a eee eee ct eee 

12,318,047.8 
11,030,644.0 
9,994,818.1 


8,872,462.1 
7,869,749.4 
6,697,841.1 
5,856,044.1 
5,084,388.3 
4,514,907.6 


Year 18 


rl < I< re I< I< I< 
Oo;o0|0 Ae ects ba 
Be Yee [ee (ee Be IB 18 |S 
NTN] — fe eee 
m1O)}\O SID] nn] & 


KK) < << 
8 |8 18 [8 |S 
BIR IS ia ia 
NINININI]NY 
DD} Ur] & |b |b 


3,995,029.4 
3,615,414.0 | 256.5 
3,330,640.2 
3,078,928.5 
2,608,536.7 
2,554,870.0 
2,336,172.2 
[rearg0° =o") 2,154,650.9.7|" 0153.4 | OO 180 ne ose | ee ee 


: Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CH4 (methane); CO2 (carbon dioxide); CO2e 
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Table E.2.8. Total (Gross and Net) Domestic Greenhouse Gas Emissions (thousand metric tons) over 
Project Duration for Each Action Alternative Based on 100-Year Time Horizon Global 
Warming Potential Values from the Intergovernmental Panel on Climate Change Fourth 
Assessment Report* 
Alternative GHG Gross COze COne Net COzxe Change 
Emissions Resulting from | in No Action from No Action 
Type Project? Alternative Alternative‘ 
(Substitute Energy 
Sources)? 


B: Proponent’s Project Direct 

B: Proponent’s Project Indirect 
|B: Proponents Project Total 

C: Disconnected Infield Roads Direct +25,326 

C: Disconnected Infield Roads Indirect +47,378 

C: Disconnected Infield Roads Total 

D: Disconnected Access Direct +23,276 

D: Disconnected Access Indirect 

D: Disconnected Access Total 
| E: Three-Pad Alternative (Fourth Pad Deferred) Direct +23,191 

E: Three-Pad Alternative (Fourth Pad Deferred) “indirect 

E: Three-Pad Alternative (Fourth Pad Deferred Total +69,395 


Note: COve (carbon dioxide equivalent); GHG (greenhouse gas) ; NA (not applicable) Project Reon would be 30 years under 
_ Alternatives B, C, and E, and 31 years under Alternative D. The global warming potential values used are carbon dioxide = 1; methane = 
W205: nitrous oxide = 298. 
| ® Indirect gross COve is from the Willow Project’s indirect GHG emissions estimated using the Bureau of Ocean and Energy 
| Management’ s (BOEM) Greenhouse Gas Life Cycle Energy Emissions Model (Wolvovsky 2021) with updates described in Appendix 
| E.2A. Numbers may not match exactly due to rounding. 
| © COve from Energy Sources Displaced by Project is estimated using the substitution rates modeled by BLM EnergySub (Appendix E.2B) 
: and in GLEEM with updates. Numbers may not match exactly due to rounding. Substitution rates from EnergySub were rounded to the 
| nearest whole percentage for use in GLEEM. 
© The net COze change is the difference between the previous columns. The + sign indicates an increase in emissions relative to Alternative 


Saaremaa ananassae recast ene eee da ena casa ae stad cater Nabe ene ene e aoa ee See a edac toe eee dee in dssdoanawtan d0 Ueeretar rere eeceateo des tebee ce tantenanewauserencueee sonctn weasedusnsesbustnbaseancunrsewesorarerseseetovbesesbonansbaauneaentasanaanasennnanenbnsrerneerseresnanennennseweeeTeSeN NaS aaaae nen ewewS 


nea 
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Table E.2.9. Total (Gross and Net) Domestic Greenhouse Gas Emissions (thousand metric tons) over 
Project Duration for Each Action Alternative Based on 100-Year Time Horizon Global 
Warming Potential Values from the Intergovernmental Panel on Climate Change Sixth 


Assessment Report* 
Alternative GHG Gross CO2e 
Emissions Resulting from 


COvze 
in No Action 


Net COze Change 
from No Action 


Type Project* Alternative Alternative “ 
(Substitute Energy 
Sources)? _ 


B: Proponent’s Project Direct 
B: Proponent’s Project Indirect 
B: Proponent’s Project Total 
C: Disconnected Infield Roads Direct 
C: Disconnected Infield Roads Indirect 
C: Disconnected Infield Roads Total 
D: Disconnected Access Direct 
D: Disconnected Access Indirect 
D: Disconnected Access Total 
|E: Three-Pad Alternative (Fourth Pad Deferred) Direct 
| E: Three-Pad Alternative (Fourth Pad Deferred) _ Indirect 
E: Three-Pad Alternative (Fourth Pad Deferred) 


pen SS 5 5555S SS SSBB SSB EN 


| 29.8; nitrous oxide = 273. 
| * Indirect gross COze is from the Willow Project’s indirect GHG emissions estimated using the Bureau of Ocean and Energy 

| Management’s (BOEM) Greenhouse Gas Life Cycle Energy Emissions Model (Wolvovsky 2021) with updates described in Appendix 

_ E.2A. Numbers may not match exactly due to rounding. 

: > COve from Energy Sources Displaced by Project is estimated using the substitution rates modeled by BLM EnergySub (Appendix E.2B) 

_ and in GLEEM with updates. Numbers may not match exactly due to rounding. Substitution rates from EnergySub were rounded to the 

| nearest whole percentage for use in GLEEM. 

| © The net COve change is the difference between the previous columns. The + sign indicates an increase in emissions relative to Alternative 
: A (No Action). 


fg ee Ee SSS 
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Table E.2.10. Total (Gross and Net) Domestic Greenhouse Gas Emissions (thousand metric tons) over 
Project Duration for Each Action Alternative Based on 20-Year Time Horizon Global 

Warming Potential Values from the Intergovernmental Panel on Climate Change Sixth 
Assessment Report* 


Alternative GHG Gross COze COvre Net COze Change 
Emissions Resulting from | in No Action from No Action 
Type Project* Alternative Alternative “ 
(Substitute Energy 
Sources)? 
B: Proponent’s Project 
B: Proponent’s Project Indirect 
B: Proponent’s Project Total +71,155 
C: Disconnected Infield Roads Direct 
C: Disconnected Infield Roads Indirect 
C: Disconnected Infield Roads Total 
D: Disconnected Access Direct 
D: Disconnected Access Indirect 
D: Disconnected Access Total 
E: Three-Pad Alternative (Fourth Pad Deferred) Direct 


E: Three-Pad Alternative (Fourth Pad Deferred) Indirect 255,359 209,046 +46,313 
E: Three-Pad Alternative (Fourth Pad Deferred) Total 


_LE: Three-Pad Alternative (Fourth Pad Deferred) =| «Coal 279,034 209,046 +69,988 é 


Note: CO2e (carbon dioxide equivalent); GHG (greenhouse gas); NA (not applicable). Project duration would be 30 years under 

| Alternatives B, C, and E, and 31 years under Alternative D. The global warming potential values used are carbon dioxide = 1; methane = 

: 82.5; nitrous oxide = 273. 

i * Indirect gross CO2e is from the Willow Project’s indirect GHG emissions estimated using the Bureau of Ocean and Energy 

| Management’s (BOEM) Greenhouse Gas Life Cycle Energy Emissions Model (Wolvovsky 2021) with updates described in Appendix 

: E.2A. Numbers may not match exactly due to rounding. 

i © COve from Energy Sources Displaced by Project is estimated using the substitution rates modeled by BLM EnergySub (Appendix E.2B) 

: and in GLEEM with updates. Numbers may not match exactly due to rounding. Substitution rates from EnergySub were rounded to the 

: nearest whole percentage for use in GLEEM. 

| © The net COve change is the difference between the previous columns. The + sign indicates an increase in emissions relative to Alternative 
i A (No Action). 


peceewennnnanennnenennnnarepenanaranecrannnnnseneace SEP -nacnenansnanncsnonseswnncanansnenansnnneennnennarnnranneeneneneeennnnnnnnnheeeemnennnnnnenenennnnen nnnnn naan nO PO PON NORCO ete O ne Oe nN NOT nena Renn e teen Rent eet nett nA Oe A An Aan nnnnt 


: The Project would increase total U.S. crude oil production which would reduce prices for oil and other energy 

_ sources and result in changes in both domestic and foreign energy consumption. The changes in domestic and 

_ foreign oil consumption as a result of Project production are estimated using the EnergySub model (Appendix 

_ E.2B). The increases in oil consumption domestically and abroad would result in GHG emissions. Emissions from 
_ the change in foreign oil consumption are estimated by applying an EPA stationary combustion emission factor to 
_ the change in foreign oil consumption estimated by the EnergySub Model. Due to the lack of information onthe | 
_ type and amount of petroleum products consumed in foreign markets, the highest emission factor (11.91 

_ kilograms of CO: per gallon, 0.47 grams of CH, per gallon, and 0.09 grams of N20 per gallon) reported by EPA 

_ across all petroleum products (EPA 2021b) was used and it was assumed that all foreign oil was combusted 

_ resulting in a conservatively high estimate of these downstream combustion emissions. The annual downstream 
_ emissions resulting from the change in foreign oil consumption over the life of the Project under Alternative B are 
| shown in Table E.2.11. | 
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Table E.2.11. Annual Downstream Greenhouse Gas Emissions (metric tons) Resulting from the Change in 
Perera on Consum baan under Aternatne B* 


COnzre CQre 
AR6 100-year GWPs AR6 20-year GWPs 
BR ee eG ae re ee 
47.0 G23 5515 6,248,239 


i VearsT in eet 6,851,698.0 D710 A tee ae 5d Somat 6,873,890 6,888,140 
Year 8-5) F/G 115.5780 DAV3 0 © [ens eee 6,135,386 6,148,105 
eS ol at 
pee 20 OR ee 


245.3 


Year9 | 5,297,644.) 209.10 [i 40s. 0] ese Ss) ge ue Jet 025 70a 
[Year 10 <7 “| 4,543,065-7 [S793 | 343 as eis ee ees 200) a 
Year 11 [40262716 |" 1158.95. 304 ee eo ee 
[Year 12 | 3,658,167. [144.40 SN 29.6) Oa Gr 0s8 a ees orn ct eee 
| Mean ee 3318 468 Oe ae a Oe fc siscrela 320,216 ae | eae ee 330,18 
Year 14. |. 2,944,280.0 [5116.25 22.29 1 E098 ohG gee gee Oe 
[Year 15°) | .2,648,164.07 [0 1045. |e 2) 20:00 a RIS Sores gee see1002 240 
[Year16 | | 2,184833/1 7 | 86.2) [5G 916.5 1a ee ee ee eee 
Year 17 “5 | 2,019,447.5'- | 7 9707 705/302 ee eee 00 eee 
[Year 182 | 1,842 561.4 2072.7 [20 7139 ae 20 ee eee 
Year 19 5 [yi,555,4113 (|) 614 uel] 10 es ee eee eee) cee 
pear Ne ns SSE ts SS 
a ee re ee a ee 

ae Se ee ee 

a a SR Oe 

FYear 24. |. 905,817.35. |. 035.7 ci ile 16.8) 17s i el OO sine eee 
seine ee 6 See 

pier Be 

ee ed 8 ae 

Oe 

eu a eee 

Beet ei 


Pee ea a ee 
ee eS Oe a 
ae CR eS 
eee ee eee 
| os Rs OSA 0G 4) OR tell ae aig 8 45 Oe eee 


Note: AR6 (sixth oe iar of the eee! Panel on Climate Change); CH4 (methane); CO2 (carbon dioxide); COze 
: (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential); N2O (nitrous oxide). 


509,660 510,716 g 
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Black carbon is a short-lived pollutant with an estimated lifetime of several days to weeks (AMAP 2011, 2015; 

_ IPCC 2021; Paris, Stohl et al. 2009). Black carbon emissions have a positive forcing effect and warm the climate 
_ both in the atmosphere and when deposited on snow or ice (Bond, Doherty et al. 2013; IPCC 2021). The IPCC 

_ (2018a) reports that black carbon emissions must fall by at least 35% across all sectors from 2010 levels by 2050 
_ to limit global warming to 1.5°C (2.7°F). 


Black carbon is a by-product of incomplete combustion. It is removed from the atmosphere through wet and dry 
deposition. Concentrations of black carbon vary depending on the season (AMAP 2015), spatial location 
(Creamean, Maahn et al. 2018), and vertical height in the atmosphere (Creamean, Maahn et al. 2018; Stohl, Klimont 
et al. 2013; Xu, Martin et al. 2017). On Alaska’s North Slope, black carbon can come from international 
transportation sources (Matsui, Kondo et al. 2011; Stohl 2006; Xu, Martin et al. 2017), biomass burning (Creamean, 
Maahn et al. 2018; Stohl 2006; Xu, Martin et al. 2017), shipping (Corbett, Lack et al. 2010; Lack and Corbett 2012), 
oil and gas exploration and production activities (Creamean, Maahn et al. 2018; Stohl, Klimont et al. 2013), and 
residential combustion (Stohl, Klimont et al. 2013). In particular, black carbon emitted from shipping can be 
deposited directly onto sea ice, and ice breakers can deposit black carbon onto the ice pack itself (Brewer 2015). 
Black carbon emitted onto ice and snow can increase melting and exacerbate warming as darker and more absorbent 
land and water surfaces are exposed as a result. With Project construction, black carbon would be emitted as part of 
particulate matter less than 2.5 microns in aerodynamic diameter (PM2,s) emissions from diesel-fired equipment, 
including engines, boilers, heaters, pumping units, and other equipment, such as aircrafts and flares. 
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Black carbon has a strong impact on Arctic regions due to its ability to change the reflective properties of ice and 
snow. When black carbon is deposited on ice or snow, it darkens the ground, decreasing the reflectiveness of the 
surface (the albedo) and warming the surface (+0.13 watts per square meter [W/m7]) (Bond, Doherty et al. 2013). 
Since black carbon emitted in the Arctic has a higher probability of being deposited onto snow or ice, this “snow- 
and ice-albedo feedback effect” is stronger when black carbon is emitted in the Arctic than when it is transported 
from lower latitudes (Sand, Berntsen et al. 2013). Black carbon that is not deposited can increase warming when it 
absorbs solar radiation in the lower troposphere and boundary layer, decreasing cloud cover and leading to 
increased melting, further enhancing the snow- and ice-albedo feedback effect as the surface turns from bright 
snow and ice into darker water. In fact, black carbon has a strong direct radiative effect, meaning it is effective at 
warming the climate through the direct absorption of radiation, and is the component of PM2.s that is most 
effective at absorbing solar energy. For the period 1750 to 2005, Bond, Doherty et al. (2013) estimated the direct 
radiative effect of black carbon to be +0.71 W/m? and the total climate forcing (including cloud, snow, and sea ice 
effects) to be +1.1 W/m?. Black carbon can also affect the formation of clouds and change their radiative 
properties, leading to increased warming (+0.23 W/m’) (Bond, Doherty et al. 2013). When black carbon mixes 
with other pollutants in the atmosphere, a coating can form around the black carbon particle, causing it to grow in 
size. It is predicted that black carbon particles that have reacted with chemical compounds in this way may have 
an increased warming effect (Kodros, Hanna et al. 2018). 


Black carbon can also cool the climate. When black carbon is lofted high into the atmosphere, it can block solar 
radiation from reaching the surface in a process called surface dimming (Flanner 2013; Sand, Berntsen et al. 
2013). Surface dimming also decreases the equatorial-polar temperature gradient, causing less heat to be 
transported to the Arctic from lower latitudes. Black carbon can also increase reflected incoming solar radiation 
by increasing high-altitude clouds that reflect solar radiation. Bond, Doherty et al. (2013) also find that black 
carbon is co-emitted with other pollutants, and these pollutants can reduce the amount of warming caused by 
black carbon alone (-0.06 W/m’). 


The effect of black carbon, although expected to be positive overall, is highly variable and dependent on the 
location and timing of the emissions, the mixing state of the atmosphere, and deposition processes. The complex 
interactions and feedbacks between black carbon and the environment all contribute to the effect of black carbon 
on the arctic climate. 


Black carbon would be emitted by sources and activities under Alternative B. For the Project, black carbon 
emissions were not explicitly quantified; however, black carbon is a component of PM2.s and black carbon 
emissions are included in PM>,; emissions that are quantified in the air quality analysis (Chapter 3.3, Air Quality). 


3.1.3 Alternative C: Disconnected Infield Roads* 


LE eT aT Ce a eh Se aN Ta naan cS cca TA 


SEIS. Annual average GHG emissions (Table E.2.12) are calculated by dividing the Project’s lifetime GHG 
_ emissions by the 30-year Project duration. The annual average emissions for Alternative C shown in Table E.2.12 | 
_ are for gross GHG emissions and do not account for the market substitution effects discussed in Section 3.2.2.2 of | 
_ the SEIS. 


| The annual direct and gross indirect emissions under Alternative C are provided in Table E.2.13 and Table E.2.14, 
_ respectively, for each year of the life of the Project. The annual downstream emissions resulting from the change 
_ in foreign oil consumption over the life of the Project under Alternative C are shown in Table E.2.15. 


_ Direct GHG emissions over the life of the Project calculated with the IPCC AR4 100-year GWPs are 9.3% higher : 
_ than Alternative B due to the increased air travel and two operations centers, 8.8% higher than Alternative D and 
_ 9.2% higher than Alternative E. The annual average direct GHG emissions (844 TMT of COvze per year) over the 
_ 30-year Project life are approximately 2.1% of the 2015 Alaska GHG inventory. The annual average total gross 

| GHG emissions of 9,537 TMT of COze per year constitute approximately 0.1% of the 2019 U.S. GHG inventory 
_ and approximately 0.3% of the U.S. net GHG emissions target for 2030. 


_ When applying the 100-year GWPs from the IPCC AR6, annual direct GHG emissions over the life of the Project | 
_ (846 TMT of COze per year) represent approximately 2.1% of the 2015 Alaska GHG inventory. The annual 
_ average total gross GHG emissions of 9,540 TMT of COze per year represents approximately 0.1% of the 2019 
_U.S. GHG inventory and approximately 0.3% of the U.S. 2030 net GHG emissions target. Thus, when applying 
_ either AR4 or AR6 100-year GWPs, total gross GHG emissions of the Project duration under Alternative C are 
0.8% higher than Alternative B, 0.7% higher than Alternative D, and 3.1% higher than Alternative E. 
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When applying the 20-year GWPs from the IPCC AR6, direct GHG emissions over the 30-year Project life are 
9.0% higher than Alternative B, 8.7% higher than Alternative D, and 9.1% higher than Alternative E. The annual 
_ average direct GHG emissions (861 TMT of COze per year) over the Project life are approximately 2.2% of the 

_ 2015 Alaska GHG inventory. The annual average total gross GHG emissions of 9,588 TMT of COze per year 

| constitute approximately 0.1% of the 2019 U.S. GHG inventory and approximately 0.3% of the U.S. 2030 net 

_ GHG emissions target. Total gross GHG emissions over the Project life under Alternative C calculated with 20- 

_ year AR6 GWPs are 0.8% higher than Alternative B, 0.7% higher than Alternative D, and 3.1% higher than 

_ Alternative E. 


| Over the Project duration under Alternative C, there would be a net increase of up to 73,327 TMT of COze from 
_ the No Action Alternative (Alternative A) to Alternative C, with the highest increase estimated with the 20-year 


_ GWPs. Regardless of the choice of GWPs, the annual average total gross GHG emissions due to the Project under 


_ Alternative C would constitute approximately 0.1% of the 2019 total U.S. GHG inventory and approximately 
- 0.3% of the U.S. 2030 net GHG emissions target. 


Black carbon would be emitted by sources and activities under Alternative C. Although black carbon is not 
explicitly quantified, it is a component of PM2.s, and PM2.s emissions would be approximately 19% greater under 
Alternative C than Alternative B (see Appendix E.3B, Air Quality Technical Support Document). Therefore, it is 
anticipated that black carbon emissions would also be greater under Alternative C than Alternative B, and the 
effects of black carbon on the environment would increase under Alternative C relative to Alternative B. 


Table E.2.12. Annual Average Gross Greenhouse Gas Emissions under Alternative C (thousand metric 
tons per year)* 


CQre 
GHG COre COvre 
‘Emissions See cm ND ann mis! (100-year AR6 GWP) | (20-year AR6 GWP) 


[Tota [9,488 | 0.9117 [0.0911 9,537 A 


Note: AR4 (fourth assessment report of the Intergovernmental Panel on Climate Change [IPCC]); AR6 (sixth assessment report of the IPCC); CH, (methane); 
CO; (carbon dioxide); COze (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential; N2O (nitrous oxide). Year 0 only included 1 
month of construction activity, and thus this year was excluded from the average annual emissions. 

“Total values may have small differences due to rounding. 


Table E.2.13. Annual Direct Greenhouse Gas Emissions (metric tons) under Alternative C* 


AR6 100-year GWPs) AR6 20-year GWPs 
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Year 19 
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AR6 100-year GWPs (AR6 20-year GWPs 
968,415 
968,415 
968,415 
968,415 
968,415 


i. Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CHs (methane); COz (carbon dioxide); CO2e 
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Table E.2.14. Annual Gross Domestic Indirect Greenhouse Gas Emissions (metric tons) under Alternative 


C* 

AR6 100-year GWPs (AR6 20-year GWPs) 
Pe re eo ee i 715 
eee a eC On el 00 ere aia P69 
po Oo e200 2,209 
OS EO? 6 a OC 2 1,689 
eee eer es iO a se a 3900 
25,137,638 
27,435,963.3 27,673,619 
24,476,621.1 
21,705,397.1 
19,147,419.1 
16,956,306.9 
15,235,555.0 
13,737,597.7 
5,128,407 
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Table E.2.15. Annual Downstream Greenhouse Gas Emissions (metric tons) Resulting from the Change in 
Foreign Oil Consumption under Alternative C* 


Pe ae, AR6 100-year GWPs AR6 20-year GWPs 

Pos Xe ee ee 
Poy ee ee Oe Le 
PYear3 [70,0 |S 0.0 Fi 0 00 See cea a eee 
RC ys Kee ee a 
Ree es OK 
6,115,578.2 
5,297,604.4 
4,543,065.7 
888,743.6 


f Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CH4 (methane); CO2 (carbon dioxide); CO2e 


bo cccSSecccasaceccadecststecssescudsedaccosnsteustcoseseccbbcecd Genesc cls sasabeausacedesddcretenstowestecusccucccsusecesuernaauacssenccentcctenadscauscutstwersnstesssudersddccessedcccsuuecsre sbededsneabassessudeasasstossuwcestewvensest bee beseeccutecencersesuscalanuesduenscunsbaceaanavassnauessenewesseusecsbsees basuasusssecuucenussesssbesseeeneeeweanaeet 


pesepnarsnnnnncannnceesenc=asnspnannensnseccsnununnnenneanavaannnsnnacsnnnansnonnmannnnranennhannasenarraanennscorenrennnnenenssrnnns onecnonsesnonnnanccnnenrsnnnnAsnansnntn>wsnsn~oneene*rennnennsns ann aeeaOnPnehS=hs=nns © S@EPna eben eanrnn Danas =SDESS=OE&== SEAS =SS5 oS S6SSSRSSS SSE RSESSSSSS=S-sh= SSSSS 5 SEEESSSSESSS SOS EA SSSEESEEASSIOS: SSSSSSSREASSEEESSSSERSSESES 


_ As mentioned in Section 2.2 of this appendix and explained in more detail in Chapter 2.0, Alternatives, 
_ Alternative D would have a 31-year Project lifetime rather than the 30-year Project lifetime for Alternatives B and 
_C. Alternative D GHG emissions estimated over the 31-year Project lifetime are shown in Section 3.2.2.3 of the 
_ SEIS. The average annual total gross emissions under Alternative D are provided in Table E.2.16; these emissions 
_ do not account for the market substitution effects discussed in Section 3.2.2.2 of the SEIS. 


| The annual direct and gross indirect emissions under Alternative D are provided in Table E.2.17 and Table 
| E.2.18, respectively, for each year of the life of the Project. The annual downstream emissions resulting from the 
_ change in foreign oil consumption over the life of the Project under Alternative D are shown in Table E.2.19. 


_ When applying the 100-year GWPs from the IPCC AR4, direct GHG COvze emissions over the 31-year Project life 
_ of Alternative D are 0.5% higher than Alternative B and 0.4% higher than Alternative E primarily due to 
_ increased air travel. The annual average direct GHG emissions (751 TMT of COvze per year) over the Project 
_ duration are approximately 1.9% of the 2015 Alaska GHG inventory. The annual average total GHG emissions of : 
9,164 TMT of COze per year constitute approximately 0.1% of the 2019 U.S. GHG inventory and approximately 
0.3% of the U.S. 2030 net GHG emissions target. When applying the 100-year GWPs from the IPCC AR6, direct 

_ GHG COvxe emissions over the Project life are 0.5% higher than Alternative B and 0.4% higher than Alternative 
_ E. The annual average direct GHG emissions (752 TMT of COve per year) over the Project life are approximately 
1.9% of the 2015 Alaska GHG inventory. The annual average total GHG emissions are 9,166 TMT of COze per | 
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_ year, which represent approximately 0.1% of the 2019 U.S. GHG inventory and approximately 0.3% of the U.S. 
_ 2030 net GHG emissions target. Thus, when applying the 100-year GWPs from either AR4 or AR6, total gross 

_ GHG emissions over the Project life under Alternative D are 0.05% higher than Alternative B, 0.7% lower than 
_ Alternative C, and 2.3% higher than Alternative E. 


_ When applying the 20-year GWPs from the IPCC AR6, direct GHG CO:e emissions over the Alternative D 

_ Project life are 0.5% higher than Alternative B and 0.4% higher than Alternative E. The annual average direct 

_ GHG emissions (767 TMT of COve per year) over the 31-year Project life are approximately 1.9% of the 2015 

_ Alaska GHG inventory, and the annual average total GHG emissions of 9,212 TMT of COvze per year constitute 
_ 0.1% of the 2019 U.S. GHG inventory and approximately 0.3% of the U.S. 2030 net GHG emissions target. Total | 
_ gross GHG emissions over the Project duration under Alternative D calculated with 20-year IPCC AR6 GWPs are 
_ 0.04% higher than Alternative B, 0.7% lower than Alternative C and 2.3% higher than Alternative E. 


_ Over the 31-year life of the Project under Alternative D, there would be a net increase of up to 71,282 TMT of 
_ COvze from the No Action Alternative (Alternative A) to Alternative D, with the highest increase estimated using 
_ the 20-year IPCC AR6 GWPs. Regardless of the choice of GWPs, the annual average total gross GHG emissions 
_ due to the Project under Alternative D represent approximately 0.1% of the 2019 total U.S. GHG inventory and 


Black carbon would be emitted by sources and activities under Alternative D. Although black carbon is not 
explicitly quantified, it is a component of PM2.s, and PM2.s emissions would be approximately 8% greater under 
Alternative D than Alternative B and emissions under Alternative D would be approximately 10% less than 
Alternative C (see Appendix E.3B, Air Quality Technical Support Document). Therefore, it is anticipated that 
black carbon emissions would be greater under Alternative D than Alternative B but reduced relative to 
Alternative C. Similarly, the effects of black carbon on the environment described in Section 3.2.1, Affected 
Environment, would increase under Alternative D relative to Alternative B. 


Table E.2.16. Annual Average Greenhouse Gas Emissions under Alternative D (thousand metric tons per 
ear)* 


CQvze 
GHG COre COrxe 
Raisin ean (100-year AR6 GWP) | (20-year AR6 GWP) 
Direct 743 0.2834 0.0017 751 752 767 y 
Indirect 8.373 0.595 0.086 8,413 | 8.414 8,445 
Total’ 9,116 0.8779 0.0879 9,164 9,166 9,212 


Note: AR4 (fourth assessment report of the Intergovernmental Panel on Climate Change [IPCC]); AR6 (sixth assessment report of the IPCC); CH, (methane); 
CO; (carbon dioxide); CO ze (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential; NO (nitrous oxide). 
“Total values may have small differences due to rounding. 


Table E.2.17. Annual Direct Greenhouse Gas Emissions (metric tons) under Alternative D* 


j CO 4 COzre COzre 
AR6 100-year GWPs AR6 20-year GWPs 
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: Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CHs (methane); CO2 (carbon dioxide); COze 
: (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential); N2O (nitrous oxide). 


Table E.2.18. Annual Gross Domestic Indirect Greenhouse Gas Emissions (metric tons) under Alternative 
D* 


Year 8 

Year 9 

Year 11 

Year 12 17,040,082 

Year 13 
‘ 874.5 : 5319) ; 


AR6 100-year GWPs AR6 20-year GWPs 
PYear 2 708.8. | 0.1) es OO eS Le ee ee ee 
{Year 200 52,611.00 01 2 000 Gl a eee 
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: Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CH4 (methane); COz (carbon dioxide); COze 
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Table E.2.19. Annual Downstream Greenhouse Gas Emissions (metric tons) Resulting from the Change in 
Foreign Oil Consumption under Alternative D* 


(AR6 100-year GWPs AR6 20-year GWPs 

Me OO er Or Mn OO ur ee mi Ole oa hole 
op EIS SNES OE ACT OS es 0 
Ore ch U0 din OO Kee be tOKe en Te OO Ou | 
ee OS OO ee OC en ce O ay tee lea On | 
ee OO OO Ti Oar i So Oy 0 
Sere Mi OO ecu a CO a he Of 
6,898,595 
6,148,662.6 6,181,365 
4,869,626 
3,949,761 
3,526,747 
2,618,229 
2,451,442 
1,973,874 
1,731,761 

Oe OA ES ie ee 1,305,853 1,308,560 

985,486 

693,421 

610,393 

668,835 

608,463.3 
es a ica ee ae 503,884 565,053 


Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CH4 (methane); CO2 (carbon dioxide); COze 


eae ee aa ae EAC a eee TEI Lada eens AIC TNR ee CLI Coa eneb ee ne cee tata e dee eWEEEeEbta cece EEenewseaetNecsesenabaTeeda=sosseasau teed -terase~CeeeeueTeb EaSadewnereeeeeananteeeuNers=eecannnsnonsnsnnwewwawWeverseenenn enn noaueen ween nrnrernaen 


a ee ae OCI ee a ee err ener erie ane ticker ac Sam NU ric Cr esa rd appt nt nied Tah an abe nga uso ee aa 


| have a 30-year Project life. 


_ Project facilities proposed for Alternative E are generally the same as Alternative B, with the exception that 
_ Alternative E would not include construction of drill site BT4, and drill site BT2 would be located farther north at — 
| the coordinates for BT2 in Alternative B. BTS would be located east of the location proposed for other action 
alternatives, which would also reduce the length of the BTS road and infield pipelines. 


| Alternative E GHG emissions estimated over the 30-year Project life are shown in Section 3.2.2.3 of the SEIS. : 
_ The average annual total gross emissions under Alternative E are provided in Table E.2.20; these emissions do not | 
_ account for the market substitution effects discussed in Section 3.2.2.2, Alternative A: No Action. | 


| The annual direct and gross indirect emissions under Alternative E are provided in Table E.2.21 and Table E.2.22, | 
_ respectively, for each year of the life of the Project. The annual downstream emissions resulting from the change 
_ in foreign oil consumption over the life of the Project under Alternative E are shown in Table E.2.23. 


| When applying the 100-year GWPs from the IPCC AR4, direct GHG COze emissions over the 30-year Project life : 
_ of Alternative E are 0.1% higher than Alternative B. In contrast, the indirect gross GHG emissions (as well as | 
_ total GHG emissions) are lower under Alternative E than Alternative B because total oil production would be 
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lower under Alternative E and total emissions are dominated by indirect emissions. This is true when applying the | 
_ other GWPs also as discussed below. The annual average direct GHG emissions (773 TMT of COsze per year) 
_ over the Project duration are approximately 1.9% of the 2015 Alaska GHG inventory. The annual average total 

_ GHG emissions of 9,253 TMT of COze per year constitute approximately 0.1% of the 2019 U.S. GHG inventory | 
_ and approximately 0.3% of the U.S. net GHG emissions target for 2030. The 100-year GWPs from the IPCC AR6 | 
_ direct GHG COsze emissions over the Project life are 0.1% higher than Alternative B. The annual average direct 
_ GHG emissions (774 TMT of COze per year) over the Project life are approximately 1.9% of the 2015 Alaska 

_ GHG inventory. The annual average total GHG emissions are 9,255 TMT of COsze per year; they represent 

_ approximately 0.1% of the 2019 U.S. GHG inventory and approximately 0.3% of the U.S. net GHG emissions 

_ target for 2030. Thus, when applying the 100-year GWPs from either AR4 or AR6, total gross GHG emissions 

_ over the Project life under Alternative E are 2.2% lower than Alternative B, 3.0% lower than Alternative C, and 

_ 2.3% lower than Alternative D. 


_ When applying the 20-year GWPs from the IPCC AR6, direct GHG COze emissions over the Project life are 
_ 0.01% higher than Alternative B. The annual average direct GHG emissions (789 TMT of COze per year) over the | 
_ 30-year Project life are approximately 2.0% of the 2015 Alaska GHG inventory, and the annual average total | 
_ GHG emissions of 9,301 TMT of COve per year constitute 0.1% of the 2019 U.S. GHG inventory and 

_ approximately 0.3% of the U.S. net GHG emissions target for 2030. Total gross GHG emissions over the Project 

| duration under Alternative E calculated with 20-year IPCC AR6 GWPs are 2.2% lower than Alternative B, 3.0% 

: lower than Alternative C, and 2.3% lower than Alternative D. 


_ Black carbon would be emitted by sources and activities under Alternative E. Although black carbon is not 
_ explicitly quantified, it is a component of PM2.s and PM2.s emissions under Alternative E would be approximately | 
_ comparable to (0.005% higher than) Alternative B while emissions under Alternative E would be approximately | 
- 16% less than Alternative C and approximately 6% less than Alternative D (see Appendix E.3B, Air Quality 

_ Technical Support Document). Therefore, it is anticipated that black carbon emissions under Alternative E would 

_ be comparable to Alternative B but less than Alternatives C and D. Similarly, the effects of black carbon on the 

_ environment under Alternative E, described in Section 3.2.1, Affected Environment, would be comparable to 

| Alternative B. 


Table E.2.20. Annual Average Greenhouse Gas Emissions under Alternative E (thousand metric tons per 


year)* 
GHG CQvre CQvre 
Emissions | ©? oe No (20-year AR6 GWP 
Direct 765 0.2816 0.0018 | 789 
indirect | 8,439 0.599 0.087 | 8,512 
Total? 9,204 0.8807 0.0886 


Note: AR4 (fourth assessment report of the Intergovernmental Panel on Climate Change [IPCC]); AR6 (sixth assessment report of the IPCC); CH, 

: (methane); CO, (carbon dioxide); COze (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential; N2O (nitrous oxide). Year 0 
: only included 1 month of construction activity and thus this year was excluded from the average annual emissions. 

: *Total values may have small differences due to rounding. 


ct Year CO2 CH4 N20 CQvre CQre 
AR6 100-year GWPs AR6 20-year GWPs 


Proje 
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AR6 100-year GWPs AR6 20-year GWPs) 
895,218 
895,218 
895,218 
866,877.3 895,218 
866,877.3 895,218 
895,218 
895,218 
8 CO Ce a eee 895,218 


a 


Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CH, (methane); COz (carbon dioxide); CO2e 
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Table E.2.22. Annual Gross Domestic Indirect Greenhouse Gas Emissions (metric tons) under Alternative 


E* 

(AR6 100-year GWPs (AR6 20-year GWPs 
er ee es I sy 
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11,596,365 
10,483,597 
8,724,446 
2,893,204 
2,443,808 
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Table E.2.23. Annual Downstream Greenhouse Gas Emissions (metric tons) Resulting from the Change in 
Foreign Oil Consumption under Alternative E* 
Project Year CO2 CHa N20 
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3,769,803.6 148.8 28.5 : 
3,433,887.8 135.5 25.9 : 
3,010,177.7 118.8 ; 3,019,928 3,026,188 


Note: AR6 (sixth assessment report of the Intergovernmental Panel on Climate Change); CH4 (methane); CO2 (carbon dioxide); COze 
: (carbon dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential); N2O (nitrous oxide). 


3.1.6 Module Delivery Options 

Project lifetime and annual average direct GHG emissions from module delivery options alone are shown in Table 
E.2.24 for Option | (Atigaru Point Module Transfer Island), Option 2 (Point Lonely Module Transfer Island) and 
Option 3 (Colville River Crossing). Note that emissions from Option 3 vary based on the action alternative it is 
paired with for analysis. Table E.2.6 also provides the differences between Options | and 2 from Option 3. 
Annual average GHG emissions for module delivery options are calculated by dividing the Project lifetime GHG 
emissions by the expected duration of module delivery emissions, which is 6 years. Direct GHG emissions from 
Option 2 are more than twice the emissions from Option | because vehicles would travel a longer distance to 
reach Point Lonely. Direct GHG emissions from Option 3 are considerably less than Options 1 and 2 (under all 
action alternatives) because Option 3 would make use of the existing Oliktok Dock and construct the least amount 
of new infrastructure to support sealift module delivery. Total GHG emissions for the Project would be the sum of 
the selected alternative and the selected module delivery option. 


Black carbon would be emitted by sources and activities as part of all module delivery options. Although black 
carbon is not explicitly quantified, it is a component of PM>2.s, and PM2.s emissions would be greatest under 
Option 2 and lowest under Option 3 (under all action alternatives). Therefore, it is anticipated that black carbon 
emissions would also be greatest under Option 2 and lowest under Option 3 (under all action alternatives), and the 
effects of black carbon on the environment described in Section 3.1.2.3 of this appendix, would be greatest under 
Option 2 and lowest under Option 3 (under all action alternatives). 
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Table E.2.24. Direct Greenhouse Gas Emissions Associated with Module Delivery Options (thousand 
metric tons) 
GHG Emissions Total Annual Total Annual Total Annual 
CQvre Average COze COvre Average COze COvze Average COze 
(100-year (100-year AR4 (100-year (100-year AR6 (20-year (20-year AR6 


AR4 AR6 AR6 
GWP) GWP) GWP) 
| Option 1: Atigaru Point MTI 140 23 140 | 23 141 23 
Option 2: Point Lonely MTI 1 341 ay) 341 et 342 ey) 
Option 3: Colville River Crossing 40 7 40 7 A000 | 7 
— Alternatives B, C, and E 
Option 3: Colville River Crossing 43 a 43 ‘| 43 7 
| — Alternative D 
Option | minus Option 3 Tah 100 Vy 100 Lg 101 17 
Alternatives B, C, and E) 
POption 1 minus Option 3 oF 16 97 16 oF 16 
(Alternative D) _ 
Option 2 minus Option 3 301 50 301 50 BO) ea 50 
(Alternatives B, C, and E) 
Option 2 minus Option 3 298 50 298 50 298 50 
Alternative D 


Note: AR4 (fourth assessment report of the Intergovernmental Panel on Climate Change [IPCC]); AR6 (sixth assessment report of the IPCC); COze (carbon 
dioxide equivalent); GHG (greenhouse gas); GWP (global warming potential); MTI (module transfer island). 


3.2 Climate Test Tool* 


RATS eT eh ae a ne] 


_ During the public comment period on the Draft SEIS, BLM received comments about a new tool, the climate test 
| tool, for evaluating the Project’s potential impact on U.S. climate commitments and goals (Earthjustice, Natural 
_ Resources Defense Council et al. 2022). The climate test tool was developed by NRDC (Bustamante, Alexander 
et al. 2022); draft unpublished manuscript submitted for review to ‘Climate Policy’) to assist agencies in 
_ determining the significance of a project’s GHG emissions. An overview of the tool is shown in Figure E.2.1 
_ (Bustamante, Alexander et al. 2022)). A rearrangement of the equation shown on the right hand side of Figure 

_ E.2.1 indicates that if the emissions intensity of a project (proj ect lifecycle emissions divided by energy supplied 
_ by project) is greater than the emissions intensity of the remaining carbon budget (remaining emissions in budget 
aes by remaining fossil energy demand if all existing fossil projects proceed without constraint), then the 


Climate Test Decision Metric: 
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impact 


Figure E.2.1. Conceptual framework, data sources, and input metrics for the climate test decision metric to 
evaluate emissions significance for individual fossil fuel projects Figure ource: Bustamante, 
Alexander et al. (2022)* 


Appendix E.2A Climate and Climate Change Page 29 


Willow Master Development Plan Final Supplemental Environmental Impact Statement 


PR nv np eC Eee ene ee 


: Sionif ek a cance: A legal term Of art derived from the statutory context of the National 

_ Environmental Policy Act (NEPA) (42 U.S.C. § 4321 et seq.) and similar state-level 

_ environmental review statutes. The term references, inter alia, the key criterion determining the 
_ level of environmental review a decision-making agency is required to undertake. 40 CFR. § 
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There are no specific NEPA guidelines to determine the significance of a particular quantity of GHG emissions. 
_ The climate test tool offers one way to assess significance for determining the level of environmental review 

_ required for a particular project, i.e., a project which is found significant under this tool would require a higher 

_ level of environmental review (namely an EIS) than one which is not. The results of the climate test tool for the 
_ Willow MDP are disclosed below as reported by the NRDC. 


_ The climate test tool was applied by the NRDC to the Willow MDP Project using GHG emissions and production | 
_ data from the Draft SEIS. The BLM has incorporated the results of this climate tool assessment of the Project 
_ under the U.S. 2050 net-zero CO2 emissions scenario into the SEIS, as this is a U.S. goal established by Executive | 
_ Order 14057, “Catalyzing Clean Energy Industries and Jobs Through Federal Sustainability” (86 FR 70935). A 
_ summary of the methods of the tool and results for the Project are provided below and the results are presented in 

_ the analysis of the action alternatives in Sections 3.2.2.3 to 3.2.2.6 of the SEIS. 


_ The climate test tool first calculates the “emissions impact” of the Project as the ratio of the total lifecycle Project 

_ emissions to the total remaining emissions under the US 2050 net-zero goal. The total remaining emissions are 

_ calculated as the difference between budgeted emissions under the climate goal and committed emissions from 

_ existing fossil fuel resources. In the analysis of the Project, NRDC applied the annual average direct and indirect 

_ domestic gross CO2 emissions under each action alternative in Tables E.2.3 to E.2.5 of Appendix E.2A of the 

_ Draft SEIS. The emissions were applied to each year of the life of the Project (i.e., 2023 to 2052 under 
_ Alternatives B, C and E, and 2023 to 2053 under Alternative D). Emissions from Module Delivery Option 1 were — 
_ applied; the other module delivery options were not evaluated. For the total remaining U.S. emissions under the 
_ 2050 net-zero CO: emissions scenario, NRDC used the difference between the estimated total remaining CO2 

_ emissions from Princeton University’s Net Zero America study (Larson, Greig et al. 2021) and committed 

_ emissions from existing fossil fuel infrastructure estimated by Tong, Zhang et al. (2019). 
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_ The tool then calculates the Project’s “energy contribution” as the ratio of the total energy supplied by the Project 
_ to the total unmet fossil energy demand under the US 2050 net-zero goal. The energy supplied by the Project was 
_ calculated by NRDC using the annual oil production in Appendix D.1 of the Draft SEIS. The total unmet fossil 

_ energy demand under the US 2050 net-zero goal was estimated using modeling of Larson et al. (2021) and data 

_ on existing fossil energy supplies from the USEIA (2019). 


_ The results of NRDC’s climate test analysis for the U.S. 2050 net-zero scenario for the gross domestic emissions 

_ of the Project under each action alternative are provided in Table E.2.25 (adapted from Earthjustice, NRDC and 

_ The Wilderness Society 2022). The Climate Test tool analysis does not specifically address the effect of Project 

_ production on foreign consumption and emissions, and thus the latter is not reported here. The results of NRDC’s 

_ climate test analysis using gross domestic Project emissions for the U.S. 2050 net-zero scenario indicate that the 

_ effect of the Project GHG emissions under all action alternatives is greater than the energy contribution (i.e., the 

_ ratio of effect to contribution is greater than 1) with Alternative B having a slightly lower ratio than the other 

_ action alternatives (i.e., Alternative B is less than Alternative E is less than Alternative C is less than Alternative 
_D). The ratio here (“emission significance metric result” in Table E.2.25) is the ratio of the Emissions Impact 

_ Submetric to the Energy Contribution Submetric, with a ratio greater than one defined as being significant, as 
_ reported by NRDC for all action alternatives for the Willow MDP. This is consistent with BLM’s development of | 
_an EIS for the Willow MDP project which is the appropriate level of environmental review under all action 
_ alternatives. 


| When considering the net emissions, i.e., the difference between the No Action and any action alternative, the 
_ climate test analysis assumes that the emissions and energy produced by substitute energy sources are only from 
existing infrastructure sources. The net emissions analysis ee NRDC found approximately 1% reduction in the 
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Table E.2.25. Results from the National Resource Defense Council’s Climate Test Analysis of the Gross 
Greenhouse Gas Emissions from the Project for a Net-zero 2050 Scenario* 


| Alternative Project Total | Project Total | Emissions Impact | Energy Contribution Emission 
CO: Emissions Energy Submetric Submetric Significance 
GtCO2 Supplied % % Metric Result (“% 
A: No Action 0 0 0 


Alternative 


B: Proponent’s Project 0.283 3,745 1.186 0.445 2.662 


C: Disconnected Infield 0.285 3.745 1.195 0.445 2.683 
Roads 


D: Disconnected Access 0.283 3.745 1.210 0.437 D966 


E: Three-Pad Alternative U2EF 3.649 1.159 0.434 2.671 
_|(Fourth Pad Deferred) | eee 


| Source: Adapted from Earthjustice, NRDC and TWS (2022). 
| Notes: CO> (carbon dioxide); EJ (exajoule); GtCO2 (gigaton CO2). Results for Module Delivery Option | plus each action alternative. 
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1. Background 


The Energy Substitution Model (EnergySub) is a tool that enables BLM to compare unobservable long-run 
market conditions with and without potential mineral production under onshore development scenarios using 
baseline energy projections developed by the U.S. Energy Administration (EIA).' The EnergySub model is not a 
national forecasting model and it was not designed to be a replacement for EIA’s integrated modeling systems or 
the annual energy projections developed by EIA’s Office of Energy Analysis. The BLM developed this model to 
assess potential market responses associated with onshore oil, gas, and coal related management actions, 
including possible substitution between various energy sources and changes in energy prices and consumption, 
given market conditions projected by the EIA. EnergySub was adapted from Bureau of Ocean Energy 
Management’s (BOEM) Market Simulation Model (MarketSim), which assesses potential market impacts of 
development of offshore oil and gas resources along the Outer Continental Shelf (OCS).* While EnergySub 
includes substantive updates to enable the model to simulate changes in onshore mineral development, it retains 
much of the overall structure and functionality of MarketSim. 


EnergySub was used to assess the potential energy market impacts attributable to onshore oil production from the 
Willow Master Project. Estimates of displaced energy substitutes and potential effects on foreign oil demand 
produced through EnergySub simulations were used as inputs in the supplemental analysis of the Project’s 
potential GHG emissions. The BLM used EnergySub to conduct a quantitative analysis for this SEIS because of 
on the specific production aspects of Willow and BLM's prior use of BOEM’s MarketSim model in the original 
Willow EIS. A quantitative analysis of market effects, however, is one of several approaches to assessing the 
impacts of BLM management decisions. 


2. Model Overview 


EnergySub is an excel-based partial equilibrium model that uses a series of supply and demand equations with a 
set of assumed long-run elasticities and partial adjustment parameters to create a mathematical representation of 
U.S. energy markets. The model simulates end-use domestic consumption of oil, natural gas, coal and electricity 
in four sectors (residential, commercial, industrial and transportation); production of primary energy fuel sources; 
and the transformation of renewable and nonrenewable fuel sources into electricity. The model primarily 
represents U.S. energy markets but captures interactions with foreign markets through its mathematical 
representation of a global oil market with aggregated foreign supply and demand and its inclusion of domestic 
imports and exports. 


EnergySub relies on baseline long-run energy projections developed by the EIA to calibrate its supply and 
demand equations to an initial market equilibrium. Production schedules for onshore federal oil, natural gas, and 
coal from a Reasonable Foreseeable Development Scenario (RFD) are used to shock the supply side of 
EnergySub’s initial market equilibrium, causing the model to solve its system of equation for new equilibrating 
prices for energy and energy sources in each year of the production scenario based own and cross price 
elasticities. > Solving for new prices yields equilibrium quantities of energy and energy sources supplied and 
demanded, accounting for substitution between energy fuel sources. | 


' The EIA is the statistical and analytical branch of the Department of Energy and operates within the U.S. Federal statistical 
system as the single federal government authority on energy statistics. Their mandate is to collect, analyze, and disseminate 
energy information to inform and promote policymaking, efficient markets, and public understanding of energy and its 
interactions with the economy and the environment. As part of this mandate, EIA’s Office of Energy Analysis develops and 
maintains the National Energy Modeling System and World Energy Projection System in order to develop annual energy 
projections widely used by Members of Congress, industry participants, government agencies, and the public. 

> See Industrial Economics, Inc. (2017). 

3 EnergySub extrapolates baseline energy projections through the life of the production scenario when the modeled time 
period extends beyond the AEO and IEO 2050 baseline projections. 
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3. Model Framework 


As mentioned above, EnergySub uses a series of equations with assumed long-run supply and demand elasticities 
and partial adjustment parameters to depict energy market. These elasticities and adjustment parameters facilitate 
the market equilibrating process that moves the simulated energy market from observable short-run conditions 
towards long-run equilibrium conditions in each year of the simulation. While these long-run conditions cannot be 
directly observed, they can be inferred from short-run market conditions and the model’s underlying parameters. 
The following sections outline EnergySub’s supply and demand equations and discuss how the model 
equilibrates. 


4. Oil Market 


EnergySub models a simplified world oil market with sector detail for the domestic market, a single supply 
equation for foreign production, and a small number of demand equations for foreign consumption. While 
EnergySub can distinguish to some degree where oil in the domestic market is produced (i.e., AK onshore, AK 
offshore, lower-48 onshore, lower-48 offshore), the foreign oil market is a single market made up of all oil 
consumed and produced outside of the U.S.. The estimation of impacts to foreign submarkets is currently beyond 
the modeling capabilities of EnergySub. 


The equations that follow below illustrate how EnergySub estimates U.S. oil demand, foreign oil demand, U.S. oil 
supply, foreign oil supply, oil imports delivered to the U.S. by tanker, U.S. crude oil exports, and U.S. exports of 
refined petroleum products. These equations estimate supply and demand for oil by the residential, commercial, 
industrial and transportation sectors. Oil use for electricity generation is represented elsewhere in the model’s 
electricity module. 


4.1 U.S. Oil Demand 
Qpoit = Aoi,t * Pie : | | se + (1 — Ypoi) Qpoit-1 
j 


for each U.S. end-use sector i; and j = g (gas), c (coal), and e (electricity) where: 


QOpoistepresents the quantity of oil demanded in sector i at time f, 

Aoitis a constant calibrated to the AEO market projections, 

Pris the price of oil at time f, 

Noi is the long-run price elasticity of oil demand in sector /, 

P;, is the price of energy source / at time f, 

Noji is the long-run elasticity of demand for oil with respect to the price of energy source j in sector J/, and 

Ypoi is the rate at which demand for oil in sector i adjusts.’ 
The four U.S. end-use sectors i are residential, commercial, industrial, and transportation. To estimate cross-price 
effects in the industrial and other sectors, EnergySub uses a single weighted average minemouth price of coal 
(instead of the separate regional coal prices described in Section 7 below).° 


4.2 Foreign Oil Demand 


Qnox,t —- Aoxt ; ee rs eI on: Ypox) Qpoxt-1 
Where: 
Opoxs represents the quantity of foreign oil demand at time ¢, 
Aox1 iS a constant calibrated to the AEO market projections, 


4 Note that this deviates from standard notation used in the empirical literature on demand and supply estimation by using 
gammas to represent adjustment rather than persistence. 
5 The model uses the weighted average price of coal, using industrial sector consumption as weights. 


Appendix E.2B - EnergySub 5 


Willow Master Development Plan Final Supplemental Environmental Impact Statement 


P.,ris the price of oil at time ¢, 
Nox is the long-run price elasticity of foreign oil demand, and 
YDox is the rate at which non-U.S. oil demand adjusts. 


Foreign oil demand is strictly a function of the oil price, and no other prices, domestic or foreign. EnergySub 
specifies three categories of foreign oil demand: (1) foreign demand for U.S. crude oil, (2) foreign demand for 
U.S. refined products, and (3) foreign demand for foreign oil. The model assumes that these three categories are 
mutually exclusive. 


4.3 U.S. Oil Supply 
Qsout = Bout * Pee alec Ysou)Qsou,t-1 


for each domestic oil source uv = lower 48 onshore non-tight oil, lower 48 onshore tight oil, lower 48 offshore, 
Alaska onshore, Alaska offshore, biofuels, natural gas plant liquids, other, or rest of world; where: 


QOsout represents the quantity of oil supplied from U.S. source wu at time f, 
Bout is a constant calibrated to the AEO market projections, 

P,,xis the price of oil at time f, 

Nou is the long-run elasticity of oil supply from source u, and 

Ysou 1S the rate at which U.S. oil supply u adjusts. 


Consistent with the EIA classification, the term “oil” includes all liquid fuels that are close substitutes for 
petroleum products (e.g., biofuels). 


4.4 Foreign Oil Supply 


1) 
Qsoy,t = Boyt ; Pee la Ysoy) Qsoy,t-1 


Where: 
Osoy,r represents the quantity of non-U.S. oil supplied at time 1, 
Boy is a constant calibrated to the AEO market projections, 
P.1s the price of oil at time ¢, 
Noy 1s the long-run elasticity of non-U.S. oil supply, and 
Y soy is the rate at which non-U.S. oil supply adjusts. 


Foreign oil supply is estimated in EnergySub’s equilibrating equations as a separate value that represents tanker 
imports and pipeline imports combined, consistent with AZO reporting. 


4.5 Oil Imports Delivered via Pipeline 


EnergySub uses the equations outlined above to find changes in oil market consumption, production, and prices 
under a given development scenario. The model’s calculation for oil imports from Canada is similar to the foreign 
oil supply formula except with its own parameter, elasticity, and adjustment rate. 
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Qsoct = Boct ’ pe + (1 — Ys0c)Qsoc,t-1 
Where: 


Osoc. represents the quantity of Canadian pipeline oil imports supplied at time ¢, 
Boer 18 a constant calibrated to the AEO market projections, 

P.,11s the price of oil at time ¢, 

Noc is the long-run elasticity of Canadian pipeline oil imports, and 

Ysoc is the rate at which the supply of Canadian pipeline oil imports adjusts. 


4.6 U.S. Crude Oil Exports 


As described above, EnergySub models oil as a global market with supply (i.e., production) and demand (i.e., 
consumption) specified separately for the U.S. and the rest of the world. To facilitate the estimation of changes in 
oil exports, EnergySub’s demand equations specify the three categories of foreign demand identified above: (1) 
foreign demand for U.S. crude oil, (2) foreign demand for U.S. refined petroleum products, and (3) foreign 
demand for foreign oil. The first of these items represents U.S. crude oil exports. Therefore, to estimate the impact 
of a given BLM development scenario on U.S. crude oil exports, EnergySub calculates the difference between 
foreign demand for U.S. crude oil between the development scenario and the AEO baseline projections. 


4.7 U.S. Exports of Refined Petroleum Products 


EnergySub estimates U.S. exports of refined petroleum products based on the specification of foreign demand for 
refined petroleum products in the model’s equilibrating equations.° For a given development scenario, the change 
in U.S. refined petroleum product exports is equal to the estimated change in foreign demand for U.S. refined 
petroleum products. This approach is similar to that outlined above for U.S. exports of crude oil, which 
EnergySub estimates based on the change in foreign demand for U.S. crude oil. 


5. Natural Gas Market 


EnergySub represents the U.S. natural gas market with exports and imports. This stands in contrast to the oil 
market, which EnergySub simulates as a global market due to the relatively low cost of transporting oil and the 
large volume of oil traded on international markets. Natural gas use for electricity generation is represented 
elsewhere in the electricity section of the model. The equations that follow specify EnergySub’s estimation of 
U.S. natural gas demand, demand for U.S. natural gas exports, and U.S. natural gas supply. 


5.1 U.S. Natural Gas Demand 


es nt 
Qngi,t Sees Agit ; Bees ; alee sh GR Yogi) Qngit-1 
j 


for each U.S. end-use sector i; and j = 0 (oil), c (coal), and e (electricity) where: 


QOpgis represents the quantity of natural gas demanded in sector / at time f, 
Agi: is a constant calibrated to the AEO market projections, 

P., is the price of natural gas at time ¢, 

Nzi is the long-run price elasticity of natural gas demand in sector 7, 

P;, is the price of energy source / at time f, 


6 As noted above, this category of foreign demand represents one of three included in the model. The other two categories are 
foreign demand for U.S. crude oil and foreign demand for foreign oil. 
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Neji is the long-run elasticity of demand for natural gas with respect to the price of energy source / in 
sector /, and 
Ybgi is the rate at which demand for natural gas in sector i adjusts. 


The U.S. natural gas demand sectors represented in EnergySub include the residential, commercial, industrial, and 
transportation sectors. As in the oil market, EnergySub uses a single weighted average minemouth price of coal 
instead of separate regional coal prices to estimate cross-price effects in the industrial sector. 


5.2 Demand for U.S. Natural Gas Exports 


1 
Qpgx,t a Agx,t ‘ ie stn (Glee Yogx) Qp9x,t-1 


Where: 
Onexs represents the quantity of U.S. natural gas exports at time f, 
Agx1 is a constant calibrated to the AEO market projections, 
Pz, is the price of natural gas at time f, 
Nex is the long-run price elasticity of export demand for U.S. natural gas, and 
YDex is the rate at which export demand for natural gas adjusts. 


U.S. natural gas exports are dependent only upon the domestic price of natural gas and no other prices, domestic 
or international. 


5.3 U.S. Natural Gas Supply 


Ngu 
Qsgut Fe Bout i ie ae (abe Ysqu)Qsgut-1 


for each domestic or imported natural gas source u = lower 48 conventional, lower 48 unconventional, lower 48 
offshore, Alaska onshore, Alaska offshore, other (e.g., synthetic natural gas and coke oven gas), pipeline imports, 
and LNG imports, where: 


Oseu,: tepresents the quantity of natural gas supplied to the U.S. market from domestic or imported source 
u at time f, 

But is a constant calibrated to the AEO market projections, 

Pz, is the price of natural gas at time f, 

Neu is the long-run elasticity of natural gas supply to the U.S. market from source u, and 

Ysgu is the rate at which natural gas from source u adjusts. 


6. Coal Market 


EnergySub represents the U.S. coal market as 14 separate sub-markets defined according to the region where coal 
is produced, with exports. The model also includes imports as exogenous to the model. The 14 coal markets in 
EnergySub correspond to the coal supply regions represented in the Coal Market Module of EIA’s NEMS, shown 
below in Figure E.2B.1. These supply regions are modeled separately to account for differences in the sulfur 
content, thermal value, rank, and production method of different coals. Because coal characteristics often differ by 
region (e.g, the Southern Powder River Basin region produces only low-sulfur, surface mined subbituminous 
coal), this approach (in most cases) implicitly captures the important differences between domestic sources of 
coal. With 14 distinct coal markets (one for each supply region), EnergySub estimates 14 equilibrium coal prices 
for each year. 
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Coal use for electricity generation is represented elsewhere in the electricity section of the model. The equations 
that follow present the model’s estimation of U.S. coal demand, demand for U.S. coal exports, and U.S. coal 


supply. 
Figure E.2B.1. EnergySub Coal Supply Regions 
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6.1 U.S. Coal Demand 


= Nci | | Neji 
Oped = ene 7 Verte f Pe Pr Gl e: Yoci) Qpcir,t-1 
J 


for each U.S. end-use sector i, for each coal supply region r; and j = g (gas), o (oil), and e (electricity) where: 


Opcir.: represents the quantity of coal demanded in sector i from coal supply region r at time f, 

Acir.: iS a constant calibrated to the AEO market projections, 

P.;1 is the minemouth price of coal from supply region r at time ¢, 

Nci is the long-run price elasticity of coal demand in sector 7, 

P;, is the price of energy source / at time f, 

nai is the long-run elasticity of demand for coal with respect to the price of energy source j in sector /, and 
Ybci is the rate at which demand for coal in sector 7 adjusts. 


Other than the electricity sector, whose coal demand is modeled separately, EnergySub’s domestic demand 
sectors for coal include industrial and other. 


6.2 Demand for U.S. Coal Exports 
Qnerx,t aa Acrx,t ‘ Be Hate Yocx) Qperx,t-1 


for each coal supply region, 7, where: 


Operx: represents the quantity of U.S. coal exports from coal supply region r at time /, 
Aor; is a constant calibrated to the AEO market projections, 
P.,. is the minemouth price of coal from supply region r at time f, 
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Ncx is the long-run price elasticity of export demand for U.S. coal, and 
YDex 18 the rate at which export demand for coal adjusts. 


Coal exports in EnergySub are only dependent upon the domestic minemouth price of coal from each coal supply 
region. No other energy prices, domestic or international, affect exports of coal. 


6.3 U.S. Coal Supply 


mt p | 
Qsert a Bye Pers rd 


for each coal supply region, r, where: 


Oscr,x represents the quantity of coal supplied to the U.S. market from coal supply region r at time f, 
B.r1 is a constant calibrated to the AEO market projections, 

Pert 18 the minemouth price of coal for coal supply region r at time ¢, 

Ner is the long-run elasticity of coal supply to the U.S. market from coal supply region r, and 

Ysc 1S the rate at which coal supply adjusts. 


As noted above, EnergySub treats coal imports as exogenous. For each BLM development scenario, imports are 
assumed to be the same as under the baseline scenario. The model makes this simplifying assumption because 
imports are projected to make up a de minimis fraction (less than 1 percent) of U.S. coal demand according to the 
AEO and imports do not align with the 14 coal markets specified in the model. 


7. Electricity Market 


Equations in EnergySub represents the U.S. electricity market and models U.S. exports and imports of electricity 
as net imports. The electricity sector in EnergySub also provides additional demand for oil, natural gas, and coal. 
The equations below present EnergySub’s approach for estimating U.S. electricity demand, U.S. electricity 
supply, and demand for fossil fuels for electricity production. 


To depict the use of coal for electricity generation with greater spatial detail, EnergySub divides the electricity 
supply market into nine regions (the U.S. Census Divisions), shown in Figure E.2B.2 below. Each electricity 
supply region is also modeled to receive coal from the 14 separate coal supply regions described above, resulting 
in a total of 126 total coal supply-electricity supply region combinations. 
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Figure E.2B.2. EnergySub Electricity Supply Regions 
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for each U.S. electricity end-use sector i; and j = g (gas), c (coal), and o (oil), where: 


Qpei: represents the quantity of electricity demanded in sector i at time f, 

Avi; is a constant calibrated to the AEO market projections, 

P.,is the price of electricity at time ¢, 

Nei is the long-run price elasticity of electricity demand in sector i, 

P;, is the price of energy source / at time ¢, 

Nei is the long-run elasticity of demand for electricity with respect to the price of energy source j in sector 
i, and 

Ybei is the rate at which demand for electricity in sector 7 adjusts. 


The U.S. demand sectors for electricity in EnergySub include (1) residential, (2) commercial, (3) industrial, (4) 
transport, and (5) other. As in the oil and gas markets, EnergySub uses a single weighted average minemouth 
price of coal instead of separate regional coal prices to estimate cross-price effects in the industrial and other 
sectors. 


7.2 U.S. Electricity Supply 


EnergySub uses separate approaches for the éstimation of electricity derived from gas and oil, coal, and electricity 
derived from other sources. While the quantity of electricity generated from gas, oil, and coal is dependent on 
fossil fuel prices, changes in these prices do not directly factor into the generation of electricity from non-fossil 
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energy sources.’ In addition, EnergySub accounts for the cost of transporting coal from each coal supply region to 
each electricity supply region by adding the coal transportation cost to the minemouth price of coal, which yields 
an estimate of the delivered price of coal. To account for this difference in the economics of electricity generation 
for different types of power producers, EnergySub specifies electricity supply separately for three classes of 
generation as follows: 


Osere — Cres(Per/Pi) i Ysej)Qsejt-1 
for 7 = oil and natural gas, where: 


Osej,, represents the quantity of electricity supplied from fossil fuel energy source / at time f, 
Cj, 1s a constant calibrated to the AEO market projections, 

P.; is the price of electricity at time ¢, 

P;1 is the price of fossil fuel energy source j at time 1, 

ne 1s the long-run elasticity of electricity supply from fuel 7, and 

Y sei 1s the rate at which electric power from fossil energy 7 adjusts. 


Oc 33 Corzt ‘ Pet! Pens lee eS + (1 ar Vsec)Osecrzt=1 


for c = coal, for each coal supply region r and each electricity supply region z, where: 


Osecrz,: represents the quantity of electricity supplied from coal supply region r to electricity supply region 
zat times, 

Cerz,: 18 a constant calibrated to the AEO market projections, 

P. 1s the price of electricity at time ¢, 

Per, 18 the minemouth price of coal from supply region r at time ¢, 

Tcrz tepresents the transportation cost of coal from coal supply region r to electricity supply region z, 

Nec iS the long-run elasticity of electricity supply from coal, and 

Y sec 18 the rate at which electric power from coal adjusts. 


As noted above, EnergySub accounts for the cost of transporting coal between each of the 14 coal supply regions 
and each of the nine electricity supply regions. The model therefore includes estimates of the per-ton cost of 
transporting coal (T;,z ) for all 126 combinations of coal supply and electricity supply regions. 


Qseie = Cit Pee"! + (1 — se) Qseit—1 
for /= nuclear, hydro, wind, solar, other electric, net imports, where: 
x p 


Qsei: represents the quantity of electricity supplied from source / at time ¢, 
Ci; is a constant calibrated to the AEO market projections, 


’ All else equal, renewable electricity generation in EnergySub will increase as fossil fuel prices rise, but the effect is indirect. 
For a given level of electricity demand, fossil fuel-based generators will supply less electricity as fossil fuel prices rise, which 
will shift generation toward renewables. 
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P. 1s the price of electricity at time f, 
Nei is the long-run elasticity of electricity supply from source /, and 
Yse1 is the rate at which electric power from source / adjusts. 


7.3 Demand for Fossil Fuels to Produce Electricity 


7.3.1 Oil and Natural Gas 


Qpje,t a Kj. : Qsejt 


for 7 = oil and natural gas, where: 


Onpje1 represents the quantity of energy source j used to produce electricity at time ¢, 
Kj, is a constant calibrated to the AEO market projections, and 
Osei: represents the quantity of electricity supplied from source / at time f 


7.3.2 Coal 


Oneeret =a Kear Cveaat 
Zz Zz 


for c= coal, where: 


dz Qpcerz,t is the sum of demand for coal from coal supply region r for electricity production across all z 
electricity production regions at time f, 


Keyt is a constant calibrated to the AEO market projections, and 


Dz Qsecrz,t is the sum of coal supplied for electricity production from coal supply region r across all z 
electricity production regions at time f. 


8. Model Calibration 


EnergySub calibrates the supply and demand equations outlined above to market conditions reflected in baseline 
long-run projections through the parameters 4, B, C, and K. These parameters are derived from the elasticities, 
adjustment parameters, market quantities, and prices in the long-run projections of energy production and 
consumption. They serve as constants in the model’s supply and demand equations and benchmark the model’s 
simulated market responses to observable market conditions in the initial baseline equilibrium. 


EnergySub has extensive data requirements and needs detailed long-run forecasts for the supply, demand, and 
prices of energy and individual energy sources to establish its initial market equilibrium and derive these the 
calibration parameters. EnergySub relies on projections developed by EIA for its baseline because annual 
projections developed using their integrated models provide the most complete data set for long-run energy 
market conditions. Data and statistics produced by EIA are widely accepted as best available information and 
regularly used by Members of Congress, industry participants, government agencies, and other interested parties. 


9. Equilibrium 
The equilibration calculation of EnergySub selects Po, Pgt, Per, and P.r, for each period t such that the quantity 


of oil, natural gas, coal (by coal supply region), and electricity supplied equals the quantity demanded in each 
period t. For coal, the national market not only needs to be in equilibrium but the quantity of coal supplied by 
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each coal supply region r at period ¢ must equal the quantity of coal demanded from coal supply region r at each 
period t. The model specifies these equilibrium conditions as follows: 


World Oil Market 


Qnoe.t + Qpoxt + Day Qpoi,t = Qsoy,t ape pr Osout 


where: 
Qpoe,t 18 the U.S. demand for oil to produce electricity at time ¢, 
Qpox,t 18 foreign demand for oil at time f, 
di Qpoi,t is the U.S. demand for oil across all other end use sectors i at time f, 
Qsoy,t 1s the oil supply from foreign sources at time ¢, and 


Yu Qsowt is the domestic oil supply from all domestic sources at time t. 


U.S. Natural Gas Market (with exports and imports 


Qnget + Dy Qpgit + Qogxt = > Qsgut 
i u 


where: 
Qpge,t is the U.S. demand for natural gas to produce electricity at time ¢, 
DAO Dgit is U.S. demand for natural gas across all end use sectors i at time ¢, 
Qpgx,t 18 the demand for U.S. natural gas exports at time f, and 


Dane Sout is the supply of natural gas from all u domestic sources at time tf. 


U.S. Coal Markets, by Supply Region 


> Once a » Oncine oF Opecnt ne Qsert 
Z i 


where: 


yz Qpcerzt is the quantity of coal demanded from coal supply region r across all electricity production 
regions z at time ¢, 


yi Qpcir,t is the quantity of coal demanded from each coal supply region r across all end-use sectors i at time 
t, 


Qpcxr,t 18 the quantity of coal demanded for exports from each coal supply region r at time ¢, and 


Qscr,t 1S the quantity of coal supplied by each coal supply region r at time f. 


US. Electricity Market (with net imports) 
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yy Qneit = » Qsejt as ». ys Osecrzt a yy Qseit 
I iI Zee l 


where: 
Yi Qpeit is the demand for electricity across all end-use sectors 7 at time /, 


yj Qse;,t 18 the supply of fossil fuel electricity (excluding coal), for all other j fossil fuel sources at time ¢, 


Yr dz Qsecrz,t is the supply of coal-fired electricity across all r x z electricity production regions at time ¢, 
and 


X12 Qsext is the supply of renewable electricity across all / renewable sources at time t. 


The equilibration process is initiated once a Reasonable Foreseeable Development Scenarios (RFD) for onshore 
federal oil, natural gas, and coal is introduced into the model. The RFD serves as a supply shock, moving the 
system of equations into a state of disequilibrium. These supply shocks can reflect an increase or decrease in the 
future supply of the corresponding energy source depending on whether production under the RFD is incremental 
to or a component of projected baseline supply. Once EnergySub’s system of equations are moved out of 
equilibrium, the model uses reduced gradient methods to solve its system of equation for new equilibrating prices 
for energy and energy sources. Solving for these new prices yields equilibrium supply and demand quantities of 
energy and energy sources, accounting for substitution between energy fuel sources. When zero disparity between 
supply and demand across all 17 fuel markets is achieved, EnergySub saves the market-clearing prices and 
proceeds to the next year in the production scenario to perform the same equilibration. 


10. Adjustment Rates and Elasticities 


All elasticities and adjustment rates in EnergySub have default values that were obtained from the literature, 
derived from NEMS supply curves, inferred from NEMS output, or obtained from BOEM’s MarketSim Model.® 
The sections below document the default adjustment rates and elasticities used in EnergySub when modeling 
production scenarios for the Willow Master Project. 


To the extent possible, EnergySub relies upon values from peer-reviewed studies in the empirical economics 
literature. Reliance on peer-reviewed data is central to ensuring that EnergySub’s simulated market responses 
reflect the best information available. In the few cases where peer-reviewed values are not available, elasticity 
estimates were derived from NEMS outputs or from expert input. 


10.1 Adjustment Rates 


EnergySub includes a series of adjustment rates in the supply and demand equations to capture the transition from 
short-run to long-run market effects. These adjustment rates account for the portion of demand or supply that is 
allowed to change from one year to the next. No data on the adjustment rates for specific energy sources are 
readily available. In the absence of such data, EnergySub’s default is to assume that adjustment rates are related to 
the retirement of energy producing and consuming capital (i.e., equipment that produces energy or consumes 
energy), as indicated by their average lifespan. When Reasonably Foreseeable Production scenarios include large 
changes in production volumes, relative to projected supplies in the previous year, adjustment rates may need to 


8 Many of the elasticities used from the BOEM MarketSim model were provided by energy economist Dr. Stephen Brown 
(2011) of the University of Nevada, Las Vegas (UNLV). See Industrial Economics, Inc. (2017). 
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be disabled to enable the model to equilibrate properly. For the modeling of Willow, adjustment rates were set to 
1 and allowed to drop out of supply and demand equations, enabling the model to accommodate the large year 
over year change in AK onshore oil production during the Project’s first year of production. 


10.2 Demand Elasticities 


EnergySub’s demand elasticities measure changes in the consumption of energy and energy sources relative to a 
percent change in price. EnergySub utilizes own-price and cross-price demand elasticities for each energy source 
included in the model to capture the complex interactions between different segments of U.S. energy markets. For 
each major energy consuming sector (e.g., the residential sector), BLM prioritized using own-price and cross- 
price demand elasticities from the same empirical study to ensure that each sector’s simulated responses were 
based on price sensitivities derived using the same methods, assumptions, and data. The selection of demand 
elasticities also considered the quality of the estimates produced by each study. BLM’s assessment of quality for 
individual elasticity estimates considered, among other factors, (1) whether they are statistically significant, (2) 
methods by which they were derived, and (3) the richness of the data supporting each estimate (e.g., whether they 
are based on a multi-year panel or reflect energy market data for a single year). 


Based on these criteria, EnergySub relies heavily on own-price and cross-price demand elasticities from Serletis et 
al. (2010) for the residential and commercial sectors and Jones (2014) for the industrial sector. Serletis et al. 
(2010) investigate inter-fuel substitution possibilities for energy demand across four fuels (i.e., oil, gas, electricity, 
and coal) using EJA data for the 1960-2007 period. Based on these data, Serletis ef al. estimated own-price and 
cross-price elasticities for the commercial, residential, and industrial sectors, using a flexible translog functional 
form. Across most sectors, Serletis et al. produced statistically significant elasticity values of the expected sign. 


Jones (2014) focuses on inter-fuel substitution in the industrial sector, using EIA data for the 1960-2011 period 
for the same fuels included in Serletis et al. (2010) plus biomass. Jones specifies a dynamic linear logit model to 
estimate own-price and cross-price elasticities, and within this framework, estimates both short-run and long-run 
elasticities. In addition, to assess the role of biomass in industrial sector inter-fuel substitution, Jones develops two 
sets of models, one including the four energy sources traditionally included in industrial sector energy models 
(i.e., natural gas, oil, coal, and electricity) and another that includes these energy sources plus biomass. Jones 
finds that the addition of biomass reduces both the own-price and cross-price elasticities of demand for the four 
traditionally modeled fuels. The effect is most significant for those values associated with electricity. In both 
models, the four traditional energy sources are found to be substitutes with each other with the exception of 
electricity and oil; the cross-price elasticities for these energy sources are not Statistically significant. 


Table E.2B.1 presents the default own-price and cross-price demand elasticities used in EnergySub for the 
residential, commercial, industrial, and transport sectors. The table also shows the default elasticity values for 
miscellaneous demand sectors included in EnergySub (e.g., natural gas demand in U.S. export markets). As 
indicated in the table, EnergySub uses results from Serletis et al. (2010) as defaults for the commercial and 
residential sectors, except for the elasticity of demand for natural gas with respect to the price of oil and the 
elasticity of demand for oil with respect to the price of natural gas. The estimates for these cross-price elasticities 
in Serletis e¢ al. were of the unexpected sign (negative) and were not statistically significant. Therefore, in lieu of 
Serletis et al., EnergySub uses results from Newell and Pizer (2008) for these values, for both the commercial and 
residential sectors. Newell and Pizer (2008) estimate these cross-price relationships for the commercial sector 
only. While EnergySub would ideally use default values specific to the residential sector, alternative values for 
these cross-price elasticities were not readily available for the residential sector. Given the similarities between 
the commercial and residential sectors, EnergySub uses these two cross-price demand elasticities from Newell 
and Pizer (2008) as a reasonable approximation of the corresponding residential sector values. 


For the industrial sector, EnergySub relies almost exclusively on demand elasticities from Jones (2014) as 
defaults. Although Serletis et al. (2010) estimate elasticity values for the industrial sector, the values in Jones 
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(2014) are based on fuel consumption data that exclude fuel use for purposes other than energy (e.g., petroleum 
products used as lubricants). As described above, Jones (2014) estimates long-run demand elasticities with two 
specifications, one including biomass as a substitute and another excluding biomass. Based on the statistical 
significance of the elasticities with biomass included, EnergySub uses the elasticities from the specification that 
includes biomass. The two exceptions to this are the cross-price elasticity of demand for oil with respect to the 
price of electricity and the cross-price elasticity of electricity in response to oil prices, as Jones’ estimates for 
these values are not statistically significant. For these values, EnergySub uses estimates from Serletis ef al. 
(2010). 


Table E.2B.1 also shows EnergySub’s default own-price demand elasticities for the transport sector and various 
miscellaneous demand categories. For these categories, EnergySub relies upon elasticity values from multiple 
sources. For oil demand in the transportation sector, EnergySub uses a U.S.-specific elasticity value obtained from 
Dahl’s (2012) review of price elasticities estimated for more than 100 countries. This value represents the average 
of the elasticity values identified in the empirical literature. For non-U.S. oil demand, EnergySub applies the 
value reported in a Huntington et al. (2019) review of crude oil demand elasticities in major industrializing 
economies. For U.S. natural gas exports, EnergySub uses estimates from Dahl’s prior (2010) review of the 
elasticity literature as defaults. 


Two categories for which appropriate demand elasticity values were not identified in the literature are 
miscellaneous coal demand and demand for U.S. coal exports. EnergySub uses the same industrial sector value 
obtained from Jones (2014) for the former and assumes a value of -1.00 for the latter. 
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Table E.2B.1. Energ 
ENERGY SOURCE 


Sub Default Demand Elasticities 
ELASTICITY ELASTICITY 


ELASTICITY WITH ELASTICITY 
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-0.264 
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0.009 
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Transport Sector? -0.300 
Oil — Rest of World -0.15 
Demand for US Crude* 
Oil — Rest of World 
Demand for US Refined 
Products* 
Oil — Rest of World -0.15 
Demand for non-US oil4 


Natural Gas — Transport> 


Natural Gas — US Export 
Markets® 

Electricity — Transport? 
Electricity — “Other” 


Coal — Other® 
Coal — US Export Markets> 


Notes: 

1. Commercial and residential sector values are from Serletis et al. (2010), except for the cross-price elasticity for gas in 
response to oil prices and the cross-price elasticity of oil in response to gas prices. For these latter two values, 
EnergySub uses demand elasticities from Newell and Pizer (2008). Also, Deryugina et al. (2020) estimate a range of 
residential elasticity values for electricity consistent with the value in Serletis et al. (2010). 

2. For the industrial sector, EnergySub uses demand elasticities from Jones (2014), except for the cross-price elasticity of 

electricity in response to oil prices and the cross-price elasticity of oil in response to electricity prices. For these values, 

EnergySub uses demand elasticities from Serletis er al. (2010). 

. Dahl (2012) 

. Huntington et al. (2019) 
. Assumed to be -1.00. 

. Dahl (2010) 


7. Assumed to be average of own-price elasticity values for industrial, commercial, and residential sectors 
Industrial sector value from Jones (2014). 


-0.15 


-1.00 
-0.89 


-1.00 
-0.18 


-1.468 
-1.00 


Nn BW 
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10.3 Supply Elasticities 


EnergySub includes default supply elasticities, summarized in Table E.2B.2, for every production category 
modeled for a given fuel (e.g., onshore tight oil production in the lower 48 states). These supply elasticities 
measure how responsive energy producers are to changes in market prices. Consistent with the demand elasticities 
summarized above, several of EnergySub’s supply elasticities were obtained from the economic literature, with 
data sources varying by fuel type. 


For tight oil and other lower 48 onshore oil, EnergySub uses elasticities from a recent study by Newell and Prest 
(2019). The paper specifically compares the price responses of conventional and unconventional (tight) oil 
drilling and production. Using micro-data for more than 150,000 oil wells in Texas, North Dakota, California, 
Oklahoma, and Colorado, Newell and Prest (2019) estimate the elasticity of well drilling and the elasticity of oil 
production, separately for conventional and unconventional wells. To estimate drilling elasticities, they use 
multiple model specifications, estimating changes in drilling activity as a function of price in some cases and as a 
function of revenue in other cases. The production elasticities estimated by Newell and Prest (2019), however, all 
represent the change in production as a function of the change in revenue, rather than price. To align the supply 
elasticities in EnergySub with the specification of supply, EnergySub uses the elasticity of well drilling with 
respect to the oil price from Newell and Prest (2019), which they estimate separately for both conventional and 
unconventional wells. 


Luchansky and Monks (2009) serves as the source for EnergySub’s default supply elasticity for domestic 
biodiesel. This paper uses monthly data for 1997 through 2006 to estimate the market supply and demand for 
ethanol at the national level. Applying these data to four specifications of supply, Luchansky and Monks (2009) 
estimated supply elasticities ranging from 0.224 to 0.258. EnergySub uses the midpoint of this range (0.24) as the 
default supply elasticity for biodiesel. 


For a number of oil supply elasticities, EnergySub relies on values included in BOEM’s MarketSim model based 
on expert input provided to BOEM by three energy economists: Dr. Charles Mason of the University of 
Wyoming, Dr. Seth Blumsack of Penn State University, and Dr. Gavin Roberts of Weber State University. 
EnergySub relies on input provided to BOEM by these experts for the oil supply elasticities related to lower 48 
offshore, rest-of-world oil production, Canadian pipeline imports, natural gas plant liquids, and other oil 
production. For oil production in Alaska, EnergySub uses supply elasticities derived from specialized simulations 
of NEMS, as described in detail below. 


For gas production, EnergySub draws on a variety of sources for elasticities, depending on the production source. 
For domestic onshore conventional and unconventional shale gas production in the lower 48, EnergySub uses 
values from Newell, Prest & Vissing (2019), who use data from approximately 62,000 gas wells drilled in Texas 
between 2000-2015 to determine price-responsiveness across the supply process. The study assesses the decision 
to drill the well, well completion, and produce gas over time and, of these, finds drilling activity to be the most 
responsive to changes in price. EnergySub makes use of the gas price response values broken out for conventional 
and unconventional wells, though the study notes that these values may not differ significantly from each other 
statistically. For offshore production in the lower 48, EnergySub uses the same 0.19 elasticity as for offshore oil 
production in the lower 48, obtained through the expert input process described above. For onshore and offshore 
production in Alaska, EnergySub uses elasticity values derived from specialized simulations of NEMS, as 
detailed below. For other gas production, EnergySub applies the supply elasticity reported in Brown (1998). 


Table E.2B.2. EnergySub Default Supply Elasticities 


SOURCE/ SUPPLY ELASTICITY 


Lower 48 Onshore Non-Tight! 0.93 Other? 
Lower 48 Onshore Tight! 0.73 Biodiesel* 0.24 
Lower 48 Offshore? 0.19 Rest of World? 0.28 
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FUEL SOURCE/ SUPPLY ELASTICITY 


Alaska Onshore? 0.42 Nae Plant 
Liquids 

Natural 

Gas 


Canadian Pipeline 


Alaska Offshore? 0.58 0.38 


Imports” 
Alaska Offshore? 
Other’ 0.51 
Pipeline Imports® 0.52 
LNG Tanker Imports? 
Hydro? 


Lower 48 Conventional° 
Lower 48 Unconventional* 0.68 
Lower 48 Offshore® 0.19 
Alaska Onshore? 
Oil! 


Natural Gas? 1.50 Wind Onshore? 0.65 
Coal!° 0.27 Wind Offshore? 0.01 
Nuclear? 0.53 Solar? 2.03 


Other Electric? Imports? 


Northern Appalachia"! WY PRB — North!° 
Central Appalachia!! 4.62 WY PRB - South!! Sala 
Southern Appalachia"! 1.50 Western Wyoming"! 0.73 
East Interior’! 7.40 Rocky Mountain!! 2.43 
Coal West Interior"! Arizona/New 

0.47 Meica BAS 
Gulf Lignite”™! Be? Alaska/Washington'" — 0.60 
Dakota Lignite!! Imports? 1.00 
Western Montana!! 


Notes: 
1. Newell and Prest (2019). 
. Expert input from C. Mason, G. Roberts, & S. Blumsack, as cited in Industrial Economics Inc. (2021). 
. Derived from AEO (2020). 
. Luchansky and Monks (2009). 
. Newell, Prest & Vissing (2019) 
. Assumed to be the same.as Oil, Lower 48 Offshore 
. Brown (1998). 
. Derived from specialized NEMS run of the AEO 2015 provided to DOI by EIA. 
9. Assumed value. 
10. Derived from AEO 2018a, as provided by BOEM (2018). 
11. Derived from NEMS 2019 Reference Case supplemental data provided to BLM by EIA. 


COND NW BW LO 


For coal supply, EnergySub uses supply elasticities unique to each of the 14 coal supply regions, as derived from 
annual supply curve data generated by NEMS’ Coal Market Module (CMM).” The annual supply curve data 
provided by EIA represent 41 distinct coals for a given year for combinations of coal supply region, sulfur 
content, mining method, and rank. For example, the Central Appalachia coal supply region has five different 
supply curves for a given year, representing a mix of low- and medium-sulfur coal, underground and surface 
mines, and premium and bituminous coals. In addition, the annual supply curve for each of the 41 coals is 
represented as 11 data points, with each data point representing production at a given price point. 


Using the EJA data, we estimated supply elasticities for each of the 41 coal types, for every year between 2019 
and 2040. To generate elasticity values, we applied the standard econometric method of regressing the log- 
transformed price on the log-transformed quantity, which yielded the elasticity of supply as the coefficient. Each 
regression was performed over the three central points of the appropriate supply curve. The following equation 
displays this regression: 


In(Qs) a Bs|n (Pst) + Bo 


° While not publicly available, EIA provided these supply curve data for the purposes of this project and provides them to 
other modelers on a regular basis. 
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Q; + represents the quantity supplied on supply curve s in year f, 
Bt represents the elasticity of supply for supply curve s in year 7°, 
P, + represents the price of coal on supply curve s in year t, and 

Bo represents the regression constant. 


Running the above regression for each of the 41 supply curves for every year between 2019 and 2040 yields an 
initial set of elasticities. To convert the year-specific and supply curve-specific results to regional supply 
elasticities, we developed a weighted average coal supply elasticity for each of the 14 coal supply regions across 
all years, using the quantity associated with the coals produced by each coal supply region as weights. Table 
E.2B.2 above displays the results of the supply elasticity calculation for each coal supply region. 


Where appropriate economic research does not exist or could not be obtained for a specific supply elasticity value, 
projections from the AEO were used to infer these values.'' Elasticity estimates may be inferred from the AEO 
projection for a given year by comparing the differences in energy prices between two scenarios with the 
differences in energy quantities. For a given energy source and fuel, an annual inferred elasticity value was 
calculated three times: (1) based on the low oil price case vs. the high oil price case, (2) the low price case vs. the 
reference case, and (3) the reference case vs. the high price case, for all AEO projection years from 2017 through 
2040. The formula for this annual inferred elasticity is as follows. 


(oe) 


in (Bet) 
Pet 


Where 77; is the inferred elasticity in year ¢, Qa,;,and Oz-represent the quantities supplied in year ¢ for cases A and 
B respectively (each case is compared with both of the other cases), and P4,and Pz, are the prices at time ¢ for 
cases A and B. The resulting series of inferred elasticities are averaged, excluding extreme outlier results derived 
from the AEO data.'? 


For a limited number of producing sectors, elasticity values were unavailable from the literature and the data 
generated by the constrained NEMS run or recent editions of the AEO yielded elasticity values that appeared 
unrealistically high or were insufficient to support estimation of a supply elasticity. In such cases, EnergySub uses 
a default supply elasticity of 1.0. 


11. Limitations 


As described above, EnergySub uses a system of equations to create a mathematical representation of complex 
energy markets. Like all mathematical models that try to simplify real world phenomenon, EnergySub is limited 
by data constraints and simplifying assumptions. The BLM designed this model to simulate potential market 
responses to changes in the price of energy and energy sources stemming from marginal changes in the supply of 
onshore oil, natural gas, and coal, given long-run market conditions projected by the EIA. EnergySub is not a 
national forecasting model or a replacement for NEMS, WEPS, or the long-run energy projections developed by 
EIA’s Office of Energy Analysis using these integrated models. 


' Coal supply elasticities are also represented as 7, in Equation 1. 
'l Tn some cases, the supply elasticities were derived from prior releases of the AEO rather than AEO 2020 when results from 
the 2020 data resulted in unrealistic elasticity values. 


'2 More specifically, elasticities were estimated based on differentials between the low-price case and reference case, the 
reference case and the high-price case, and the low-price case and the high-price case. They then were averaged across these 
three variants and across years. 
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Energy markets are influenced by a number of factors that change over the long run, including consumer 
behavior, technological innovation, system constraints, and government regulations. However, these factors are 
not explicitly modeled in EnergySub. Instead, EnergySub relies on its calibration to EIA’s long-run market 
projections to capture market constraints and patterns of energy production and consumption implicitly through 
the baseline equilibrium and calibration parameters. Since the majority of the model’s assumptions about the 
underlying structure of energy markets are adopted implicitly, EnergySub users have limited ability to fine tune 
individual assumptions about factors that can reshape energy markets. This makes EnergySub a suitable model for 
simulating market responses to supply shocks, but a less useful tool for developing independent energy forecasts 
that can serve as benchmarks for other simulations. 


The EnergySub model cannot use multiple sets of elasticities and parameters within a single simulation. This 
causes simulated market responses to remain static each year of the simulation. Although real-world energy 
consumption patterns and market responses may change over the long run, elasticities and parameters within 
EnergySub must remain constant because simulations for long-run production scenarios cannot be broken into 
shorter segments with model specifications adjusted manually between different segments of time. 


EnergySub’s equations for oil and natural gas markets contain limited regional detail. Although EIA develops 
long-run projections for domestic oil and gas production based on detailed estimates of technically recoverable 
resources (regardless of mineral ownership) and play-level analysis, AEO results reflect aggregate production 
across mineral ownership and producing basins for a limited number of domestic supply sources (e.g., Alaska 
onshore, Alaska offshore, Lower 48 offshore, Lower 48 onshore conventional, Lower 48 onshore 
unconventional). While this level of detail is sufficient to estimate substitution effects across aggregated markets 
for these different supply sources, it provides insufficient regional detail to assess substitution effects across 
regional submarkets within Alaska or the Lower 48 or identify specific geographic shifts in production that could 
occur. In addition, because mineral ownership is not directly captured in the EIA data or in EnergySub’s model 
specifications, EnergySub cannot capture substitution among federal, state, and private minerals within energy 
markets or across regional submarkets. 


EnergySub currently models foreign oil supply and demand using a limited number of equations that represent 
non-U.S. sources for oil and non-U.S. demand for U.S. crude oil, non-U.S. demand for U.S. refined products, and 
non-U.S. demand for non-U.S. oil. Estimating impacts to foreign submarkets is currently beyond the modeling 
capabilities of EnergySub. While EnergySub is able to estimate how changes in global oil prices affect the supply 
and demand for U.S. and non-U.S. sources of oil, it is not able to determine which foreign countries may be 
increasing or decreasing their consumption and production of oil based on these changes in prices. 


BLM acknowledges there are limitations to EnergySub and not all potential market responses and impacts were 
quantified within the model. Nevertheless, results from EnergySub modeling still provide valuable insights to 
how mineral related management decisions affect energy market conditions. 


12. Application of EnergySub to Willow 


EnergySub was used as a tool to compare unobservable long-run market conditions with and without production 
from Willow to gain insights to how production from the Project may affect projected energy prices, production 
and consumption. As discussed above, EnergySub relies on long-run projections of energy production and 
consumption developed by the EIA to establish initial market equilibrium conditions. EnergySub’s system of 
equations was calibrated to projections from the 2022 Annual Energy Outlook for the modeling of Willow. This 
baseline largely reflected modeled equilibrium market conditions from the AEO Reference case, however, AEO 
data was supplemented with additional data on international supply and demand. The AEO Reference case 
reflects EIA’s best assessment of how U.S. and foreign energy markets will operate through 2050 based on key 
assumptions. EIA considers the Reference case to be a reasonable baseline from which to measure the impact of 
alternative scenarios and assumptions (EIA, 2022). The BLM considered using alternative energy projections 


’) 
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developed as AEO side cases but decided to refrain from introducing its own assumptions about long-term market 
conditions and calibrated the model to supply and demand projections from the Reference case.!° 


AEO projections for domestic oil and gas production are developed within NEM’s Oil and Gas Supply Module 
based on estimates of total technically recoverable resources within the U.S. AEO 2022 supply projections 
include production in Alaska between 2022 and 2050 from both existing fields, including expansion fields around 
the Prudhoe Bay and Alpine Fields for which companies have already announced development schedules, and 
undiscovered fields that most likely exist based on the region’s geology (EIA, 2022c). Since AEO projections 
reflect what market conditions may look like with production from Willow, EnergySub was used to simulate 
market responses to foregoing oil from the Project in order to create a counterfactual of what energy market 
conditions may look like in the absence of Willow. 


Results from the EnergySub simulations include estimates of the energy substitutes that would replace forgone oil 
production from Willow if the Project was not approved, and changes in foreign oil consumption (i.e., foreign 
demand for U.S. and non-U.S. produced oil) stemming from changes in onshore oil production from Willow. 
Conversely, these energy substitutions can be interpreted as the alternative energy sources displaced by oil 
produced from Willow. 


12.1 Energy Substitute Effects 


EnergySub estimates displaced energy substitutes by converting all consumption of energy into barrel of oil 
equivalents to enable the comparison of energy consumption across fuel sources and uses. Estimates of energy 
and energy sources potentially displaced by oil produced under Willow’s alternatives are reported below in Table 
E.2B.3. Simulated substitution effects for Alternatives B and C are identical because the timing and level of 
production anticipated under these alternatives is the same. Although production scenarios under Alternatives D 
and E differ slightly in volumes and timing from those under Alternatives B and C, market substitution effects 
across the four action alternatives are similar. EnergySub’s market simulations showed that oil produced from 
Willow would primarily displace oil produced from other domestic locations (i.e., AK offshore and lower 48 
onshore and offshore) or foreign sources under all action alternatives. Approximately 30% of oil produced from 
Willow was simulated to displace oil produced elsewhere domestically, while 52% of Willow’s production under 
the Alternatives B, C, D, and E displaced oil that would have otherwise been imported via tankers or pipelines 
from foreign producers. 


Oil is a global commodity, with prices determined by global supply and demand. Increased Alaska oil production 
from Willow was simulated to lower global oil price. At peak production, production under the four action 
alternatives was simulated to reduce the price of a barrel of oil by about 20¢. Lower oil prices relative to prices for 
substitute energy sources is likely to increase demand for oil and be a catalyst for fuel switching from alternative 
fuels sources. EnergySub simulations estimated that between 9 and 9.5% of production under the action 
alternatives displaced energy consumed by the residential, commercial, transportation, and electricity sectors from 
other fuel sources. As shown in Table E.2B.3, between 7.4 and 7.7% of total production under the action 
alternatives was simulated to displace biofuels, natural gas liquids, and electricity generated from nuclear and 
renewable energy sources. Of these displaced energy substitutes, 0.5% of the oil from the Project was simulated 
to displace electricity that would have otherwise been generated without fossil fuels. Simulations also showed that 
1.4% of Willow’s production under the action alternatives would displace consumption of energy from natural 
gas. The remaining 0.4% of production under the four alternative was simulated to displace coal. 


In addition to displacing energy from substitute energy sources, oil production over the life of Willow was shown 
to increase overall energy consumption during the life of the Project. While reductions in the consumption of 

electricity and energy from coal and natural gas would be the energy equivalency to 2% of oil produced under the 
action alternatives, nearly 11% of total oil production under the alternatives would represent new oil consumption 


'3 When developing annual energy projections, EIA runs side cases to show how NEMS and WEPS responds to changes in 
key input variables compared with the modeled results from the Reference case. 
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that would be unrealized at higher prices without production from the Project. The net effect of Willow’s 
production on overall energy consumption would be positive under the four action alternatives. Net increases 
would be slightly lower under Alternative D relative to those under Alternatives B, C, and E because of 
differences in how much and when oil from the Project was projected to enter the energy market. 


Table E.2B.3: Displaced Fuels and Changes in Demand (Substitution Effects 2023 through 2052 


PERCENT OF WILLOW MASTER OIL ALT B ALT C ALT D ALT E 
THAT: . 


Displaces Domestic Oil 30.1% 30.1% 30.1% 
Displaces Oil Imports 52.0% 52.0% 
Displaces Natural Gas 1.4% 1.4% 
Displaces Coal 0.4% 0.4% 
Displaces Biofuels and NGL 7.2%, 6.9% 
Displaces Electricity from Renewable Sources 0.5% 0.5% 
Changes in Demand* Alt B Alt E 
Oil 10.8% 10.8% 
Natural Gas -0.8% -0.9% -0.8% 
Coal -0.3% -0.3% -0.3% 
Electricity -0.9% -0.9% -0.9% 


*Change in demand does not represent displaced energy or energy sources. 


12.2 Changes in Foreign Oil Consumption 


As outlined in Section 4 (Oil Market), EnergySub models a single foreign oil market using a limited number of 
supply and demand equations. Foreign oil consumption is equal to non-U.S. demand for U.S. crude oil, non-U.S. 
demand for U.S. refined products, and non-U.S. demand for non-U.S. oil. Changes in foreign oil consumption are 
strictly a function of global oil prices, where demand for both U.S. and non-U.S. oil increases as global oil prices 
decrease. As discussed in Section 11 (Limitations), assessing impacts to foreign submarkets is beyond the current 
modeling capabilities of EnergySub. 


Production from Willow will increase the global supply of U.S. crude oil and U.S. refined products and bring 
down global oil prices. Simulations showed that peak production from Willow may decrease the global price of 
oil by approximately 20¢ per barrel. Since oil is a global commodity, lower oil prices is beneficial to domestic and 
foreign consumers and spurs additional demand for both U.S. and non-U.S. sources of oil. Relative to market 
conditions where oil from Willow is foregone, total foreign oil consumption was simulated to be slightly higher 
when global oil supplies included production from Willow. Ina global oil market where annual foreign demand is 
projected to exceed 38,000 million barrels (MMb) by the end of 2050 (IEO 2021), peak production under 
Alternatives B, C, D, and E was simulated to increase total foreign oil demand by approximately 13.7 MMb, more 
than 12 MMb of which constituted new foreign demand of oil from non-U.S. sources. Over the producing life of 
the Project, production from Willow was simulated to increase overall foreign consumption of oil by 124.4 MMb 
under Alternatives B and C, 130.4 MMb under Alternative D, and 120.1 MMb under Alternative E. 


12.3 Uncertainty 


The EnergySub results presented above reflect modeled market responses and energy substitutes likely to be 
displaced by oil from Willow. These energy substitutes also reflect the energy source most likely to replace 
foregone oil if the Project is not approved. They were derived using baseline projections of what energy markets 
will look like through 2050, elasticities which provide measures for how supply and demand between alternative 
energy sources may change in response to changes in prices, and production schedules provided by the proponent. 
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Results from the modeling are an estimation of what may happen in the future based on key assumptions about 
market conditions and production under alternative development scenarios. 


Energy markets are dynamic and projections about future market conditions are inherently uncertain because 
many of the events that shape patterns of energy consumption and production cannot be predicted with certainty. 
The baseline projections used for this modeling reflect EIA’s best assessment of how markets will operate 
through 2050 under a simplifying assumption that current regulations and consumption patterns will not change 
over the long term (EIA, 2022).'* BLM acknowledges that new laws and policies governing energy production, 
efficiency, and GHG emissions are likely to be enacted, and that these regulations may have significant 
implications for energy markets and substitutes in the coming decades. EIA will continue to incorporate new 
legislature and regulations into their modeling as funding and implementing regulations for them are enacted, and 
BLM will continue to evaluate the suitability of new data for future EnergySub calibration. 


The EnergySub modeling for Will does not account for structural changes that would have to occur within current 
energy markets to meet climate commitments and achieve net-zero emission goals. As the U.S. works towards 
achieving net-zero, energy production and consumption patterns will change. Energy markets may become 
increasingly electrified through greater deployment of renewable energy sources, enabling sectors that have 
historically been heavily reliant on fossil fuels to reduce their demand and consumption for these fuel sources. 
Technological innovation will also play a significant role in transforming how energy will be produced and 
consumed, though its implications for specific fuel sources and uses is not known at this time since many of the 
technologies have yet to be developed or economically scaled for widespread adoption. 


Even in a low carbon future, fossil fuels are likely to continue to play a role in the U.S.’s energy portfolio. 
Princeton’s Net-Zero America Project has been developing pathways to achieve net-zero emissions by 2050 using 
existing technologies. Four of their five pathways projected that oil and gas consumption would continue beyond 
2050, and that carbon capture and sequestration technology would play an important role in offsetting emissions. 
Under their fifth scenario, oil and gas are phased out by 2050 but oil continues to account for more than 20% of 
the energy fuel mix until the late 2030’s (Larson et al. 2020). Researchers and industry experts are continuing to 
explore potential pathways for decarbonization and the role of fossil fuels and other energy sources in a low 
carbon economy is still uncertain. Specific data on how the energy transition will affect demand for fossil fuels 
and alternative energy sources is not yet available. 


BLM acknowledges that energy substitutes for Willow may look significantly different in a low carbon future, 
and that modeling substitution effects using data that depicts current energy markets and historical market 
responses produces results with inherent uncertainty. As the energy transition progresses, reliance on other supply 
sources of oil (including other domestic production in the Lower 48 and foreign imports) to replace energy 
associated with oil from Willow may wane over the producing life of the Project. However, the timing and degree 
to which domestic energy markets may become less reliant on these alternative sources of supply is highly 
uncertain. BLM will continue to evaluate and update its methods for estimating energy substitutes as new energy 
statistics and information becomes available. 


'4 The version of NEMS used by EIA to produce the AEO 2022 included current legislation and environmental regulations 
for which implementing regulations were available as of the end of November 2021. The potential effects of proposed or 
hypothetical federal and state legislation, regulations, and standards—or sections of legislation that have been enacted but 
lacked funds to execute or did not have the required implementing regulations in place as of the end of November 202 1— 
were not reflected in NEMS when the AEO 2022 projections were developed. Additional information on the assumptions 
underlying the AEO 2022 Reference case are available at www.eia.gov/outlooks/aeo/assumptions/ 
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LIST OF ACRONYMS 


AAAQS Alaska Ambient Air Quality Standards 

ADEC Alaska Department of Environmental Conservation 
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AQRVs air quality related values 

BLM Bureau of Land Management 

CAP criteria air pollutant 

CASTNET Clean Air Status and Trends Network 
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dv deciview . 

EPA Environmental Protection Agency 
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HAP hazardous air pollutant 
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1.0 AIR QUALITY 

The U.S. Environmental Protection Agency (EPA) has determined that 50 kilometers (km) (31 miles) is sufficient 
to determine whether an emissions source will cause or contribute to exceedances of ambient air quality standards 
and is the approved distance for regulatory near-field air quality models (40 CFR 51, Appendix W). The far-field 
(regional) modeling domain is more than 300 km (186 miles) from the Willow Master Development Plan Project 
(Project) in all directions except south of the Project, where the closest point is approximately 250 km (155 
miles). 


1.1 Affected Environment 


1.1.1 Regulatory Framework 
In Alaska, the Alaska Department of Environmental Conservation (ADEC) has the authority to implement and 


enforce the Alaska Air Quality Control Regulations (18 AAC 50) through an EPA-approved State Implementation 
Plan. The Alaska Ambient Air Quality Standards (AAAQS) were promulgated in 18 AAC 50.010. The National 
Ambient Air Quality Standards (NAAQS) and AAAQS are provided in Table E.3.1. 


Table E.3.1. National and Alaska Ambient Air Quality Standards 
Pollutant? Averaging NAAQS> NAAQS? AAAQS“ = Form 
Secondary 


12 ug/m? 
35 ug/m? 


Not to be exceeded more than once per year on 
average over three years 


3 hours 


Note: AAAQS (Alaska Ambient Air Quality Standards); CO (carbon monoxide); N/A (not applicable); NAAQS (National Ambient Air Quality Standards); 
NO) (nitrogen dioxide); O3 (ozone); PM2 5 (particulate matter less than 2.5 microns in aerodynamic diameter); PM yo (particulate matter less than or equal to 
10 microns in aerodynamic diameter); ppb (parts per billion); ppm (parts per million); SO, (sulfur dioxide); .g/m* (micrograms per cubic meter). 

* Lead and ammonia are not shown as they are not pollutants of concern in the analysis area. 

® Source: 40 CFR 50 

© Source: 18 AAC 50.010 

4 All AAAQS are primary except for 3-hour SO). 


EPA designates geographic areas demonstrating compliance with the NAAQS as “attainment,” while areas that 
exceed the NAAQS are designated as “nonattainment.” If there is insufficient data to designate an area as 
“attainment” or “nonattainment,” the area will be designated as “unclassifiable.” The analysis area for air quality 
is designated as “attainment/unclassifiable” for all criteria air pollutants (CAP). 


The closest Class I area to the Project is Denali National Park, which is located more than 700 km (435 miles) 
south of the Project and is not in the analysis area for air quality. The three assessment areas within the far-field 
analysis area for air quality are Gates of the Arctic National Park, Noatak National Preserve, and the Arctic 
National Wildlife Refuge (Figure E.3.1). These three areas were selected following input from Federal Land 
Managers during discussions in the initial stages of the Willow MDP EIS process. The Class II prevention of 
significant deterioration (PSD) increments are presented in Table E.3.2. 


The air quality related values (AQRVs) are resources that may be affected by a change in air quality (NPS 2011). 
The Federal Land Managers’ Air Quality Related Values Work Group identifies AQRVs as “visibility or a 
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specific scenic, cultural, physical, biological, ecological, or recreational resource identified by the FLM [federal 
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Figure E.3.1. Analysis Areas for Air Quality and Regional Ambient Air Quality Monitors, Three Federally 
Managed Assessment Areas, and the Far-Field (Regional) Modeling Domain 
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Table E.3.2. Prevention of Significant Deterioration Increments for Class II Areas 
Pollutant Averaging Time Class II PSD Increment Form 


(ug/m?) 
Annual mean, not to be exceeded 
Not to be exceeded more than once per year 
9 
2 : 
3 


Oz 
SO. | 24hours | 9 Ss Not to be exceeded more than once per year 

SO7- | Annual | 20S [Anal mean, not to be exceeded 

PMos_ | = 24hours | 9 S_—S Not to be exceeded moree than once per year 

Annual mean, not to be exceeded 
eee Ue re us 


Not to be exceeded more than once per year 
Annual mean, not to be exceeded 
Source: 40 CFR 52.21 


Note: NO; (nitrogen dioxide); PM> ; (particulate matter less than 2.5 microns in aerodynamic diameter); PM o (particulate matter less than or equal to 10 
microns in aerodynamic diameter); PSD (prevention of significant deterioration); SO (sulfur dioxide); ug/m* (micrograms per cubic meter). 


D 
l 
0 
0 
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Visibility is a measure of how far and well we can see into the distance and is sensitive to changes in air quality. 
Visibility impairment (i.e., haze) occurs when sunlight is absorbed or scattered by tiny particles (e.g., sulfates 
[SO4*], nitrates [NO3’]) and gases (e.g., nitrogen dioxide [NO>]) (EPA 2017). The absorption and scattering of 
light impairs visibility conditions (i.e., visual range, contrast, coloration). Haze causing pollutants can be directly 
emitted or formed through the reaction of precursor gases emitted into the atmosphere (e.g., formation of SOy 
from sulfur dioxide [SO2]). The Regional Haze Rule (RHR) was promulgated in 1999 to improve and protect 
visibility in Class I areas (40 CFR 51.308). The Project area is not a Class I area; however, the RHR can inform 
current conditions and assessment of progress in required visibility improvements. The RHR defines reasonable 
progress goals to improve visibility on the most impaired days and ensure no degradation on the least impaired 
days, with the goal of attaining natural conditions (i.e., estimated visibility conditions in the absence of human- 
made air pollution) in each Class I area by 2064. Under the RHR, visibility is quantified using the deciview (dv) 
haze index, which is derived from light extinction. An incremental change in dv corresponds to a uniform and 
incremental change in visual perception for the entire range of visibility conditions. Single-source impacts on 
visibility are assessed by comparing the 98th percentile of the source contribution to the haze index to defined 
thresholds. A source that exceeds 0.5 dv (approximate 5% change in light extinction) is considered to contribute 
to visibility impairment, while a source that exceeds 1.0 dv (approximate 10% change in light extinction) is 
considered to cause visibility impairment (FLAG 2010). 


Atmospheric deposition can negatively affect ecosystems and other AQRVs. Dry deposition is continuous while 
wet deposition can only occur in the presence of precipitation. Potential deposition impacts include, but are not 
limited to, acidification of soils and waterbodies and nutrient enrichment (FLAG 2010). Wet or dry deposition of 
acidic pollutants formed from emitted SO2 and nitrogen oxides (NOx) is referred to as acid rain (EPA 2018b). 
There are currently no federal standards for atmospheric deposition, but FLMs use critical loads and Deposition 
Analysis Thresholds for assessing both cumulative impacts and source-specific impacts from new or modified 
PSD sources. A critical load is the level of deposition below which no harmful effects to an ecosystem are 
expected. Deposition Analysis Thresholds are screening thresholds that define the additional amount of deposition 
within an FLM’s area below which impacts are considered negligible. 


The National Emission Standards for Hazardous Air Pollutants defines maximum achievable control technology 
(MACT) standards that are technology-based standards for each regulated source category. MACT is applicable 
to all major sources (potential to emit more than 10 tons per year of a single hazardous air pollutant [HAP] or 25 
tons per year of any combination of HAPs) and some area sources (any stationary source of HAPs not classified 
as a major source) in specific source categories. 


1.1.1.1 Flaring Regulations 

Flaring in Alaska is regulated by three agencies: the Alaska Oil and Gas Conservation Commission (AOGCC), 
the ADEC, and the Bureau of Land Management (BLM). Flares are important safety devices that are used to 
ensure controlled combustion of natural gas to avoid a potentially explosive environment if the gas were to be 
vented to the atmosphere rather than flared. Flares would be used for gas released to prevent over pressurizing 
piping and equipment, to handle gas removal from systems during maintenance, and to address gas released 
during an emergency rapid depressurization of Willow Processing Facility gas handling systems (SLR 2022). 
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AOGCC prohibits the waste of oil and natural gas in accordance with the Alaska Oil and Gas Conservation Act 
(Section 31.05.170 (15)(H)). The Act specifies that the release, burning, or escape of oil or natural gas from an oil 
or gas producing well is prohibited unless authorized by AOGCC (USDOE 2019). Any wasted oil or natural gas 
must also be reported to AOGCC with a statement of compliance actions (USDOE 2019). The State of Alaska 
also prohibits flaring except in the case of emergencies or system testing (20 AAC 25.235). This regulation 
authorizes flaring under several conditions, including for periods less than one hour if resulting from emergencies, 
operational upsets, or planned lease operation. For flaring longer than 1 hour, AOGCC would consider 
authorization if flaring was necessary for safety in emergencies, in which case operators must report the volume 
of gas flared. In addition, if the Willow Processing Facility is subject to “major” source permitting requirements, 
any flares planned to be constructed at the facility would be subject to best available control technology 
requirements to minimize emissions from flares, as well as any other applicable equipment (SLR 2022). 


ADEC regulations applicable to flaring are included in Standard Permit Condition LX — Visible Emissions and 
Particulate Matter Monitoring Plan for Liquid Fuel-Burning Equipment and Flares, revised on July 22, 2020, 
which is adopted by reference in 18 Alaska Administrative Code 50.346. 


BLM also has flaring provisions in-the Notice to Lessees and Operators of Onshore Federal and Indian Oil and 
Gas Leases, Royalty or Compensation for Oil and Gas Lost, commonly referred to as NTL-4 (44 Federal Register 
76600 [1979]), that are applicable to operators of federal oil and gas leases. Currently, the provision requires 
payment of royalties for oil or gas that is flared without authorization or if it is determined to be “avoidably lost.” 


1.1.2 Characterization of Existing Air Quality in the Analysis Area 


Regional air quality is affected by a variety of factors, including climate, meteorology, and the magnitude and 
location of air pollutant sources. This section provides descriptions of the regional climate, meteorology, and 
existing regional sources of air pollution that affect air quality in the analysis area. Existing air quality in the 
analysis area is assessed through a review of recent ambient air quality monitoring data and AQRVs. 


1.1.2.1 Climate and Meteorology 

The Project is located on the North Slope within the National Petroleum Reserve in Alaska (NPR-A). Several 
monitoring stations were used to characterize climate and meteorology in the analysis area. Monthly average 
precipitation and temperature data were acquired from the National Oceanic and Atmospheric Administration 
National Weather Service (NWS) stations at Umiat, Kuparuk, Utqiagvik (Barrow), and Nuigsut (Figure E.3.2). A 
monitoring station operated by ConocoPhillips Alaska, Inc. (CPAI) at Nuiqsut was used to characterize prevailing 
wind patterns. 


Table E.3.3 provides summaries of the average monthly temperature and precipitation from the NWS stations 
shown in Figure E.3.2. The annual average temperature in the NPR-A is approximately 10 degrees Fahrenheit (F), 
with monthly average maximum temperatures below freezing from October to May (BLM 2012). The coldest 
temperatures (usually in February) range from -10 degrees to -15 degrees F at the maximum and -25 degrees to - 
30 degrees F at minimum on average (see Table E.3.3). Summer temperatures rise above freezing, with the 
highest temperatures typically occurring in July. The average maximum and minimum temperatures in July range 
from 45 degrees F to 65 degrees F and 35 degrees F to 40 degrees F, respectively. 


Precipitation in the analysis area is low, with Nuiqsut receiving 2.74 inches of precipitation per year on average 
(see Table E.3.3). Precipitation is highest during summer, with over three-fourths of the total annual precipitation 
falling between June and September. Although snowfall is sparser during the summer months, it can occur during 
any month; the highest average snowfall rates occur in October. Snow is generally on the ground from October to 
May (BLM 2012). 


The wind rose in Figure E.3.3 shows the distribution of wind direction and speeds measured at the CPAI Nuiqsut 
monitoring station, located approximately 46 km (28.5 miles) east-northeast of the Project, from 2016 to 2020. 
The prevailing wind direction at Nuiqsut was from the northeast with wind speeds averaging 4.9 meters per 
second (m/s) (11.0 miles per hour). The maximum observed wind speed was 22.4 m/s (50.1 miles per hour) and 
calm winds were infrequent, occurring for less than 1.5 % of hours during the 5-year period. Figures E.3.4 
through E.3.7 provide seasonal wind patterns for the winter, spring, summer, and fall seasons, respectively, for the 
5-year period. 
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Figure E.3.2. Monitoring Stations Used to Characterize Climate and Meteorology in the Project Area 
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Table E.3.3. Monthly Climate Summary Data at Monitoring Stations in the Air Quality Analysis Area 
Utqiagvik (Barrow)* Jan Feb Mar Apr May Jun Jul Aug Sept Oct 


Dec Annual 


51.6 40. 


Sis 
S 

= 

ic} 

+ 


-22.9 41.6 |°387 Up se 
Average Total Precipitation (in 0.10 0.19 1.04 | 1.04 | 0.40 | 0.04 2.74 


Note: F (Fahrenheit); in (inches); Max. (maximum); Min. (minimum). The sum of the monthly precipitation totals may not equal the annual total because of different data completeness requirements for monthly 
and annual data. 

* Source: National Oceanic and Atmospheric Administration (NOAA) National Weather Service (NWS) data, obtained from the Western Regional Climate Center(https://wrcec.dri.edu/summary/Climsmak.html). 
Period of record: Utgiagvik (1901 to 2016); Umiat (1945 to 2001); Kuparuk (1983 to 2016). Historical records are under Utqiagvik’s former name of Barrow. 

> Units of total precipitation are inches of liquid water equivalent. 

© Source: NOAA NWS data obtained from NOAA National Centers for Environmental Information (https://www.ncdc.noaa.gov/cdo-web/datatools/normals). Period of record: 1981 to 2010. As of January 6, 2022, 
the 1981-2010 period is the most recent climate normal (i.e., 3 decades) available. 

4 Source: NOAA NWS data obtained from Natural Resources Conservation Service (http://agacis.rec-acis.org/?fips=02185). Period of record: 1998 to 2021. Months within each year with > 1 missing day are 
omitted from averages. Annual data with > 1 missing day is also omitted from averages. Due to this, the sum of monthly averages does not equal the annual average. The annual value is based on 2002, 2004. 2009, 
and 2011 years only, since only those years satisfied the data completeness criteria. 
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Figure E.3.3. Wind Rose Data from the ConocoPhillips Alaska, Inc. Nuiqsut Monitoring Station for the 
Period 2016 to 2020* 
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Figure E.3.4. Wind Rose Data from the ConocoPhillips Alaska, Inc. Nuiqsut Monitoring Station for the 
Winter Months (December, January, and February) during 2016 to 2020* 
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Records: 13248 
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Figure E.3.5. Wind Rose Data from the ConocoPhillips Alaska, Inc. Nuiqsut Monitoring Station for the 
Spring Months (March, April, and May) during 2016 to 2020* 
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Figure E.3.6. Wind Rose Data from the ConocoPhillips Alaska, Inc. Nuiqsut Monitoring Station for the 
Summer Months (June, July, and August) during 2016 to 2020* 
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Figure E.3.7. Wind Rose Data from the ConocoPhillips Alaska, Inc. Nuiqsut Monitoring Station for the 
Fall Months (September, October, and November) during 2016 to 2020* 


1.1.2.2 Existing Regional Sources of Air Pollution 

A summary of existing regional emissions for the North Slope and adjacent waters (Beaufort Sea and Chukchi 
Sea Planning Areas) is available from the 2012 baseline scenario of the Bureau of Ocean Energy Management 
Arctic Air Quality Modeling Study: Emissions Inventory, Final Task Report (Fields Simms, Billings et al. 2014). 
Existing emissions from onshore sources (e.g., oil and gas production and exploration, airports, pipelines, non-oil- 
and gas-related stationary and mobile sources) comprise the majority of the total existing emissions, and 
emissions from offshore sources (e.g., drilling rigs, survey/drilling vessels and aircraft, commercial vessels) are 
small in comparison (Fields Simms, Billings et al. 2014). Overall, onshore oil and gas sources comprise the 
largest fraction of existing emissions for all CAPs except particulate matter less than or equal to 10 microns in 
aerodynamic diameter (PMjo) and particulate matter less than 2.5 microns in aerodynamic diameter (PM2s) for 
which dust from unpaved roads comprises the largest fraction (Fields Simms, Billings et al. 2014). The major 
existing sources of HAPs in the region are onshore oil and gas, other nonroad vehicles and equipment, on-road 
vehicles, and waste incineration, landfills, and other combustion sources. 
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1.1.3 Air Quality Monitoring 


1.1.3.1 Criteria Air Pollutants* 

CPAI operates the Nuigsut Monitoring Station, which is the most representative station in the region of the 
Project (see Figure E.3.1) (BLM 2018). Monitoring data from the CPAI Nuiqsut monitoring station are provided 
in Table E.3.4 for 2018 through 2020. All CAPs are monitored except for lead, for which there are no monitoring 
stations in the analysis area. All of the monitored concentrations are well below the NAAQS and AAAQS. This is 
consistent with the existing air quality of the larger analysis area, which is designated as 
“attainment/unclassifiable” for all CAPs. 


Table E.3.4. Measured Criteria Air Pollutant Concentrations at the Nuic g 
Averaging Rank 2018 2019 2020 Avg NAAQS/ _ Below 
Period AAAQS' NAAQS/ 
3 | AAAQS? 


3 
2. 5 
gm’) [Annual _ [Average 
: [Shours _|4th highest daily max 


Note: AAAQS (Alaska Ambient Air Quality Standards); Avg. (average); CO (carbon monoxide); max (maximum); NAAQS (National Ambient Air Quality 
Standards); NO> (nitrogen oxides); O3; (ozone); PM jo (particulate matter less than or equal to 10 microns in aerodynamic diameter); PM2 s (particulate matter 
less than 2.5 microns in aerodynamic diameter); ppb (parts per billion); ppm (parts per million); SO (sulfur dioxide); g/m* (micrograms per cubic meter). 
NAAQS/AAAQS for ozone (O;) were converted from ppm to ppb and sulfur dioxide (SO>) 24-hour and annual standards were converted from g/m’ to ppb. 
Data used in the table has not been reviewed by the Alaska Department of Environmental Conservation for Prevention of Significant Deterioration quality; 
however, the selection of the Nuiqsut station for monitoring data was made during the development of the Willow Environmental Impact Statement 
modeling protocol, which was reviewed by air specialists at the Alaska Department of Environmental Conservation and other agencies. 


Ol 
2 |S 
o 


1.1.3.2 Visibility* 


_ Visibility and air pollutant concentration data is collected by Interagency Monitoring of Protected Visual 

_ Environments at monitoring sites close to Class I areas across the country. The three closest monitors to the 

_ Project with available data are Toolik Lake Field Station, Gates of the Arctic National Park and Preserve (a Class 
_ IL area), and Denali National Park (a Class I area) (see Figure E.3.1). Data from these monitors are presented in 

_ Figures E.3.8 through E.3.13 and Table E.3.5. Denali National Park is outside the analysis area for air quality but 
_ is included here as it is the closest Class I area. Denali National Park has the longest visibility data record from 
1989 through 2019. Gates of the Arctic National Park has available visibility data from 2010 through 2014, and 

_ Toolik Lake Field Station only has data for 2019 because it is a new Interagency Monitoring of Protected Visual 
_ Environments (IMPROVE) site that became operational in November 2018. Data is shown for the 20% haziest 
_and 20% clearest days. The 20% haziest days include anthropogenic and natural influences following the 

_ algorithm of EPA (2003) as revised by IMPROVE in December 2019 and is influenced by natural emission 


_ sources such as wildfires. At Gates of the Arctic, the haze index on the haziest days shows a consistent downward 


_ trend (through the years of the plot available from IMPROVE) that is near estimated natural visibility conditions! 
_ of 7.7 dv (visual range of approximately 129 miles), while the haze index on the clearest days has consistently 

_ been between 3 and 4 dv, which is slightly above the estimated natural conditions of 2.8 dv (visual range of 

| approximately 349 km [217 miles]). At Denali National Park, the haze index shows generally decreasing trends 

_ for both the haziest days and the clearest days, but the haziest days have some outlier years, most notably 2004, 

_ likely due to wildfires. Estimated natural visibility conditions' at Denali National Park are 7.3 dv (visual range of 


_ approximately 209 km [130 miles]) and 1.8 dv (visual range of approximately 360 km [224 miles]) for the haziest 


! http://vista.cira.colostate.edu/IMPROVE/Data/NaturalConditions/nc2 12 2019 2 
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| and clearest days, respectively. In recent years, the haze index values approach those estimated for natural 
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Figure E.3.8. Visibility Data for Gates of the Arctic National Park 
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Source: FED 2020 
Figure E.3.9. Visibility on the Haziest Days for Gates of the Arctic National Park 
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Figure E.3.10. Visibility on the Clearest Days for Gates of the Arctic National Park 
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Source: https://views.cira.colostate.edu/fed/SiteBrowser/Default.aspx?appkey=SBCF VisSum 
Figure E.3.11. Visibility Data for Denali National Park* 
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Figure E.3.12. Visibility on the Haziest Days for Denali National Park* 
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Figure E.3.13. Visibility on the Clearest Days for Denali National Park 


Note: AK (Alaska); dv (deciview); IMPROVE (Interagency Monitoring of Protected Visual Environments) 
Source: https://views.cira.colostate.edu/fed/SiteBrowser/Default.aspx?appkey=SBCF_VisSum 
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| monitored wet deposition fluxes of ammonium (NH), NO3, and SO,” at each site are provided in Figures 
_E.3.14, E.3.15, and E.3.16, respectively. The blue dots on the graphs indicate yearly concentrations that have met 
_ the annual completeness criteria, while the red dots indicate that yearly concentrations have not met the annual 

_ completeness criteria. Trendlines are also shown in black and represent a 3-year moving average where the 

_ minimum data completeness criteria are met for that 3-year period. The wet deposition fluxes of NHs, NOs’, and 
_ SO,* are small at all monitors (most annual values below 1.0 kilogram per hectare per year) with no apparent 


_ trend in most cases. However, the wet deposition fluxes of NOx at Poker Creek have shown an upward trend over 


the last decade, and 2019 and 2020 had the two highest measurements in over two decades. 


_ The NADP also provides estimates of total (wet and dry) sulfur and nitrogen deposition for critical load analysis 

_and other ecological studies using a hybrid approach with modeled and monitoring data (NADP 2014). Wet 

_ deposition data from NTN, along with air concentration data from networks such as the Clean Air Status and 

_ Trends Network (CASTNET), is used (EPA 2018a). The estimated total deposition flux of nitrogen and sulfur is 
_ provided in Figure E.3.17 for Denali National Park for 1999 through 2020, which is the only monitor in Alaska 

_ with recent CASTNET data (DEN417 in Figure E.3.1). The highest monitored total deposition fluxes of nitrogen 
_and sulfur occurred in 2002 and were 0.741 kilograms of nitrogen per hectare per year (kg N/ha/year) and 0.601 

_ kilograms sulfur per hectare per year (kg S/ha/year), respectively. The mean deposition fluxes of nitrogen and 

_ sulfur are 0.297 kg N/ha/year and 0.287 kg S/ha/year, respectively. The total deposition flux of nitrogen was well 
_ below the critical load for nitrogen deposition defined by the FLMs for the tundra ecoregion of Alaska (1.0 to 3.0 
_kg N/ha/year) in all years. 
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Figure E.3.14. Trends in Wet Deposition of Ammonium (NH;) at Poker Creek (NTN Site AK01), Denali 
National Park (NTN Site AK03), Gates of the Arctic National Park (NTN Site AK06), and 
Toolik Lake Field Station (NTN Site AK96)* 
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igure E.3.15. Trends in Wet Deposition of Nitrates (NO;) at Poker Creek (NTN Site AK01), Denali 
National Park (NTN Site AK03), Gates of the Arctic National Park (NTN Site AK06), and 
Toolik Lake Field Station (NTN Site AK96)* 
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Source: https://nadp.slh.wisc.edu/sites/ntn-AK03/ 

Figure E.3.16. Trends in Wet Deposition of Sulfates (SO,”) at Poker Creek (NTN Site AK01), Denali 
National Park (NTN Site AK03), Gates of the Arctic National Park (NTN Site AK06), and 
Toolik Lake Field Station (NTN Site AK96)* 
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Source: EPA 2018a 
Figure E.3.17. Total Nitrogen and Sulfur Deposition Flux at Denali National Park* 
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